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AVERTISSEMENT. 



J'ai réuni dans cet ouvrage les divers écrits publiés 
par Gauss sur la Méthode des moindres carrés. 

L'illustre géomètre, que les sciences viennent de 
perdre, attachait une grande importance à cette partie 
de ses travaux, et la meilleure manière de combiner 
les observations était, à ses yeux, tin des problèmes les 
plus importants de la philosophie naturelle. 

Gauss n'ignorait pas les critiques dont sa théorie 
a été l'objet; mais son opinion bien arrêtée était que 
les géomètres adopteraient entièrement ses idées lors- 
que ses Mémoires, aujourd'hui fort rares, se seraient 
répandus davantage. C'est pour cela, sans doute, 
qu'il a bien voulu m'écrire qu'il me verrait avec le 
plus grand plaisir en publier la traduction. Ces 
Mémoires forment un Traité complet de la combinai- 
son des observations , qui n'exige ni commentaires ni 
annotations. Les questions de priorité sur lesquelles 
se sont engagées des discussions assez vives, y sont 
traitées brièvement, mais de la manière la plus nette 
et la plus loyale. J'ai donc dû me borner au rôle de 
traducteur : c'était le seul qui fût utile, et le seul 
d'ailleurs que Gauss m'eût autorisé à prendre. 

.1. BERTRAND. 
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1. 

Quelque soin que l'on apporte aux observations qui con- 
cernent la mesure des grandeurs physiques , elles sont for- 
cément soumises à des erreurs plus ou moins considérables. 
Ces erreurs, dans le plus grand nombre des cas, ne sont 
pas simples, mais découlent à la fois de plusieurs sources 
distinctes qu'il est bon de distinguer en deux classes. 

Certaines causes d'erreurs dépendent, pour chaque obser- 
vation, de circonstances variables et indépendantes du ré- 
sultat que Ton obtient : les erreurs qui en proviennent 
sont nommées irrégulières ou fortuites, et de même que les 
circonstances qui les produisent, leur valeur n'est pas sus- 
ceptible d'être soumise au calcul. Telles sont les erreurs 
qui naissent de l'imperfection de nos organes et toutes celles 
qui sont dues à des causes extérieures irrégulières , comme , 
par exemple, les trépidations de l'air qui rendent la vision 
moins nette 5 quelques-unes des erreurs dues à Timperfec- 
o. 1 
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lion inévitable des meilleurs instrumenta appartiennent à 
la même catégorie. Nous citerons, par exemple, la rugosité 
de la partie intérieure du niveau, le défaut do rigidité 
absolue, etc. 

Il existe, au contraire, d'autres causes qui, dans toutes 
les observations de même nature, produisent une erreur 
identique, ou dépendant de circonstances essentiellement 
liées au résultat de l'observation. Nous appellerons les er- 
reurs de cette catégorie, des erreurs constantes ou régulières. 

Il est du reste évident que cette distinction est jusqu'à un 
certain point relative et dépend du sens plus ou moins 
large que Ton veut attacher à l'idée d'observations de même 
nature. Par exemple, si l'on répète indéfiniment la mesure 
d'un même angle, les erreurs provenant d'une division 
imparfaite du limbe appartiendront à la classe des erreurs 
constantes. Si, au contraire, on mesure successivement 
plusieurs angles différents, les erreurs dues à l'imperfec- 
tion de la division seront regardées comme fortuites tant 
que l'on n'aura pas formé la table des erreurs relatives à 
chaque division. 

2. 

■ 

Nous excluons de nos recherches la considération des 
erreurs régulières. C'est à l'observateur qu'il appartient 
de rechercher avec soin les causes qui peuvent produire 
une erreur constante , pour les écarter s'il est possible, ou 
tout au moins apprécier leur effet, afin de le corriger sur 
chaque observation, qui donnera alors le môme résultat 
que si la cause constante n'avait pas existé. Il en est tout 
autrement des erreurs irrégulières : celles-là, par leur na- 
ture, se refusent à tout calcul, et il faut bien les tolérer dans 
les observations. On peut cependant, par une combinaison 
habile des résultats, réduire autant que possible leur in- 
fluence. C'est à cette question importante que sont consa- 
crées les recherches suivantes. 
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3. 

Les erreurs qui , dans des observations de même espèce, 
proviennent d'une cause simple et déterminée se trouvent 
renfermées entre certaines limites que Ton pourrait sans 
aucun doute assigner, si la nature de cette cause était 
elle-même parfaitement connue. Dans la plupart des cas, 
toutes les erreurs comprises entre ces limites extrêmes 
doivent être regardées comme possibles. Une connais- 
sance approfondie de chaque cause apprendrait si toutes 
ces erreurs ont une facilité égale ou inégale, et, dans le 
second cas, quelle est la probabilité relative de chacune 
d'elles. La même remarque s'applique à Terreur totale qui 
provient de la réunion de plusieurs erreurs simples. Cette 
erreur sera, elle aussi, renfermée entre deux limites dont 
Tune sera la somme des limites supérieures, l'autre celle 
des limites inférieures, correspondant aux erreurs simples. 
Toutes les erreurs comprises entre ces limites seront possi- 
bles, et chacune pourra résulter, d'une infinité de manières, 
de valeurs convenables attribuées aux erreurs partielles. On 
comprend néanmoins , en écartant les difficultés purement 
analytiques , qu'il y a possibilité d'apprécier la probabilité 
plus ou moins grande de chaque résultat, si Ton suppose con- 
nues les probabilités relatives à chacune des causes simples. 

Certaines causes pourtant produisent des erreurs qui ne 
peuvent pas varier suivant une loi continue, mais qui, au 
contraire , sont susceptibles d'un nombre fini de valeurs : 
nous pouvons citer, comme exemple, les erreurs qui pro- 
viennent de la division imparfaite des instruments (si tou- 
tefois on veut les classer parmi les erreurs fortuites), car 
le nombre des divisions, dans un instrument donné, est 
essentiellement fini. Il est clair néanmoins que, si toutes les 
causes qui concourent à produire l'erreur totale ne sont pas 
supposées dans ce cas , leur somme formera une série sou- 
mise à la loi de continuité, ou, tout au moins, plusieurs 

i . 
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séries distinctes , s'il arrive qu'en plaçant par ordre de gran- 
deur toutes les valeurs possibles des erreurs discontinues, 
la différence entre deux ternies consécutifs de la série soit 
moindre que la différence entre les limites extrêmes des 
erreurs soumises à la loi de continuité. Dans la pratique, un 
pareil cas ne se présentera presque jamais 5 il supposerait 
des défauts trop grossiers dans la construction de l'instru- 
ment. 

4. 

Désignons par la notation <p (x) la facilité relative d'une 
erreur x : on doit entendre par là , à cause de la conti- 
nuité des erreurs, que <p (x) dx est la probabilité que Ter- 
reur soit comprise entre les limites x etx •+- dx. Il n'est pas 
possible, en général, d'assigner la forme de la fonction©, 
et Ton peut même affirmer que cette fonction ne sera jamais 
connue dans la pratique. On peut néanmoins établir plu- 
sieurs caractères généraux qu'elle doit nécessairement pré- 
senter : <?(x) est évidemment une fonction discontinue; 
elle s'annule pour toutes les valeurs de x non comprises 
entre les erreurs extrêmes. Pour toute valeur comprise entre 
ces limites, la fonction est positive (en excluant le cas in- 
diqué à la fin du paragraphe précédent); dans la plupart 
des cas, les erreurs égales et de signes contraires seront 
également probables, et l'on aura 

j(x) =<p( — x). 
Enfin, comme les petites erreurs sont plus facilement com- 
mises que les grandes, <p (x) sera en général maximum pour 
x = o et diminuera sans cesse lorsque x croîtra. 
L'intégrale 

J <j> (x) dx 

exprime la probabilité pour que Terreur, encore inconnue, 
tombe entre les limites a et b. On en conclut que la valeur 
de cette intégrale prise entre les limites extrêmes des erreurs 
possibles sera toujours égale à l'unité. Et comme © (x) est 
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nulle pour les valeurs non comprises enlre ces limites , on 
peut dire, dans tous les cas, que 



L 



oc 

y(x)dx = I 



5. 

Considérons l'intégrale 



X 



xo(x) dx 

et représentons sa valeur par À\ Si les causes d'erreur sont 
telles, qu'il n'y ait aucune raison pour que deux erreurs 
égales et de signes contraires aient des facilités inégales, on 
aura 

?(*) = <?{— *)> 

et, par suite, 

k =r o. 

Nous en conclurons que, si k ne s'évanouit pas et a, pat- 
exemple, une valeur positive, il existe nécessairement une 
cause d'erreur qui produit uniquement des erreurs positives 
ou qui , tout au moins , les produit plus facilement que les 
erreurs négatives. Cette quantité k, qui est la moyenne de 
toutes les erreurs possibles , ou encore la valeur moyenne 
de x, peut être désignée commodément sous le nom de partie 
constante de V erreur. Du reste , on prouve facilement que 
la partie constante de Terreur totale est la somme des par- 
ties constantes des erreurs simples qui la composent. 

Si la quantité A est supposée connue et qu'on la retranche 
du résultat de chaque observation , en désignant par x' l'er- 
reur de l'observation ainsi corrigée, et la probabilité cor- 
respondante par <p' (x 1 ), on aura 

— ?'(*') = ?(•*) 

et, par suite, 

I x'y'[x') dx'=. l xy[x)dx — / I y (x)dx = k — k = o ; 
J—oo J—so J— x> 



- 
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en sorte que les erreurs des observations corrigées n'ont pas 
de partie constante. Ce qui, du reste, semble évident à 
priori. 

6. 

La valeur de l'intégrale 

x«p (x) dx , 

c'est-à-dire la valeur moyenne de fait connaître l'exis- 
tence ou la non-existence d'une erreur constante, ainsi que 
la valeur de cette erreur; de même l'intégrale 




c'est-à-dire la valeur moyenne de a:*, paraît très-propre à 
définir et à mesurer, d'une manière générale, l'incertitude 
d'un système d'observations ; de telle sorte qu'entre deux 
systèmes d'observations inégalement précises, on devra re- 
garder comme préférable celui qui donne à l'intégrale 

x 1 <$> (x) dx 

- 00 

une moindre valeur. Si Ton objecte que cette convention 
est arbitraire et ne semble pas nécessaire, nous en conve- 
nons volontiers. La question qui nous occupe a, dans sa 
nature môme, quelque chose de vague et ne peut être bien 
précisée que par un principe jusqu'à un certain point arbi- 
traire. La détermination d'une grandeur par l'observation 
peut se comparer, avec quelque justesse, à un jeu dans lequel 
il y aurait une perte à craindre et aucun gain à espérer : 
chaque erreur commise étant assimilée à une perte que Ton 
fait, la crainte relative à un pareil jeu doit s'exprimer par la 
perte probable, c'est-à-dire par la somme des produits des 
diverses pertes possibles par leurs probabilités respectives. 
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Mais quelle perte doit-on assimilera une erreur détermi- 
née? C'est ce qui n'est pas clair en soi 5 cette détermination 
dépend en partie de notre volonté. Il est évident, d'abord, que 
la perte ne doit pas être regardée comme proportionnelle à 
Terreur commise 5 car, dans cette hypothèse, une erreur 
positive représentant une perte, l'erreur négative devrait 
être regardée comme un gain : la grandeur de la perte doit, 
au contraire, s'évaluer par une fonction de Terreur dont la 
valeur soit toujours positive. Parmi le nombre infini de 
fonctions qui remplissent cette condition, il semble natu- 
rel de choisir la plus simple, qui est, sans contredit, le 
carré de l'erreur, et, de celte manière, nous sommes con- 
duit au principe proposé plus haut. 

Laplace a considéré la question d'une manière analogue, 
mais en adoptant, pour mesure de la perte , l'erreur elle- 
même prise positivement. Cette hypothèse, si nous ne nous 
faisons pas illusion , n'est pas moins arbitraire que la nôtre : 
faut-il, en effet, regarder une erreur double comme plus 
ou moins regrettable qu'une erreur simple répétée deux 
fois, et faut-il, par suite, lui assigner une importance 
double ou plus que double? C'est une question qui n'est 
pas claire, et sur laquelle les arguments mathématiques 
n ont aucune prise; chacun doit la résoudre à son gré. 
On ne peut nier pourtant que l'hypothèse de Laplace ne 
s'écarte de la loi de continuité et ne soit, par conséquent, 
moins propre à une étude analytique; la nôtre, au contraire, 
se recommande par la généralité et la simplicité de ses con- 
séquences. 

Posons, en conservant les notations précédentes, 

y (jt) jc* dx =r m* : 

nous appellerons m Y erreur moyenne à craindre ou, sim- 
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plement, V erreur moyenne des observations considérées. 
Nous ne limitons pas, du reste, cette dénomination au ré- 
sultat immédiat des observations , mais nous l 'étendons , au 
contraire, à toute grandeur qui peut s'en déduire d'une 
manière quelconque. Il faut bien se garder de confondre 
cette erreur moyenne avec la moyenne arithmétique des 
erreurs , dont il est question dans l'art. 5. 

Si nous comparons plusieurs systèmes d'observations ou 
plusieurs grandeurs résultant d'observations auxquelles on 
n'accorde pas la même précision, nous regarderons leur 
poids relatif comme inversement proportionnel à m 1 , et leur 
précision comme inversement proportionnelle à m. Afin de 
pouvoir représenter les poids par des nombres, on devra 
prendre, pour unité, le poids d'un certain système d'obser- 
vations arbitrairement choisi. 

8. 

Si les erreurs des observations ont une partie constante, 
en la retranchant de chaque résultat obtenu , l'erreur 
moyenne diminue, le poids et la précision augmentent. En 
conservant les notations de l'art. 5, et désignant par m' 
l'erreur moyenne des observations corrigées , on aura en 
effet 

m' 3 = f * x'^'{x')dx f = f °° {x — &)*it{x)dx 

j — » «/_ » 

X r x> /■» ce s\ oc 

x 7 ? (x) dx — 2 k I jç(x) dx -f- h 1 I y (x) dx 
•OC J — X J — oo 

— m - — 2 + V — m- — f.\ 

Si , au lieu de retrancher de chaque observation le nom- 
bre Ar, on retranchait un autre nombre /, le carré de 
Terreur moyenne deviendrait 
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9. 

Soient X un coefficient déterminé et /x la valeur de l'in- 
tégrale 

X m 
X m 



f 



{jl sera la probabilité que Terreur d'une certaine observa- 
tion soit moindre que X m en valeur absolue; i — fx sera, 
au contraire, la probabilité que cette erreur surpasse X m. 

Si, pour p = -> Xm a la valeur p, il y aura probabilités 

égales pour que l'erreur soit plus petite ou plus grande 
que p : p pourra donc être appelé Y erreur probable. La re- 
lation qui existe entre X et p. dépend de la nature de la 
fonction q>, qui est inconnue dans la plupart des cas. Il est 
intéressant d'étudier cette relation dans quelques cas parti- 
culiers. 

I. Si les limites extrêmes des erreurs possibles sont 
-ha et — a, et si, entre ces limites, toutes les erreurs 
sont également probables, la fonction <p (x) sera constante 

entre ces mêmes limites , et, par conséquent, égale à 

Par suite, on aura 

tant que X sera inférieur ou égal à enfin, 



m 



= m — 0 r 8660254 rn y 



et la probabilité pour que Terreur ne surpasse pas Terreur 
moyenne est 

y/^ = 0,5773503. 

II. Si — a et H- a sont encore les limites des erreurs 
possibles , si Ton suppose , de plus , que la probabilité de 
ces mêmes erreurs aille en décroissant à partir de Terreur o 



"s 
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comme les termes d'une progression arithmétique, on aura 

ci - x 

<y (x) — : 

pour les valeurs de x comprises entre o et -h a , et 

a -f- r 

t(*) = — 

pour les valeurs comprises entre o et — a : de là on déduit 



m 



tant que X est compris entre o et \ ë ; 

X = — — °V> 
tant que a est compris entre o et i; et, enfin, 

p = //i (yt) — V3) = 0,7 17438g m. 

Dans ce cas, la probabilité que Terreur restera inférieure 
à Terreur moyenne sera 

y/|_£ = 0,6498299. 
III. Si nous supposons la fonction <j> (x) proportionnelle 



ne A * [ce qui, en réalité, n'est vrai qu approximative- 
ment (*)], elle devra être égale à 



on en conclut 



(voir Disqutsitiones générales circa seriem injitiitam, 
art. 28, Mémoires de Gotfingue, tome II). 

( * ) Il faul se reporter, pour comprendre cette remarque, à un chapitre 
du Theoria Motus Corporum cœlestium , dans lequel M. Gauss montre que 
celte loi de probabilités est la plus vraisemblable que l'on puisse adopter. 
A la fin du volume nous reproduisons ce chapitre, dans lequel l'illustre au- 
teur a fait connaître pour la première fuis la méthode des moindres carrés. 

J. B. 
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Si Ton désigne par &z la valeur de l'intégrale 

— I <> <h, 

Le tableau suivant donne quelques valeurs de cette quan- 



on aura 



tite 



0,6744897 
o, 84i6su 3 
1 ,0000000 
1 ,o364334 
1 ,2815517 
1,6448537 
2,5758293 
3,2 9 i83oi 
3,8905940 



o,5 
0,6 

0,6826895 

°>7 
0,8 

o»9 
°>99 
°>999 
°>9999 



10. 

Quoique la relation qui lie X à p dépende de la nature 
de la fonction <p, on peut cependant établir quelques résul- 
tats généraux, qui s'appliquent à tous les cas dans lesquels 
cette fonction ne sera pas croissante avec la valeur absolue 
de la variable x \ alors on aura les théorèmes suivants : 

X ne dépassera pas [â v3 toutes les fois que sera infé- 

» r 

, 2 

rieur a ^: 

X ne dépassera pas 



3^1- 



, toutes les fois que jtx sur- 



passera — 
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Lorsque p= ^ les deux limites coïncident et / ne peut 

pas être supérieur à 

Pour démontrer ce théorème remarquable, représentons 
parj- la valeur de l'intégrale 

alors sera la probabilité pour qu'une erreur soit com- 
prise entre — x et -hx. 
Posons 

* = +(r)> rf+(/) = f(r)<r, rf+'(7) = r(j) rf /. 

on aura 

+ (o) = o etf (»= 1 ; 

H*) + f (— *) 

on en conclut, en ayant égard aux hypothèses qui ont été 
faites, que, depuis y=o jusqu'à y=\, (y) est tou- 
jours croissant, ou du moins n'est pas décroissant, et que, 
par suite, Y (y) est toujours positif, ou du moins n'est 
pas négatif. Or nous avons 

d.y V (y) = dy -4- y f ' {y) dy , 

par suite, 

r+'(r)- + (r) = f'rf'W^; 

(^) — (j^) a donc une valeur constamment positive, 
ou du moins cette expression ne sera jamais négative. 

Il suit de là que i — y y (y) Sera tou j ours P os * 1 ^ el moin- 
dre que l'unité. Soit f la valeur de cette différence pour 
y = fx ; à cause de <Jj (p) = X/rc, on a 

r— ^ m 
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d'où Ton conclut 

Cela posé, considérons la fonction 

"km - v 

que nous désignerons par F (jr), et posons 

m 

on aura évidemment 

F( f t) = Xnf = +(p), 

Or, puisque <J/ (y) croît continuellement (ou du moins ne 
décroît pas, car c'est ainsi qu'on doit toujours l'entendre) 
lorsque y croît, et que, d'un autre côté, F' (y) est con- 
stant, la différence 

sera positive pour toutes les valeurs de y plus grandes que 
et négative pour les valeurs de y plus petites que p. On en 
conclut que la différence ^ (y) — F (y) est toujours posi- 
tive, et, par suite, (y) sera certainement plus grand que 
F (y) en valeur absolue, tant que la fonction F (y) sera 
positive, c'est-à-dire depuis y = pf jusqu'à y — i . La 
valeur de l'intégrale 

sera donc inférieure à celle de l'intégrale 

f [+(r)]**s 

cl à fortiori moindre que 



fui 



I 
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c'est-à-dire moindre que m 5 . Or la première de ces inté- 
grales a pour valeur 

"3u<(i-/;- " 

donc 

^ (<-(*/)• 

désignant, on ne Ta pas oublié, un nombre compris 
entre o et 1 . 

Si nous considérons f comme variable, la fraction 
aura pour diirérenticlle 

cette fraction sera donc continuellement décroissante lors- 

que y croîtra de o à i, et que Ton aura, en outre, p<C g* 

sa valeur maximum correspondra h f = o et sera égale à 
3/ut% de sorte que, dans ce cas, le coefficient X sera cer- 
tainement inférieur, ou du moins ne sera pas supérieur 

à ^ y^. Ce qu'il fallait démontrer. 

Lorsqu'au contraire jtx est plus grand que ^' la valeur 
de la fonction sera maximum lorsque 

2 — 3 p -+- itf = o , 

c'est-à-dire lorsque 

v- 

et cette valeur maximum sera 

4 . 
9('-c)' 
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par conséquent, dans ce cas, le coefficient X n'est pas plus 

grand que ■- ■,-- > comme nous lavions annoncé. 

3^1— p 

Faisons, par exemple, 

i 

P = 2 ^ 

alors X ne peut pas surpasser c'est-à-dire que l'er- 

reur probable ne peut pas surpasser o,866o2f>4 m , à la- 
quelle elle devient égale dans le premier cas examiné 
(art. 9) : on conclut facilement de notre théorème que p 

n'est pas moindre que ^y/^' lanl *l ue ^ est moindre 
que et qu'au contraire, il n'est pas inférieur à 

î y 2 , lorsque X est plus grand que i/ 1- 

L'intégrale 

y (.r) r/.r , 

se présentant dans plusieurs problèmes que nous aurons à 
traiter, il ne sera pas inutile de l'évaluer dans quelques cas 
particuliers. 
Posons 



! x* y (x) dx = n k . 

«' — 00 

I. Lorsque 



«p (x) = — , 



pour les valeurs de a: comprises entre — a et -ha, on a* 

«« 9 

5 5 
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II. Lorsque 

ce (x 1 = — -2— - 

' k ' « 

(Il p cas, art. 9), or étant encore compris entre — « et -ha, 
on a 

/?*=—= <r = — m*. 
io 5 

III. Dans le troisième cas, lorsque 

X* 

on trouvera, d'après les résultats obtenus dans le Mémoire 
cité plus haut , 

4 

On peut d'ailleurs démontrer qu'en restant dans les hypo- 
thèses admises au paragraphe précédent, le rapport ~ 

n'est jamais inférieur à 2» 

12. 

Désignons par a:, x 1 , x" , etc., les erreurs commises 
dans des observations de même espèce, et supposons que 
ces erreurs soient indépendantes les unes des autres. Soit , 
comme plus haut, y(x) la probabilité relative de l'erreur 
x \ considérons une fonction rationnelle des variables 

f ff 

L'intégrale multiple 

(ï) j < f (x)< 9 {x')y[x")...dxdx'dx".,., 

étendue à toutes les valeurs des variables .r, x* , x", etc., 
pour lesquelles la valeur de y tombe entre les limites don- 
nées o, yj. représente la probabilité que la valeur de y soit 
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comprise entre o et yj. Or cette intégrale est évidem- 
ment une fonction de yj. Représentons sa différentielle par 
tfi (y?) drt , de sorte que l'intégrale considérée soit égale à 

jf%(»)rf,, 

et que, par conséquent, représente la probabilité 

relative d une valeur quelconque de y. x pouvant être re- 
gardé comme une fonction des variables y, x\ x'\ etc., 
que nous désignerons par f( y, x', x", . . .), l'intégrale (i) 
prendra la forme 

où y doit varier depuis y = o jusqu'à y =■ yj, et les au- 
tres variables reçoivent toutes les valeurs pour lesquelles 
f (yiX' ,x" ...) est réelle. 
On aura donc 

+(/)= f p[/(r,x^y^■01 rf /(/, ^ ,J '' -^ ) T(xOT(^)■^•^'rfx^■■ ; 

l'intégration , dans laquelle y doit être regardé comme une 
constante, s étendant à toutes les valeurs des variables 
x', jr", etc., pour lesquelles /(y, x\ x",...) prend une 
valeur réelle. 

13. 

L'intégration précédente exigerait, il est vrai, la con- 
naissance de la fonction <p, qui est inconnue dans la plupart 
des cas» Lors même que cette fonction serait connue, le 
calcul surpasserait, le plus souvent, les forces de l'ana- 
lyse. Dès lors il sera impossible d'obtenir la probabilité de 
chacune des valeurs de y, mais il en sera autrement si Ton 
désire seulement la valeur moyenne de y, qui sera donnée 
g. . 2 
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par l'intégrale 

étendue à toules les valeurs possibles de y. 

Si, par la nature de la fonction, ou à cause des limites 
imposées à jr, X 7 , a:", etc., y n est pas susceptible de recevoir 
toutes les valeurs, on devra supposer que $(y) s'annule 
pour toutes les valeurs que y ne peut atteindre, et Ton 
pourra alors étendre l'intégration de — oo à -f- oo . 

Mais l'intégrale 



prise entre des limites déterminées yj et y/, est égale à 

prise depuis j- = yj jusqu'à y = r/ et étendue à toutes les 
valeurs des variables x', x", etc. , pour lesquelles/est réelle. 
Cette intégrale est é^ale, par conséquent, à l'intégrale 



y y (x) <p(j^) <?(•*•")• • • dxdx' dx" . . . , 



dans laquelle y sera exprimé en fonction de x, x', x", etc., 
la sommation s'étendant à toutes les valeurs des variables 
qui laissent y compris entre r, et r/. D'après cela, F inté- 
grale 

/ fit {?)<*/ 

J — OC 

peut se mettre sous la forme 



yy{x) <p(.r') <p(x"). . . dxdx' dx" . . . , 



l'intégration s'étendant à toutes les valeurs réelles de 
x, X*, x", c'est-à-dire depuis x — — co à .r = -h oo , 
x' = — oo à i' = + « , etc. 
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H. 

Si la fonction y se réduit à une somme de termes de la 
forme 

A *"■/..., 
la valeur de l'intégrale 

r+(r)<ri 

étendue à toutes les valeurs de y, c'est-à-dire la valeur 
moyenne de^, sera égale à une somme de termes de la forme 

x*v(x)dx. f x'to[r')dx'. f * x"'<y{x")dx"..., 

c'est-à-dire que la valeur moyenne de y est égale à une 
somme de termes déduits de ceux mêmes qui composent y, 
en y remplaçant ,r a , or'' 3 , x f,y , etc., par lturs valeurs 
moyennes. La démonstration, de ce théorème important 
pourrait facilement se déduire d'autres considérations. 

15. 

Appliquons le théorème précédent au cas où Ton a 

x 7 -f- x fi -f- x" 2 -H . . . 

y~ ^ , 

* y >• -, 

<7 désignant le nombre des termes du numérateur. 

On trouve tout de sui te que la valeur moyenne de y est égale 
à m* (la lettre m ayant toujours la signification de l'art. 5). 
La véritable valeur dey peut être inférieure ou supérieure 
à sa moyenne, de même que la vraie valeur de x* peut, dans 
chaque cas, être inférieure ou supérieure à m*-, mais la pro- 
babilité pour que la valeur fortuite de y ne diffère pas sen- 
siblement de m s , s'approchera sans cesse de la certitude à 
mesure que a deviendra plus grand. Pour le montrer plus 
clairement, comme il est impossible de chercher exactement 
cette probabilité, nous chercherons l'erreur moyenne à 
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craindre en faisant jr = m*. D'après ce qui a été dit (art. 6), 
cette erreur sera la racine carrée de la moyenne de la fonc- 
tion 

[ ; ">') ■ 

Pour la trouver, il suffit d'observer que la valeur moyenne 

x* , \ n K 
d'un terme tel que — est égale à — j (n ayant la même signi- 
fication que dans l'art. H ), et que la valeur moyenne d'un 

. IX 1 *" 1 , , , 2/71* , . 

ternie tel que — — — est égale a — r ; par conséquent, la 
valeur moyenne de cette fonction sera 

De là nous concluons que si le nombre des erreurs irré- 
gulières esfsuffisamment grand, la valeur de m sera repré- 
sentée, avec une grande certitude, par la formule 



/.r- + x' 3 4- . . . 

"' = v — * — 



et V erreur moyenne à craindre dans la détermination du 
carré de m, sera égale à 



m' 



Comme cette dernière formule contient la quantité // , 
si I on veut seulement se faire une idée du degré de préci- 
sion de cette détermination, il suffira d adopter pour la 
fonction <p une hypothèse particulière. 

Si nous prenons , par exemple , la troisième hypothèse des 

art. 9 et 11, cette erreur sera égale à m 2 l/^- Si on le 

préfère, on pourra obtenir une valeur approchée de n k au 
moyen des erreurs elles-mêmes, à l'aide de la formule 

x K -f- x'* -f- x"' -h • • • 

— — . — — — ~ « 
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On peut affirmer généralement qu'une précision deux fois 
plus grande dans cette détermination exigera un nombre 
d'erreurs quadruple, c'est-à-dire que le poids de la déter- 
mination est proportionnel au nombre <r. 

On verrait de la même manière que si les erreurs des 
observations renferment une partie constante, on déduira de 
leur moyenne arithmétique une valeur de la partie constante, 
et cette valeur sera d'autant plus approchée que le nombre 
des erreurs sera plus grand. Dans cette détermination, 

Terreur moyenne à craindre sera représentée par ^ - » 

k désignant la partie constante, cl m Terreur moyenne des 
observations non corrigées de leur erreur constante. Elle 

tn 

sera représentée simplement par si m représente Ter- 
reur moyenne des observations corrigées de la partie con- 
stante (voyez art. 8). 

lt>. 

Dans les art. 12 à 15 nous avons supposé que les 
erreurs .r, x\ x", etc., appartenaient au même genre d'ob- 
servations, de sorte que la probabilité de chacune de ces 
erreurs était représentée par la même fonction. Mais il est 
évident que les principes généraux exposés dans les ar- 
ticles 12 à 14, peuvent facilement s'appliquer au cas plus 
général où les probabilités des erreurs x, x\ x", etc., sont 
représentées par des fonctions diflérentes 

c'est-à-dire lorsque ces erreurs appartiennent à des obser- 
vations qui n'ont pas le même degré de précision. Suppo- 
sons quer désigne Terreur d'une observation dont Terreur 
moyenne à craindre soit m\ x\ a;", etc., celles d'autres 
observations dont les erreurs moyennes à craindre soient 
respectivement //?', m", etc. : alors la valeur moyenne de 
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la somme 
sera 

Or, si l'on sait, par ailleurs, que les quantités m, m', m", etc., 
sont respectivement proportionnelles aux nombres, 
f, etc., la valeur moyenne de l'expression 

OC^ "\ JC * ~v~ JC " |— • • • 



sera égale à m*. Mais si nous adoptons pour m* la valeur que 
prendra cette expression, en y substituant les erreurs 
x y x \ x", etc., telles que le hasard les offrira, l'erreur 
moyenne qui affecte cette détermination sera , d'après l'ar- 
ticle précédent, 

yjn* -f- „'« -|- ... _ „,< — — w "< _ . . , 

où n ti", etc., ont la même signification, par rapport à la 
seconde et à la troisième observation , que n par rapport à la 
première ; et si l'on peut supposer les nombres n , n\ n", etc. ? 
proportionnels à m, m', m /; , etc. .celte erreur moyenne à 
craindre sera égale à 

1 + • • . 

mais cette manière de déterminer une valeur approchée 
de m n'est pas la plus avantageuse. 

Considérons l'expression plus générale 

dont la valeur moyenne sera aussi m% quels que soient le;, 
coefficients a', a", etc. L'erreur moyenne à craindre lors- 
qu'on substitue la valeur m* à une valeur de j-, déterminée 
d'après les erreurs fortuites x, x", etc., sera, d'après 
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les principes précédents, donnée par la formule 

V^~' w «) H- a' 1 {n'* — /;/«) + a"' (/?"' — w"«7+~~ 
I H- aV* -f- a" a" 2 -4- . . 

Pour que celte erreur soit la plus petite possible, il faudra 
poser 

„» _ w * 

a = " t - t v, f* » 

«* — /«* 

a = ri"' 



Ces valeurs ne pourront pas s'évaluer tant qu'on ne eon- 

naîtra pas les rapports — , — ■ ? etc. Dans l'ignorance où l'on 

est de leur valeur exacte (*)-, le plus sûr sera de les sup- 
poser égaux entre eux [ voyez art. 11), et l'on aura alors 

c'est-à-dire que les coefficients a',»' 7 , etc., doivent être sup- 
posés égaux aux poids relatifs des diverses observations, en 
prenant pour unité le poids de celle à laquelle correspond 
l'erreur x. Ceci posé, désignons, comme ci-dessus, par a le 
nombre des erreurs proposées; la valeur moyenne de l'ex- 
pression 

x 7 H- a'jc' 7 -h a"*"' -4- . . . 
rs 

sera égale à m% et lorsque nous prendrons, pour la vraie 
valeur de m 5 , la valeur de cette expression déterminée 
au moyen des erreurs fortuites x, jr', x'\ etc., l'erreur 
moyenne à craindre sera 

y'/*' -h «'V -h a' V/"' . . . — g/;/ 1 

( * ) On ne conçoit la possibilité de déterminer exactement « , //, etc., 
que dans le seul eus où , par la nature de lu fonction ç>, les erreurs x , j\ 
jt", etc.,- proportionnelles à «, /a", etc., seraient également probables, 
c est-à-dire le cas où 

( A r ci«c </c- M. Gal'ss.) 
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et, enfin, s'il est permis de supposer les quantités n 9 n\ 
n", etc., proportionnelles à m, m', m'\ etc., cette ex- 
pression se réduira à 




résultat identique à celui que nous avons trouvé dans le cas 
où les observations sont toutes de même espèce. 

17. 

Lorsqu'une observation dont la précision n'est pas abso- 
lue, fait connaître une certaine quantité liée analytique- 
ment à une grandeur inconnue, le résultat de cette obser- 
vation peut fournir pour l'inconnue une valeur erronée, 
mais dans la détermination de laquelle il n'y a rien d'ar- 
bitraire qui puisse donner lieu à un choix plus ou moins 
vraisemblable. 

Mais si plusieurs fonctions de la même inconnue sont 
données par des observations imparfaites, chaque observa- 
lion fournira une valeur de l'inconnue, et l'on pourra éga- 
lement obtenir des valeurs, par la combinaison de plusieurs 
observations. Il y a évidemment une infinité de manières 
d'y parvenir; le résultat sera soumis, dans tous les cas, à 
une erreur possible. Selon la combinaison adoptée, Terreur 
moyenne à craindre pourra être plus ou moins grande. 

La même chose aura lieu si plusieurs quantités obser- 
vées dépendent à la fois de plusieurs inconnues. Selon que 
le nombre des observations sera égal au nombre des incon- 
nues, ou plus petit ou plus grand que ce nombre, le pro- 
blème sera déterminé, ou indéterminé, ou plus que déter- 
miné (du moins en général), et, dans ce troisième cas, les 
observations pourront être combinées d'une infinité de 
manières pour fournir les valeurs des inconnues. Parmi ces 
combinaisons, il faudra choisir les plus avantageuses, c'est- 
à-dire celles qui fournissent des valeurs pour lesquelles 
l'erreur moyenne à craindre est la moindre possible. Ce 
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problème est certainement le plus important que présente 
l'application des mathématiques à la philosophie naturelle. 

Nous avons montré, dans la. Théorie du Mouvement des 
Corps célestes, comment on trouve les valeurs les plus pro- 
bables des inconnues lorsque Ton connaît la loi de proba- 
bilité des erreurs des observations, et comme, dans presque 
tous les cas, cette loi par sa nature reste hypothétique, 
nous avons appliqué cette théorie à l'hypothèse très-plau- 
sible, que la probabilité de Terreur x soit proportionnelle 
à e~ h%x% \ de là cette méthode que j'ai suivie, surtout dans 
les calculs astronomiques, et que maintenant la plupart 
des calculateurs emploient sous le nom de Méthode des 
moindres carrés. 

Dans la suite, Laplace, considérant la question sous un 
autre point de vue, montra que ce principe est préférable 
à tous les autres, quelle que soit la loi de probabilité des 
erreurs, pourvu que le nombre des observations soit très- 
grand. Mais lorsque ce nombre est restreint, la question 
est encore intacte; de sorte que, si Ton rejette notre loi 
hypothétique , la méthode des moindres carrés serait préfé- 
rable aux autres, par la seule raison qu'elle conduit à des 
calculs plus simples. 

Nous espérons donc être agréable aux géomètres en dé- 
montrant, dans ce Mémoire , que la méthode des moindres 
carrés fournit les combinaisons les plus avantageuses des 
observations, non-seulement approximativement, mais en- 
core d'une manière absolue, et cela quelle que soit la loi de 
probabilité des erreurs et quel que soit le nombre des ob- 
servations, pourvu que I on adopte pour l'erreur moyenne, 
non pas la définition de Laplace, mais celle que nous avons 
donnée dans les art. 5 et B. 

Il est nécessaire d'avertir ici que, dans les recherches 
suivantes, il ne sera question que des erreurs fortuites dimi- 
nuées de leur partie constante. C'est à l'observateur qu'il 
appartient d'éloigner soigneusement les causes d'erreurs 
constantes. 
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Nous réservons pour une autre occasion l'examen du cas 
où les observations sont affectées d'une erreur constante 
inconnue, et nous traiterons cette question dans un autre 
Mémoire. 

18. 

PROBLÈME. 

Soit U une Jonction donnée des inconnues Y, V, V", etc. ; 
on demande l'erreur moyenne M à craindre dans la dé- 
termination de la valeur de U , lorsque, au lieu des véri- 
tables valeurs de V, V, Y", etc., on prend les valeurs 
déduites d'obseivations indépendantes les unes des autres y 
m, m', m", etc., étant les erreurs moyennes gui corres- 
pondent à ces diverses obseivations . 

Solution. — Désignons par e, e', e", etc., les erreurs des 
valeurs observées Y, Y', Y", etc. 5 l'erreur qui en résultera, 
pour la valeur de la fonction U , pourra s'exprimer par la 
fonction linéaire 

+ a'c'+ a'V 4- . . • ~- E , 

. -, V / ' 1 j - • ' f/Ll 

ou A, A , À , etc. , représentent les dérivées —, — etc. , 

lorsqu'on y remplace V, Y', Y", etc., par leurs vraies valeurs. 

Cette valeur de E est évidente si l'on suppose les obser- 
vations assez exactes pour que les carrés et les produits des 
• ri eurs soient négligeables. 11 résulte de là que la valeur 
moyenne de E est nulle, puisque Ton suppose que les er- 
reurs des observations n'ont plus de partie constante. Or 
Terreur moyenne M, à craindre dans la valeur de U, sera 
la racine carrée de la valeur moyenne de E 2 , c'est-à-dire 
que M* sera la valeur moyenne de la somme 

mais la valeur moyenne de A s e* est A* m 8 , celle de A'V 
est A"m /S , etc., enfin les valeurs moyennes des produits 
?.AA'eV sont toutes nulks; donc on aura 
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11 est bon d'ajouter plusieurs remarquas à celle solution. 

I. Puisqu'on néglige les puissances des erreurs qui sont 
supérieures à la première, nous pourrons , dans notre for- 
mule, prendre pour A, X', etc., les valeurs des coefficients 

différentiels etc., déduites des valeurs observées V, V, 

V", etc. Toutes les fois que U est une fonction linéaire, 
cette substitution est rigoureusement exacte. 

IL Si, au lieu des erreurs moyennes, on préfère intro- 
duire les poids des observations, supposons que p^p', 
p'\ etc., soient les poids respectifs, l'unité étant arbi- 
traire , et P le poids de la valeur de U; on aura . 




III. Soit T une autre fonction de V, V, V", etc. ; posons 

cVT_ dj _ , dj __ „ 
rfV"*' r/V'" - *' HT' ~ 7 ' 

L'erreur commise surT, en adoptant pour V, V, Y 7 ', etc., 
les résultats fournis par l'observation , sera 

x«?h-*V -+-xV\..=E', 

et l'erreur moyenne à craindre dans cette détermination sera 

11 est évident que les erreurs E et E' ne seront pas indépen- 
dantes l'une de l'autre, et que la valeur moyenne du pro- 
duit EE' ne sera pas nulle comme la valeur moyenne de et*': 
elle sera égale à 

xlm 5 ~h -/l'm" + x" À" ai" 2 -h 

IV. Le problème comprend le cas où les valeurs des quan- 
tités V, V, \ 7/ , etc., ne sont pas données immédiatement 
par l'observation , mais sont déduites de combinaisons quel- 
conques d'observations directes. Pour que cette extension 
soit légitime , il faut que les déterminations de ces quantités 
soient indépendantes , c'est-à-dire (ui'ellcs soient fournies 
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par des observations différentes. Si cette condition d'indé- 
pendance n'était pas remplie, la formule qui donne la valeur 
de M ne serait plus exacte. Si , par exemple, une même ob- 
servation était employée, à la fois, dans la détermination 
de V et dans celle de V, les erreurs e et e ne seraient plus 
indépendantes, et la valeur moyenne du produit ee' ne se- 
rait plus nulle. Si l'on connaît, dans ce cas , la relation qui 
lie V et V aux résultats des observations simples dont ils 
dérivent, on pourra calculer la valeur moyenne du pro- 
duit ee', comme il est indiqué dans la remarque III, et dès 
lors corriger la formule qui donnera M . 

19. 

Soient V, V, V', etc., des fonctions des inconnues x, 
y, etc.; soient u le nombre de ces fonctions, p le nom- 
bre des inconnues; supposons que des observations aient 
donné, immédiatement ou médiatement, L, L/, L", etc., 
pour valeurs des fonctions V, V, V', etc., de manière ce- 
pendant que ces déterminations soient absolument indé- 
pendantes les unes des autres. Si p est plus grand que u, la 
recherche des inconnues est un problème indéterminé. Si 
p = ci, chacune des inconnues jr, y, z, etc., peut être re- 
gardée comme calculée en fonction de V, V, V, etc. ; de 
sorte que les valeurs des premières peuvent être déduites 
des valeurs observées de ces dernières, et l'article précédent 
nous permettra de calculer la précision relative de ces di- 
verses déterminations. Si p est plus petit que gt, chaque 
inconnue x, y, z, etc., pourra être exprimée d une infinité 
de manières, en fonction de V, V, V, etc., et, en général , 
ces valeurs seront différentes; elles devraient coïncider si 
. les observations étaient, contrairement à nos hypothèses, 
d'une exactitude rigoureuse. Il est clair, d'ailleurs, que les 
diverses combinaisons fournissent des résultats dont la pré- 
cision sera, en général , différente. 

D'ailleurs si, dans le deuxième et le troisième cas, les 
quantités \ , V, V", etc., sont telles que ct — p -+- i d'entre 
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elles, ou davantage, puissent être regardées comme des fonc- 
tions des autres, le problème est plus que déterminé relative- 
ment à ces dernières fonctions et indéterminé relativement 
aux inconnues x, jr, r, etc. ; et Ton ne pourrait même pas 
déterminer ces dernières inconnues , quand bien même les . 
fonctions V, V, V", etc. , seraient exactement connues : mais 
nous excluons ce cas de nos recherches. 

Si V, V, V", etc., ne sont pas des fonctions linéaires des 
inconnues, on pourra toujours leur attribuer cette forme, 
en remplaçant les inconnues primitives par leur différence 
avec leurs valeurs approchées, que Ton suppose connues; 
les-erreurs moyennes à craindre dan3 les déterminations 
V = L, V' = L', V'=L",... 

étant désignées respectivement par m, m', m", etc., et les 
poids de ces déterminations, par p,p', />", etc , de telle sorte 
que 

pm* = p'm'* = p"m' h — 

Nous supposerons connus les rapports des erreurs 
moyennes ainsi que les poids, dont l'un sera pris arbitrai- 
rement. Si nous posons enfln 

(▼-L)VÊ = «S (V'-L')v£"' = 

les choses se passeront ensuite comme si des observations 
immédiates, également précises et dont l'erreur moyenne 

aurait pour valeur msjp, avaient donné 

(' = O, I»' = O , V " — o , . . . . 

20. 

PROBLÈME. 

Désignons par v, v', v", etc., les fonctions linéaires 
suivantes des indéterminées x, y, z, etc., 

v = ax -h by -h cz -+-...-+-/, 
</ = a x -f b'y -h c'z -4- . . . -h 
v" = a"x+ b"y + c"z-h . . . 4- t\ 
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Parmi tous les systèmes des coefficients •/.', v. ! \ elc. , 
qui donnent identiquement 

k étant indépendant de x, y, z, etc., trouver celui pour 
lequel x 5 H- x' ? -+- -f- . . est minimum. 
Solution. — Posons 

m- -h « V « V -t- . . = £, 
; 4- />V ■+- b"v" 4- ... — « , 
"'4- r'.-'-f- r'V 4-.. = Ç , 



ç, y;, f seront de» fonctions linéaires de .r, y, z, et I on 
aura 



(3) 



r t = xj^ab b 7 H- z ^ bc 4- . . . H- ^ 



on 



^ ^ =: a 2 4- 4- fl"* 4- • 

et de même pour les autres 2. 

Le nombre des quantités etc., est égal au nom- 

bre rj des inconnues jp, s , etc.; on pourra donc obtenir, 
par élimination, une équation de la forme suivante (*), 

x = A 4- (cca) l 4- (af) rj 4- (a 7 ) Ç -f- . . . , 

qui sera satisfaite identiquement lorsqu'on remplacera 
par leurs valeurs (3). Par conséquent, si l'on 



(*) On verra plus loin la raison qui nous a conduit à désigner les coeffi- 
cients de cette formule par la notation (««) , etc. 

(IVo<e rfe M. Gjuss.) 
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pose 

l a (a a) 4- b (ap) 4- c (ay) -f- . . . = a, 

(4) ) a' ( aa ) -h ô'(af) (a 7 ) 4-. . .= a', 
( (aa) 4- (aft) 4- r" (a 7 ) + . . . = a", 

on aura identiquement 

(5) xi' H- aV a"/' -h ... — x — A. 

Cctte équation montre que parmi Jes diilërcnls systèmes 
de coefficients x, y.', x", etc., on doit compter le système 

x — a , x — a , x — a , . . . . 

On aura d'ailleurs, pour un système quelconque, 

(x — a) c + a') c/-+- (x" — a") f"4-. . .= A — /-, 

et cette équation, étant identique, entraine les suivantes : 

(x — «)« H- (x' — a')*'-h(x" — a")*" 4-. . .= o, 
(x — a)* A' +-(•/'- a 7 ') A" + ... = o, 

(x-a) c 4-(x' — a ')c'4-(x"- a w ) c" . . = o. 

Ajoutons ces équations après les avoir multipliées, respec- 
tivement, par («a), («(3), (07)1 etc., nous aurons, en vertu 
du système (4) , 

(x — «) f -4- (x'— «') «' 4- (x" - a") a" 4- . .= o, 
c'est-à-dire 

x'+x' s + ...= a'+a' î + ...+ (>'.-a) 5 +(x'- a')' + . . . ; 
par conséquent, la somme 

• x l 4- x" 4- *" J 4- • • • 
aura une valeur minimum, lorsque l'on aura 

i t if " 

x = a , x = a , x = a , . . . . 

D'ailleurs cette valeur minimum s'obtiendra de la manière 
suivante. 
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L'équation (5) montre que l'on a 

a « -h a'ot! H- a" -4- . . - = î , 
boc -h £V H- fc"a" -h. . .= o, 
fz+f ï -h c a +. ..= o, 



Multiplions ces équations, respectivement, par (aa), 
(«(3), («7), etc., et ajoutons; en ayant égard aux rela- 
tions (4), on trouvera 

y} -h a' 3 -4- a"' 4-. . .= (aa). 

21. 

Lorsque les observations auront donné des équations 
approximatives 

V — O , v = o , v" = O , . . . , 

il faudra, pour déterminer l'inconnue x, choisir une com- 
binaison de la forme suivante, 

y. v -h y.' / -f- y" v" -f- . . . =rz O , 

telle que l'inconnue x acquière un coefficient égal à 1 , et 
que les autres inconnues se trouvent éliminées. 

Le poids de cetie détermination sera , d'après l'art. 18, 

1 

D'après l'article précédent y on obtiendra la détermina- 
tion la plus convenable , en prenant 

y. — a , •/.' — a', x" =r a", ... ; 

alors x aura la valeur A. On obtiendrait évidemment la 
même valeur sans connaître les multiplicateurs a, a', 
etc., en effectuant l'élimination sur les équations 

; = o, » = o, ç = 0, . 
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le poids de cette détermina lion sera 

i 

et Terreur moyenne à craindre 

m y£ (a«) = m' s[j? (aa) = m" ^(aa) =. . . . 

Une marche analogue conduirait aux valeurs les plus 
convenables des autres inconnues y, etc., qui seront 
celles que l'on obtiendrait en effectuant 1 élimination sur 
les équations 

Ç = o, q = o, S ~ o, 

Si nous désignons par SI la somme 

ou, ce qui revient au même, 

p( V - L)' + P' (V - L' )' + y/' ( V" - L*)' + . . . , 
on aura évidemment 

du da do. 

2 Ç = — , 2 ïj = — , 2 ç = — - î . • • i 

rfj rfz ' 

par conséquent, les valeurs des inconnues, déduites de la 
combinaison la plus convenable, et que nous pouvons 
appeler les valeurs les plus plausibles, sont précisément 
celles qui donnent à SI une valeur minimum. Or V — L 
représente la différence entre la valeur observée et la valeur 
calculée j donc les valeurs les plus plausibles des inconnues 
sont celles qui rendent minimum la somme des carrés des 
différences entre les valeurs calculées et observées des 
quantités V, V, V", etc., ces carrés étant respectivement 
multipliés par le poids des observations. J'avais établi 
depuis longtemps ce principe par d'autres considérations 
(Theoria Motus Corporum cœlestium). 

Si Ton veut assigner la précision relative de chacune des 
déterminations, il faut déduire des équations (3), les valeurs 
o. 3 
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de or, /, s, etc., qui se présenteront sous la forme suivante : 
/ x = A -+- (aa) Ç -f- (af) rj H- (ay) Ç -f- . . ., 

(7) J J = BH-(pa)Ç-f-(Pp), 4- (py)Ç -»-..., 

( 8 = C-h(7a)Ç + (7p)ii4-(ry)Ç + --.- 

Les valeurs les plus plausibles des inconnues x,y, z, etc., 

seront A, B, C, etc. Les poids de ces déterminations seront 

i i i 

(17)' (Jf)' M'"" 

et les erreurs moyennes à craindre 

pour x, m')Jp (aa) = m' V^'(aa), . . . , 

pour 7, m sfr (pp) = m' ^(PP), 

pour z, w (77) = m ' {F{rt)i • • • » 

ce qui s'accorde avec les résultats obtenus antérieurement 
( Theoria Motus Corporum cœlestium). 

22. 

Le cas où il n'y a qu'une seule inconnue est le plus fré- 
quent et le plus simple de tous. On a alors 

il sera utile d'en dire quelques mots. 
On aura 

et, par conséquent, 

5 = (/>-*-#>' + f»"+...)*-(^ + /»'L' + ^i: + ...)î 

d'où 

(««) = — : — r~ — » . 

p+p-hp'... 

A== pl+p'V + p"V-h... - 

~ p + p' + p* -+-••• 
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Ainsi , si , par plusieurs observations qui n'ont pas la 
même précision et dont les poids respectifs sont p,p f , 
p", etc., on a trouve*, pour une même quantité, une pre- 
mière valeur L, une deuxième L', une troisième L", etc., 
la valeur la plus plausible sera 

^L+/) , L / + /»T + ... 

et le poids de celte détermination sera 

1> 4-//-h/î"-h 

Si louies les observations sont également plausibles , la 
valeur la plus probable sera 

L-4-L' + L"-h... 
, 

vs 

c'est-à-dire la moyenne arithmétique entre les valeurs ob- 
servées; en prenant pour unité le poids d'une observation 
isolée, le poids de la moyenne sera xs. 



SECONDE PARTIE, 

raKSF.yrtK le a février i8a3, a la société royale î»e cormcui \ 



23. 

Il reste encore à exposer quelques recherches destinées à 
étendre et à éclairer la théorie précédente. 

Cherchons d'abord si l'élimination qui fournit les varia- 
bles x, y y z, etc. , en fonction de £, >j, etc. , est toujours 
possible. Puisque le nombre des équations est égal à celui 
des inconnues, on sait que cette élimination sera possible 
si Sï^? clc * sont indépendants les uns des autres; dans 
le cas contraire, elle serait impossible. 

3. 
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Supposons, pour un instant , que £, » , f , clc* , ne soient 
pas indépendantes , niais qu il existe entre ces quantités l'é- 
quation identique 

o = FH -h G» -h HÇ -\ -f-K; 

nous en conclurons 

F^fl , + G^fl6 + H^flc+.,. = o > 



Posons 

la F -h & G + cH+...= 0, 
J a' F -H V G -f- ✓ H -*-...= G', 

( I} j a"F-+-$"G-+-c"H-h... = G", 



il viendra 

fle + û'0'H-fl"9"-f.., = o, 
^ + ^ + W + ...= o, 
cG -h SB' H- c"G"-K . .= o, 



tQ+l'B' -f. l" Q" + . . ,z= — K. 

En multipliant les équations (i), respectivement par 0, & y 
0", etc., et ajoutant , il vient 

0i + ô' î -+-G" î -f-...= o, 

«t cette équation entraîne les suivantes : 

6 = o, G' = o, 6" = o, 

De là nous concluons, en premier lieu, K = o. En second 
lieu, les équations (i) montrent que les fonctions t>, v', 
t>", etc., sont telles, que leurs valeurs ne changent pas lors- 
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que les variables :r, 7, 2, etc. , prennent des accroissements 
proportionnels à F, G, H, etc. Il en sera évidemment de 
même des fonctions V, V, V, etc. : or cela ne peut avoir 
lieu que dans le cas où il serait impossible de détermi- 
ner jr, y, z, etc. , à Taide des valeurs de V, V, V", etc. . 
lors même que celles-ci seraient exactement connues \ mais 
alors le problème serait indéterminé par sa nature , et nous 
exclurons ce cas de nos recherches * 

24. 

Désignons par j3, fi\ |3", etc., des multiplicateurs qui 
jouent le même rôle relativement à l'inconnue^, que les 
multiplicateurs a , a', a", etc., relativement à l'inconnue x, 
cVst-à-dirc tels, que l'on ait 

a ($*) + b (ft) + c (Pv)H-...= p, 
^(p«)-{-^(PP)4-c'(p 7 )-+-...= P' 1 
« ir (Ppt) + ^pP)+c ,r (P7)+...= P\ 

on aura identiquement 

p^-t-pv-hpv.. .=r — b. 

Soient y, y', /', etc. , les multiplicateurs analogues relatifs 
à la variable z tels, que l'on ait 

a ( 7 a) + 6 (7P)+* (77)-+--™7t 
a' (7a) 4- b' ( 7 p) -f- c' (77) . . .= 7'. 

^( V a)+^(7P)-4-r"( 7 7)H-... = 7' , r 
..........* •••* 

et, par suite, 

7PH _ 7 V-f- /',/'-+-• . .= z — C. 
De la même manière que Ton a trouvé (art. 20) 
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nous trouverons ici 

2p«=o, 2>*='> 2p«-°> 2>'=- B > 

et ainsi de suite. 

On aura aussi , comme dans l'art. 20, 

2>' = (**), = 2>' = (77),---- 

Multiplions les valeurs a, a', a", etc. (art. 20), respecti- 
vement, par |3, JS', etc., et ajoutons \ nous aurons 

c'est-à-dire 

2«P = («P)- 

En multipliant 0, etc., respectivement, par a, a', 

a", etc., et ajoutant, on trouvera 

donc 

(«P) = (P«). 
On trouverait, de la même manière, 

(« y )=( 7 «) = 2«7. (P7) = (7P) = 2PV».«.- 

25. 

Désignons par X, X', X", etc., les valeurs que prennent les 
fonctions f, ^, ^, etc., quand on y remplace x,y, z, etc., 
par leurs valeurs les plus plausibles, A, 13, C , etc., c'est-à- 
dire posons 

aA + iB + cC + ...+ /=)i, 
a'A+i'B+f'C+. .-|-/'=à', 
* 'A + i"B + c"C-f... + /" = r, 



Si nous faisons 

♦ 

y -h V» + V> + . ,. = M, 
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de telle sorte que M soit la valeur de la fonction £1, qui 
répond aux valeurs les plus plausibles des variables -, M 
sera (art. 20) la valeur minimum de £1. 
Par suite, 

a\ + a'V + a" V-h . . . 

sera la valeur que prend £ , lorsque 

x — A t y = B , z = C , • • . . 

Cette valeur est nulle, d'après la manière même dont A , B, 
C, etc., ont été obtenus. On aura donc 

on obtiendrait de même 

^^ôX = o, = o, . . . , 

et 

V a X=o, 2? X = °' 27* = °> 

Enûn, en multipliant les valeurs de X, X', X", etc., respec- 
tivement, par X, X', X", et ajoutant, il viendra 

/; + n'+rri — = v -h x" +v" + ..., 

c'est-à-dire 

- 

26. 

Remplaçons , dans l'équation 

v = ax by -+- cz -+■ . . . ■+• /, 

x> j, z, par les expressions (7) [art. 21], on trouvera, en 
employant des réductions faciles, 

» 

v = a £ -h P » -h 7 Ç-K..4-X, 
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Multipliant respectivement, ou ces équations ou les équa- 
tions (i) de l'art. 20, par X", etc., et ajoutant en- 
suite, on obtient l'identité 

27. 

La fonction A peut se présenter sous plusieurs formes 
qu'il est important d'indiquer. 

Élevons au carré les équations (i) [art. 20], et ajoutons- 
les membre à membre, nous trouverons 

a = x * 2 ** + j2 2 b * + 2 ' 2 c ' + " ' * + 2 ^ 2 ah 
nxz^ac+ lyz ^ bc-\-, . .H- a* ^) fl/ + 2 / 2 

c'est la/; remière forme. 

Multiplions les mêmes équations, respectivement, par ^, 
^, t>", etc. , et ajoutons, on aura 

remplaçons y, ^, etc. , par les valeurs indiquées dans 
l'article précédent , nous trouverons 

ft = |x-4-«r -4- ÇzH-.. . — Aç — Ryj — CÇ — . . . -f- M, 

ou 

& — Ç — A) + n(.r-R)+ç(i-C) + ...+ M; 

c*est la seconde forme. 

Enfin, remplaçons , dans cette seconde forme, x — A, 
y — B, z — C, etc., par les expressions (7) [art. 21 J, nous 
obtenons la troisième forme : 

ft = (aa) ¥-h(W)ri 7 -h(yy)V-h. . .-f-2(ap)Ç„ 
+ 2(«7)ÇÇH-a(P7)uÇ-h...-+-M. 

On peut donner une quatrième forme qui résulte évi- 
demment de la troisième et des formules des articles précé- 
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dents, 

n = + (</ - v)> + (</'- yy + . . (+ m, 

c'est-à-dire 

ft = MH- 2 ("->)'• 

Sous cette dernière forme on voit clairement que M est la 
valeur minimum de ft. 

28. 

Soient e, e / , e", etc., les erreurs commises dans les ob- 
servations qui ont donné 

V = L, V' = L', V"=L",..- 

Les vraies valeurs des fonctions V, V, V", etc., seront alors 

L — €y L' — e* , L" — <?",..., 

et les vraies valeurs de f", etc. , seront respectivement 

par conséquent, la véritable valeur de x sera 

A — a<? v(/t? — aV \fp r —a"e"ffî—..., 

et Terreur commise dans la détermination la plus conve- 
nable de l'inconnue x sera , en la désignant par Ex, 

E* r= <xe -h ftV v9" -H <*V' vf*" + • • • • 

De même, Terreur commise dans la détermination la plus 
convenable de la valeur de j- sera 

Ej = sfc 4- fV y 7 ? -+- P"'" V / ? / + - • • • 
La valeur moyenne du carré (Ex)* sera 

tri 1 p (a 1 -f- a'» -h a" 3 + . . . ) = /«'/>(««)• 
La valeur moyenne de (Er) s sera de même 



Digitized by Google 



( 4* ) 

connue nous l'avons déjà reconnu. On peut également 
assigner la valeur moyenne du produit Ex. Ey, qui sera 

m 7 p ( a p + a' p' -f- a" (J" -H . . . ) = mjp (ap). 

On énonce ces résultats plus brièvement de la manière sui- 
vante : 

Les valeurs moyennes des carrés (Er)', (E^)*,etc, sont 

respectivement égales aux produits de par les quo- 
tients différentiels partiels du second ordre 

d : ïi d'à 

et la valeur moyenne d'un produit tel que Ex . E/ est le 

produit de ^m*p par ^ ? en regardant H comme fonc- 
tion de £, »7, f , etc. 

29. 

Soit t une fonction donnée et linéaire des quantités x, 
y, z, etc., par exemple, 

t =/x -h gj -h /<z -h . . . -h * ; 

la valeur de £ déduite des valeurs les plus plausibles de 
x, y, s, etc., sera 

nous la désignerons par K. En désignant par Et Terreur 
commise en l'adoptant, on aura 

Et =/Ex-h gEj -+- A Es -h. . .; 

la valeur moyenne de cette erreur sera évidemment nulle, 
c'est-à-dire que Terreur ne contiendra pas de partie con- 
stante, mais la valeur moyenne de (Et)*, c'est-à-dire de la 
somme 

/ , [Ej) 5 + *fgEx.Ey -h 2///Ex.Ez-i-.. . 

+ A , (E«) , +... I 
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sera, d'après l'article précédent , égale au produit de m f p 
par la somme 

/»(aa) •+■ 2/^(ap ) + 2/// («y) -h . . . 

c'est-à-dire au produit de m*/; par la valeur de la fonction 
il — M, lorsqu'on y fait 

l =f, i = s, ç = /f, 

Désignons parw cette valeur de II — M ; l'erreur moyenne 
à craindre, lorsque l'on prend £ = K, sera m^ptù, et le 

poids de cette détermination sera 
Puisque Ton a identiquement 

a — M = (x — A) Ç -h (y — B) îj -4- (a — C)Ç+..., 
w sera égal à la valeur de l'expression 

(x - A)/-*- ( j- - B ) £ ( t - C) h -4- . . . 

[qui représente (f — K)]î dans laquelle on remplacera 
x 9 y, z, etc., par les valeurs correspondantes à £ = 
Y) = g, £ = ^, etc. 

Enfin , observant que f , exprimé en fonction des quan- 
tités £, >}, f , etc., aura K pour partie constante, si l'on sup- 
pose 

t = FÇ -f- G* -h HÇ+... + R, 

on aura 

w = /F -h g G -H /ï H -f- 

30. 

Nous avons vu que la fonction il acquiert son minimum 
absolu M, lorsque l'on y fait 

■r = A, J =B, s — C, . . . , 
ou , ce qui revient au même, 

% j=z O , » = o , Ç -= O , . . . . 
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Si lou attribue à l'une des inconnues une autre valeur ; que 
l'on fasse, par exemple, 

x = A -h A, 

les autres inconnues restant variables, ft pourra acquérir 
une valeur minimum relative, qui s'obtiendra à l'aide des 
équations 

da dCl 
X==AH ~ A ' 7fy = ° y dz = °"" 

et, par suite, 

» = o , Ç = o, . . . ; 

or, puisque 

* = À + (aa)Ç+(ap)nH-(«7)Ç-*-.. M 

on en conclut 



(«0 

On trouvera de même 

(«7) 

La valeur minimum relative de ft sera 

A» 



HP + H = M-f- 



(««) 



Nous en conclurons, réciproquement, que si il ne doit pas 
surpasser M -f- /i a , la valeur de x est nécessairement com- 
prise entre les limites A — ^^(aa), et A -h fz \/(aa). Il 

est important de remarquer que p <J(*ol) devient égal à 
l'erreur moyenne à craindre dans la valeur la plus plau- 
sible de x, si Ton pose 

c'est-à-dire si ^ est Terreur moyenne d'observations telles, 
que leur poids soit l'unité. 

Plus généralement, cherchons la plus petite valeur de la 
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fonction il qui puisse correspondre à une valeur donnée 
de t , t désignant , comme dans l'article précédent , l'expres- 
sion linéaire 

% f* + gy hz+ . . . + x, 

dont la valeur la plus plausible est K ; désignons par K-t-y. 
la valeur donnée de t. D'après la théorie des maximum et 
minimum, la solution du problème sera donnée par les 
équations 

d£i _ de 
dx dx ' 
da dt 
djr dy 
d£l _^dt 
dz dz* 



0 désignant un multiplicateur encore indéterminé. 

Si, comme dans l'article précédent, nous posons identi- 
quement , 

/ = F5 + G«-l-HÇ + ... + K, 

nous aurons 

K + x = G(/F+ffG + AIi + ...) + £; 

d'où 

w 

w ayant la même signification que dans l'article précédent. 

Puisque — M est une fonction homogène du second 
degré, par rapport aux variables £, yj, etc., sa valeur 
pour 

sera évidemment 

6' 6), 
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et, par conséquent, la valeur minimum de H, lorsque 

sera 

/? # 



Réciproquement, si doit rester inférieure à une valeur 
donnée M + |x. ! , la valeur de t sera nécessairement com- 

. prise entre les limites R — u v'o), K-h^y/c»), et fW^ sera 
Terreur moyenne à craindre dans la valeur la plus plau- 
sible de t , si {A représente Terreur moyenne d'observations 
dont le poids serait l'unité. 

3). 

Lorsque le nombre des inconnues x, j', etc., est un 
peu grand , la détermination des valeurs numériques de 
A, B, C, etc., par l'élimination ordinaire, est assez pénible. 
C'est pourquoi nous avons indiqué, dans la Théorie du 
Mouvement des Corps célestes, et développé plus tard, 
dans le Mémoire sur les éléments de l'orbite de Pallas 
[Commentaires de Gottingue, t. I), une méthode qui sim- 
plifie ce travail autant que possible. 

La fonction £1 doit être ramenée à la forme suivante : 



u 03 m" 2 u m > 



a>° " ui/ ' z" ' or M ' 

où les diviseurs <A>°, ifi/, CD'", etc., sont des quantités 
déterminées; m 0 , m', il", etc., sont des fonctions linéaires 
de x,y y z, etc., telles que la seconde ii ne contient pas x* 
la troisième u" ne contient ni a: ni j-, la quatrième ne con- 
tient ni x, ni y, ni z, et ainsi de suite, de sorte que la der- 
nière u(° — 0 ne contient que la dernière des inconnues 
j, etc.; enfin les coefficients de x, y, z, etc., dans 
u°, m', zi", etc., sont respectivement égaux à «Ao 0 , atë/, O^etc. 
Alors on pose 

u° — o , «' = o , «" — o, = o , . . . , 

et Ton aura très-facilement les valeurs de x> y, z, etc., 
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en résolvant ces équations et commençant par la dernière. 
Je ne crois pas nécessaire de répéter de nouveau l'algorithme 
qui conduit à cette transformation de la fonction II (*). 

Mais l'élimination qu'il faut effectuer pour trouver les 
poids de ces déterminations exige des calculs bien plus 
longs encore. Nous avons montré, dans la Théorie du Mou- 
vement des Corps célestes, que le poids de la détermination 
de la dernière inconnue, qui entre seule dans m(*~0, 
est égal au dernier terme de la série des diviseurs JU, 
ift/, G", etc. Cette recherche est facile ; aussi plusieurs cal- 
culateurs, voulant éviter , une élimination pénible, ont eu 
l'idée, faute d'autre méthode, de répéter la transformation 
indiquée en considérant successivement chaque inconnue 
comme la dernière. J'espère donc que les géomètres me 
sauront gré d'indiquer, pour calculer les poids des déter- 
minations, une méthode nouvelle qui ne me semble plus 
rien laisser à désirer sur ce point. 

32. 

Posons 

/ uo = x 9 x -f- ni>° y H- £° z -h . . . 4- -C, 

on aura identiquement 

dïl z= \âx -f- * dr -f- Ç dz -+- . . . 

du 0 u' du' u'du" 

-*-«' (* r -*-i/ </a5+ ---) 

-f- a" (rfs 4- • .) 



. ( * ) On trouvera ces calculs dans une Note à la fin du volume. J. B. 
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d'où nous déduirons : 



X 9 l)b 



Les valeurs de w°, u\ u", etc., déduites de ces équations, 
se présenteront sous la forme suivante : 

( «° = ti 
, 3) 1 «' = A'? -h *», 

«* = A"Ç-hB"»,H-Ç, 



De la différentielle complète de l'équation 

n = Ç(x - A) + r ( (/-B) + !;(3-C)+...+ M, 
retranchons l'équation 

irffï = Ç <ir -f- yj dy Çrfz -H. . . , 

il viendra 

Cette expression doit être identique avec celle que Ton 
obtient à l'aide des équations (3), c'est-à-dire 

d\ 4- 4r ( A'rfÇ 4- + ~ ( AVÇ + B"*ï 4- rfÇ) -K . . ; 



on aura donc 



» ->// " 



(4) {' = S7 + B>+-.. + l, 

a" 

1 — 7777 4-. . . -f- C, 



En substituant dans ce3 expressions, les valeurs de i/% 
u', u", etc., tirées des équations (3), on aura effectué l'éli- 
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mination. Pour déterminer les poids, nous aurons 



(5) 



JL ¥1 — £1 

~~ JL« Di>' G" + 

tJÏ? G 7 ' (D* 

g" ~*~ <eT 



(77) = 



La simplicité de ces formules ne laisse rien à désirer. On 
trouverait des formules également simples pour exprimer 
les autres coefficients (a/3), (ay), (|3y), etc. ; mais, comme 
leur usage est moins fréquent, nous nous dispenserons de 
les exposer. 

33. 

• 

L'importance du sujet nous a engagé à tout préparer 
pour le calcul et à former les expressions explicites des 
coefficients A', A", A'", B", B'", etc. 

Ce calcul peut être abordé de deux manières : la pre- 
mière consiste à reporter, dans les équations (2), les valeurs 
de n°, a', m", etc., déduites du système (3), qui doivent 
rendre ces équations identiques \ et la seconde à exprimer, 
au contraire, que le système ( 2) devient identique lorsqu'on 
y substitue les valeurs de £ , *j, f , déduites du système (3). 

La première méthode conduit aux formules suivantes : 

Kl 

G 9 G' 



të° CE)' ., 00" ... .... 

— H ; A' H 7, A" -f- A'" = o, 



Ces formules feront connaître A', A", A'", etc. 
o. 
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On aura ensuite, 
-Ob' 

_, + ër ,B''+B" = o, 



qui donneront B", B'", etc. -, puis 

CÔ" 

_i_ C'" n 



qui feront connaître C", etc.; et ainsi de suite. 
La seconde méthode donne le système suivant : 

X° A'-hrë>° = o, 

d'où Ton tire A' ; 

X° A" -+- i)b° B" -h C° = o , 

B ï + e'=o, 

d'où l'on tire B" et A" -, 

X° A'" -h ift» 0 B'" -h C° C" -h (D° = o , 
KÎ,'B" , + e'C" / + (D / = o, 
G"C" -f-(D" = o, 

d'où Ton déduira C w , B w , A w ; et ainsi de suite. 

Les deux systèmes de formules offrent des avantages à peu 
près égaux, lorsque Ton veut les poids des déterminations 
de toutes les inconnues .r, y, z, etc.; mais lorsqu'on ne 
cherche qu'une seule des quantités (aa), (|3(3), (yy), etc., 
le premier système est bien préférable. 

D'ailleurs, la combinaison des équations (1) et (4) con- 
duit aux mêmes formules, et fournit, en outre, un 
second moyen d'obtenir les valeurs les plus plausibles 
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A, tt, C, etc., qui sont 

pO (J 1>H ttf/ 

a À'J^ A"^ - A w 



C — ^ C/" — 



L'autre calcul est identique avec le calcul ordinaire dans 
lequel on suppose 

«° = o, #i' = o, a" = o,.... 
34. 

Les résultats obtenus dans l'art. 32 ne sont que des cas 
particuliers d'un théorème plus général qui peut s'énoncer 
de la manière suivante : 

Théouème — Si t représente la fonction linéaire sui- 
vante des inconnues x, y, z, etc., 

t=fx 4- &r 4- hz-h. . .+ k, 
dont l expression en fonction des variables u°, u', u", etc., 
soit 

t = A°u 9 4- k' u! 4- A"u" 4- . . . 4- K, 

K sera &i valeur la plus plausible de l, et le poids de cette 
détermination sera 

i 

Démonstration, — La première partie du théorème est 
évidente, puisque la valeur la plus plausible de t doit cor- 
respondre aux valeurs 

»■ = o , «' = o , V = o , . . . . 

Pour démontrer la seconde partie, remarquons que l'on a 

i r/fl = Idx+ndy 4- Çdz-h. . . , 
dt =fdx + gdy + hdz + . . , 

4- 
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cl par conséquent, lorsque 

£=/» *=g> Ç = /<,... , 

on a 

quelles que soient les différentielles dx y dy y dz^eie. Il suit 
de là qu'en supposant toujours , 

l —f> * = g > ç = A , . . . , 

on aura 

= k°du°+k t du'+h"du"-h.... 

Or on voit facilement que si les différentielles dx, dy, 
dz, etc. , sont indépendantes les unes des autres, il en sera 
de même de du 0 , du\ du f \etc. \ nous aurons, par conséquent, 
pour 

l =f, r =g, f =-. A,. . 
U»z=X°k\ tt' = lft>'A', u"= e"A", 

Par conséquent, la valeur de Q. correspondant aux mômes 
hypothèses, sera 

X° A 0f -f- ifï/A'* -h 0"h nt -h ... -H M : 

ce qui , d'après l'art. 29, démontre l'exactitude de notre 
théorème. 

Si d'ailleurs on désire effectuer la transformation de la 
fonction f, sans avoir recours aux formules (4) (art. 32), on 
a immédiatement les relations 

g = lfl, 0 Â 0 +lfï/A', 

qui permettront Tie déterminer A 0 , A', A", etc., et nous au- 
rons enfin 

K = — £>A° — £A' — £'A" — . ... 
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35. 

Nous traiterons particulièrement le problème suivant , 
tant à cause de son utilité pratique, que de la simplicité de 
la solution : 

Trouver les changements que les valeurs les plus plau- 
sibles des inconnues subissent par V adjonction d'une nou- 
velle équation, et assigner les poids de ces nouvelles 
déterminations. 

Conservons les notations précédentes. Les équations pri- 
mitives , réduites à avoir pour poids l'unité, seront 

v = o , f' = o , e" = o , . ; 

on aura 

Çî etc -î seront les dérivées partielles 

♦ 

d£l dSl di\ 

lîdx Tdy 2dz y " 9 

et enfin on aura , par l'élimination , 

x = A-f- (aa) \ 4- («P)u + (a 7 ) . . , 

y = B + (ap) Ç -f- (pp) « + (p 7 ) Ç +. . , 

« = C + («7)Ç-4-(P7)n + (77) Ç-H. • • » 



Supposons maintenant que l'on ait une nouvelle équa- 
tion approximative, 

dont nous supposerons le poids égal à l'unité. Cherchons 
les changements que subiront les valeurs les plus plausibles 
A , B, C , etc., et celles des coefficients (aa), (j3/3 ), etc. 
Posons 

a 4- t/*' = n*, 
1 du* 1 rfn* É 1 dn* 

et soit 

x = A* -h (aa*) + (ap*} n* + (ay*)Ç* 4- • . . 
le résultat de l'élimination. 
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Soit enfin , 

=fx -h gy 4- hz 4- . . . -h A , 
qui deviendra, en ayant égard aux équations (i), 
p*=F$+ G«+HÇ4-..4-K, 

et posons 

F/4- G# 4- HA -f-. . .= »; 

K sera évidemment la valeur la plus plausible de la fonc- 
tion u*, telle qu'elle résulte des équations primitives, sans 
avoir égard à la valeur o fournie par la nouvelle observa- 

tion , et - sera le poids de cette détermination. 

Or nous avons 

et , par suite , 

F g * 4- G n* 4- H Ç * + . . . H- K = ( i -h F/H- G g 4- H A -h . . . ) ; 

d'où Ton déduit 

r *_ ^ + 0"' + Hg > +...-HK 

I 4- w 

On a , en outre, 

x = A-f-(aa)Ç*H-(ap)>ï*+(a 7 )Ç*4- 

-**[/(«*)] + *(«(0 + M«7) + -..] 
= A 4- (aa) **4- . • . - F <>* , 

F 

= À + (aa)É*+(«PK-K..— -(F$*4- 0**4- II Ç* 4-...K). 

1 -t- w 

Nous déduirons de là , 

A* = A-™, 
I 4- w 

qui sera la valeur la plus plausible de x, déduite de /ou/es 
les observations. 



Digitized by Google 



( 55 ) 

On aura aussi 

F* 

(cwt*) — (aa) ; 

par conséquent, 



F 3 

(««) - - 



sera le poids de cette détermination. 

On trouvera de la même manière , pour valeur la plus 
plausible dey, déduite de toutes les observations, 

_ . _ GK 

B* = B \ 

l -h w 

le poids de celte détermination sera 



(pp; « + » 

et ainsi de suite. 

Le problème est donc résolu. 
Ajoutons quelques remarques. 

I. En substituant les nouvelles valeurs A*, B*, C*, etc., 
Ja fonction v* obtiendra la valeur la plus plausible 

K (F/+G^ + HA + ...) = JL., 

et, puisque Ton a, identiquement, 

. .p = -ï- 5 * + »* + — - 5 * 

l + w l + w I -h w I -h w 

le poids de cette détermination sera (art. 29) 

I -h w I 

= -+ !• 



F/ 4- G g +HA + .'.. '« 

Ces résultats pourraient se déduire immédiatement des rè- 
gles exposées à la fin de l'art. 21 . L'ensemble des équa- 
tions primitives avait, en effet, fourni la détermination 

^*=K, dont le poids était ~j une observation nouvelle 
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donne une autre détermination v>* = o, indépendante de la 
première, ayant pour poids l'unité ; leur combinaison don- 
nera la détermination 

K 



qui aura pour poids 



I -4- eu 



I 



IL On conclut de ce qui précède que, pour 
* = A*, r = B*, * = CV , 
on devra avoir 

V = o 7 n* = o, Ç* = o, .., 

et y par suite, 

r g* r 

s — ; — 5 » = j ç = — , . . . . 

1 -h w l-hw i-f- w 

Puisque d'ailleurs 

n = Ç(x - A) 4- n(jr — B) -h — C) 4-. . .H- M, 
£1*= fi -h p*% 

on devra avoir, pour ces mêmes valeurs, 



EL En comparant ces résultats avec ceux de l'ar- 
ticle 30, nous voyons ici que la fonction D. a la plus petite 
valeur qu'elle puisse obtenir lorsqu'on s'impose la condition 



i -f- « 
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36. 

Nous donnerons seulement ici la solution du problème 
suivant, qui a beaucoup d'analogie avec le précédent; mais 
nous nous abstiendrons d'indiquer la démonstration , à la- 
quelle le lecteur suppléera facilement en s'aidant de ce qui 
précède. 

Trouver les changements des valeurs' les plus plausi- 
bles des inconnues et les poids des déterminations nou- 
velles, lorsque Von change le poids de l'une des observa- 
tions primitives. 

Supposons qu'après avoir achevé le calcul on vienne à 
remarquer qu'on a attribué à une observation, à la pre- 
mière par exemple qui a donné V = L , un poids trop fort 
ou trop faible, et qu il serait plus exact de lui attribuer 
le poids p*, au lieu du poids p : il ne sera pas nécessaire 
de recommencer le calcul , mais il sera plus commode de 
former les corrections à l'aide des formules suivantes. 

Les valeurs les plus plausibles des inconnues deviendront 



x = A 




(p*- P )*\ 


P + (P* 


— p) (a* -f- bp H- cy -+■ . . .) 






(p' - p) p 


P + (l>* 


— p) (a a H- b p -f- cy -+• , . , ) 


z = C 




( P '-P)y\ 




P + (P* 


— p)(aot + +cy +. . .) 



les poids de ces déterminations s'obtiendront en divisant 
l'unité, respectivement, par 

(aa) (P*-P)«> •_ 

P+ {P*—P) («« + *P + n +-..)' 



(77)- 



P *+* (P* — P) (* a a P . .) 

(P*-P)'f 



p -h (p* — p) («a + ^ + cy-f-...) 
Cette solution convient au cas où, après avoir achevé le 
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calcul , il faudrait rejeter tout à fait Tune des observations, 
puisque cela revient à faire p* = oj de même, p* = ao 
conviendra au cas où l'équation V = L , qui dans le cal- 
cul avait été regardée comme approchée, serait rigoureu- 
sement exacte. 

Si, après le calcul terminé, plusieurs équations nouvelles 
venaient s'adjoindre aux proposées, ou si les poids attribués 
à plusieurs d'entre elles étaient erronés , le calcul des cor- 
rections deviendrait trop compliqué, et il serait préférable 
de tout recommencer. 

37. 

Nous avons donné, dans les art. 15 et 46, une méthode 
pour déterminer, approximativement, la précision d'un 
système d'observations (*); mais cette méthode suppose 
connues exactement les erreurs réelles que l'on a effective* 
ment rencontrées dans une suite nombreuse d'observa- 
tions; or cette condition n'est remplie que bien rarement, 
pour ne pas dire jamais. 

Si les quantités dont l'observation fournit les valeurs 
approchées dépendent, suivant une loi donnée, d'une ou de 
plusieurs inconnues, on pourra trouver, par la méthode 
des moindres carrés, les valeurs les plus plausibles de ces 
inconnues; si, dès lors, on calcule les valeurs correspon- 
dantes des grandeurs observées, ces dernières pourront être 
regardées comme différant peu des véritables : de sorte que 
leurs différences avec les valeurs observées, représenteront 
les erreurs commises avec une certitude d'autant plus 
grande , que les observations seront plus nombreuses. 
Telle est la marche suivie dans la pratique par les calcu- 

(*) Les recherches sur le môme sujet insérées par nous (Zeitschrift /ur 
Astronomie und verwandle Wissenschqfton , vol. I, page i85) sont fondées 
sur l'hypothèse relative à la probabilité des erreurs à laquelle nous avions 
été conduit dans la Théorie du Mouvement des Corps célestes. 

(Note de M. GaissO 

On trouvera ce Mémoire à la fin du volume. J. R. 
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la.teurs, qui ont essayé, dans des cas compliqués, d'évaluer 
à posteriori la précision des observations. Quoique suffi- 
sante dans bien des cas, cette méthode est , théoriquement, 
inexacte et pourrait quelquefois conduire à de graves er- 
reurs ; c'est pourquoi il est très-important de traiter la 
question avec plus de soin. 

Conservons les notations de l'art. 19. La méthode dont 
il s'agit consiste à regarder A, B,C, etc., comme les véri- 
tables valeurs des inconnues jr, etc., et X, X', X", etc., 
comme celles des fonctions y, v etc. Si toutes les obser- 
vations ont une égale précision et que leur poids commun 

p = p' = p"=... 

soit pris pour unité, ces mêmes quantités, changées de signe, 
représentent, dans celte supposition, les erreurs des obser- 
vations, et, par conséquent, d'après Fart. 15, 

sera Terreur moyenne des observations. Si les observations 
n'ont pas la même précision , — X, — X', — X", etc. , re- 
présenteront les erreurs des observations, respectivement 
multipliées par les racines carrées des poids , et les règles 
de l'art. 16 conduiront à la même formule, 

qui exprime déjà l'erreur moyenne de ces observations 
lorsque leur poids est égal à l'unité. 

Mais le calcul exact exigerait évidemment que l'on rem- 
plaçât X, X', X", etc. , par les valeurs de \>\ v", etc. , dé- 
duites des véritables valeurs des inconnues x,y, z, etc., 
et la quantité M par la valeur correspondante de fl. 
Quoique l'on ne puisse pas assigner cette dernière valeur, 
nous sommes certain pourtant qu'elle est plus grande 
que M qui est son minimum : elle n'atteindrait cette li- 
mite que dans le cas, infiniment peu probable, où les va- 
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leurs véritables des inconnues se confondraient avec les plus 
plausibles. Nous pouvons donc affirmer, en général, que Fer- 
reurmoyenne calculée par la pratique ordinaire est plus pe- 
tite que l'erreur moyenne exacte, et que l'on attribue , par 
conséquent , aux observations une trop grande précision. 
Voyons ce que donne une théorie rigoureuse. 

38. 

Avant tout, il faut chercher comment la quantité M dé- 
pend des véritables erreurs des observations. Désignons ces 
erreurs, comme dans Fart. 28, pare, e 7 , e", etc., et po- 
sons , pour plus de simplicité , 

et 

m sjp = m! \Fp' = m" s/p 7 ' = . . . = p. 

Soient 

A — B — j», C— z°, 

les vraies valeurs des inconnues J?,/, z, etc., pour les- 
quelles £ , m , £, etc. , soient, respectivement, — £° , — yj°, 
— £°, etc. Les valeurs correspondantes de v y v*, t^', etc., se- 
ront évidemment 

— i, — s', —«",...; 

de sorte qu'on aura 

H°=fl6H-«V-4-rt"«"-}-..., 
r° = be + b's! -h b" s" -+-. . . , 

Ç°=rC£ -+- cY -h C"ï" -h. . . , 



z*~ as -f- aV -h a" s" -+- . . . , 
/•=pt-f. pV 4- p" f" -+-'..., 



enfin 

sera la valeur de la fonction ft, correspondant aux vraies 
valeurs des variables a:, y, z y etc. 
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Puisque Ton a identiquement' 
a = M + (x — A ) Ç H- (y - B) * -h [z — C) K -+- . . . , 

on aura aussi 

M = n« — — y* — z 4 ? — 

De là , il résulte évidemment que M est une fonction homo- 
gène du deuxième degré des erreurs e, e\ e", etc. -, cette fonc- 
tion pour diverses valeurs des erreurs pourra devenir plus 
ou moins grande. Dans l'ignorance où nous sommes de ces 
valeurs, il est bon d'examiner attentivement la fonction M, 
et de calculer d'abord sa valeur moyenne d'après les prin- 
cipes du calcul des probabilités. Nous obtiendrons cette 
valeur moyenne en remplaçant les carrés e f , etc., par 
m*, m'*, etc., et en omettant les termes en eé ', ee", etc., dont 
la valeur moyenne est zéro*, ou, ce qui revient au même, 
en remplaçant chaque carré e*, e'% e"*, . . . , par jx*, et en né- 
gligeant ce', ee", .... D'après cela, letermell 0 fournira gj^*; 
le terme — .r° £° donnera 

— (a a a! a! + a" a" . . .) p* = — f ; 

chacune des autres parties donnera également — , de sorle 
que la valeur moyenne totale sera 

gj désignant le nombre des observations, et p le nombre des 
inconnues. La vraie valeur de M pourra , suivant les cas 
que le hasard présentera, être plus grande ou plus petite 
que cette valeur moyenne , mais la différence sera d'autant 
moindre que le nombre des observations sera plus grand ; 
de sorte que 

V ° — p 

pourra être regardé comme une valeur approchée de f-t : par 
conséquent la valeur de p , fournie par la méthode erronée 
dont nous avons parlé dans l'article précédent, devra être 

augmentée dans le rapport de \Jn — p à y gt. 
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« 

39. 

Afin de faire voir plus clairement jusqu'à quel point il es! 
permis de regarder la valeur de M, fournie par les obser- 
vations, comme égale à la valeur exacte, il faut chercher 
quelle est Terreur moyenne à craindre lorsque l'on fait 

» M 

u 7 = 

r CT p 

Cette erreur moyenne est la racine carrée de la valeur 
moyenne de la quantité 



que nous écrirons ainsi : 

et comme la valeur moyenne du second terme est évidem- 
ment nulle, la question se réduit à chercher la valeur 
moyenne de la fonction 

*- * 

Désignons cette valeur moyenne par N, l'erreur moyenne 
cherchée sera 



Si l'on développe la fonction ¥, on voit qu'elle est une 
fonction homogène des erreurs e, e', e", etc., ou, ce qui re- 
vient au même, des quantités e, e', e", etc.; on trouvera 
donc la valeur moyenne : 

i°. En remplaçant les quatrièmes puissances e 4 , e'*, 
e fn , etc., par leurs valeurs moyennes; 
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2°. En remplaçant les produits e*<;'% e*c ni , etc., par 
leurs valeurs moyennes, c'est-à-dire par m* m", m* m"% etc.; 

3°. En négligeant les produits tels que e 8 e', c* e'e'\ etc. 

Nous supposerons (art. 46) les valeurs moyennes de e 4 , 
e", e ,n y etc., proportionnelles à /«*, w /4 , m"*, etc, de sorte 

que les rapports des unes aux autres soient — ? v 4 désignant 

la valeur moyenne des quatrièmes puissances des erreurs 
pour les observations dont le poids serait Puni té. 

Les règles précédentes pourraient se traduire de cette 
autre manière : * 

Remplacer chaque quatrième puissance e*, e' 4 , s"*, etc., 
par v 4 -, chaque produit eV*, cV", etc., par ^ 4 , et négliger 
tous les termes, tels que eV ou eVe", eeW". 

Ces principes étant compris, on verra facilement que : 

I. La valeur moyenne de £l oi est 

GTV 4 (d 1 o) |A 4 . 

IL La valeur moyenne du produit e* o:^ 0 est 

auv* + («V-4- tf"a"-t-. . .) p 4 :=«a(v 4 — p 4 ) + p 4 , 

car 

fla + û'a' -f- fl"a" -+-... = i . 
De même, la valeur moyenne de t ft x° £° est 

♦ 

la valeur moyenne de e'^x 0 J° est 

et ainsi de suite. 

Donc la valeur moyenne du produit 

(e' + fi "+6"»4- ...)*• i* ou n°x*Ç f 

sera 

v 4 — p 4 ■+- ©|i 4 . 

Les produits {l°/°y} 0 ou U°a 0 Ç°, etc., auront la même 
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valeur moyenne; donc le produit 

ft^x-p + jV -H »•?+....) 
aura pour valeur moyenne 

pv 1 -4- p (ct — i) ii*. 

III. Afin d'abréger les développements qui vont suivre, 
nous adopterons la notation suivante. Nous attribuerons à la 
caractéristique S un sens plus étendu que nous ne l'avons 
fait jusqu'ici, en lui faisant désigner la somme des termes 
semblables, mais nçn identiques, qui proviennent de toutes 
les permutations des observations. Nous aurons, d'après 
cette notation , 

Calculant par parties la valeur moyenne du terme x°* £ 0 *, 
nous aurons d'abord, pour valeur moyenne du produit 

«'aV 4- a' («" -+- a" 9 -h . . . ) 

De même, la valeur moyenne du produit a" e'* £ os , est 

„'»a" ( v < — _|_ a'^^^ J 

et ainsi de suite. 

Par conséquent, la valeur moyenne du produit 

sera 

Or la valeur moyenne de aa'ee' £°* est 

2 aa' «a' p\ 
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La valeur moyenne de aaV£ 0 * est 

■ 

et ainsi de suite. D'où l'on conclut facilement que la valeur 
moyenne du produit 

est 

Ceci posé, nous aurons pour valeur moyenne du pro- 
duit x oî J°% 

! 

IV. On trouvera d'une manière analogue, pour valeur 
moyenne du produit x°j° £°y; 0 , 

Or on a 

2"«*'P'=2 fla -2*p^2 <ia *P' 

^b^'W^ab^^-^ab^, 

2 rta==i ' 2*p = i ' 2 flpr=o ' 2* a=o; 

cette valeur moyenne sera donc 

(v*-3fx')2^^4-^(«-+-2^-2 a ?)- 

V. On trouverait, par un calcul semblable, que la valeur 
moyenne dear 0 3 0 £°f 0 est 

(V — 3f*«)2«««7 +^( l +2 /?c, 2 a7 ) ; 
o. 5 



Digitized by Google 



( ™ ) 

et ainsi de suite. En additionnant, on obtient la valeur 
moyenne du produit 

cette valeur est 
(v 4 — 3p 4 )2 + *P + <"7 + •)] 

M. On trouverait de la même' manière 

( v « — 3p«) V [6|J( fl « +*f + f 7 +.. .)] -h (p H- 2)^ 

pour valeur moyenne du produit 

jV(x^ 9 4-/ < V -h 2°ç°-h.. .), 

et 

(v 4 — 3 u*)2 [ r 7 [ a<x + + <"7 + . . .)] (p 2) p 4 
pour valeur moyenne du produit 

sPÇi^çi 4- + **ç , -t-. ..); 
et ainsi de suite. 

Nous aurons donc , par l'addition , la valeur moyenne du 
carré 

elle sera 

( v « _ 3fx 4 }]£ [(*a b$ +q 4-. . .)'] (p* 4- 2p)f*'. 
VII. Nous concluons enfin de tous ces préliminaires, 
N = (gt — 2p) v 4 4- (o* — o — 2 erp 4- 4 P 4- p : ) a 4 
4- (•/— 3p 4 )2[(r?a 4- +• cy 4-.. )»] 

— (v 4 — 3/x 4 ) [p — 2 ("* 4~ *P 4- ry 4-. . J- 



- 
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Donc l'erreur moyenne à craindre, lorsqu'on prendra 




sera 

V^^-^t-Sf^ + ^^H- ..)■} 

40. 

La quantité 

qui entre dans l'expression précédente, ne peut pas géné- 
ralement se réduire à une forme plus simple. Cependant 
on peut assigner deux limites entre lesquelles sa valeur doit 
nécessairement être comprise : 

i°. On déduit facilement des relations précédentes, 

(fla + ^4-f 7 -.. + (aat? + £ p' -h cy' H ) J 

-4. («a" -h M" -h <*v" ■+-• . + • • • = * a H" h ? + c 7 -H - • -ï 

* 

d'où nous concluons que 

«a + -f q +. .. 

est une quantité positive plus petite que l'unité, ou du 
moins qu'elle n'est pas plus grande. Il en sera de même de 
la quantité 

nV + A'f -+- r' V ' 
qui est égale à la somme 

-H + ^'p"-+- c'y" 4-... )'-+-••• i 

de même 

5. 
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sera plus petit que l'unité ; et ainsi de suite. Donc 

2[(«« + *pH-* 7 + •••)'] 

est nécessairement plus petit que n. 
a°. On a 

V(rta-hip + C7-4-...) = P» 

car 

d'où Von déduit facilement que 

est plus grand, ou du moins n'est pas plus petit, que 



Par conséquent le terme 
est nécessairement compris entre les limites 



v* — 3 v 4 — 3 p 

— C_ et — —7 

ou bien, entre les limites plus étendues, 

v «_3»« v » —3^ 
i— ~et — • 

cr — p nj — p 

Donc le carré de l'erreur moyenne à craindre pour la va- 
leur 

M 

r cr — p 



est compris entre les limites 



av < - 4f* < et J±L_. 
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de sorte qu'on pourra atteindre un degré de précision aussi 
grand que l'on voudra, pourvu que le nombre des observa- 
tions soit suffisamment grand. 

Il est très-remarquable que dans l'hypothèse de l'art. 9 
(III) , sur laquelle nous nous étions autrefois appuyé pour 
établir la théorie des moindres carrés, le second ternie 
du carré de l'erreur moyenne disparaît complètement (car 
on a v* — 3|m* = o)$ et comme, pour trouver la valeur 
approchée jx, de l'erreur moyenne des observations, il faut, 
dans tous les cas, traiter la somme 

comme si elle était égale à la somme des carrés des w — p 
erreurs fortuites, il en résulte que, dans cette hypothèse , 
la précision de cette détermination devient égale à celle que 
nous avons trouvée, art. lo, pour la détermination déduite 
de gj — p erreurs vraies. 
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SUPPLÉMENT 

A LA 

THÉORIE DE Li COMBINAISON DES OBSERVATIONS 

QUI EXPOSE AUX MOINDRES ERREURS. 



PRÉSENTÉ LE iG SfcPTESMRE l8>f> A LA SOCIÉTÉ ROYALE DE I.OTTISÛIE. 



Dans le Mémoire précédent, on a supposé que les quan- 
tités à déterminer, à l'aide d'observations imparfaites, 
fussent dépendantes de certains éléments inconnus, en fonc- 
tion desquels on sut les exprimer : le problème consistait 
alors à déduire des observations, le plus exactement pos- . 
siblc, la valeur de ces éléments. 

Dans la plupart des cas, c'est effectivement ainsi que la 
question se présente; mais quelquefois il en est un peu 
autrement, et l'on pourrait même se demander, au premier 
abord, si le problème se ramène au précédent. 11 n'est p;is 
rare, en effet, que les quantités auxquelles se rapportent 
les observations ne soient pas explicitement exprimées en 
fonction de certains éléments, et qu'elles ne paraissent ré- 
ductibles à une telle forme que par des opérations difficiles 
ou ambiguës. Il arrive souvent, d'autre part, que la nature 
du problème fournit certaines conditions auxquelles les 
valeurs observées doivent rigoureusement satisfaire. 

Cependant, en y regardant de plus près, on aperçoit que 
ce cas ne diffère pas essentiellement du précédent , et qu'il 
peut s'y ramener. Si l'on nomme, en effet, le nombre des 
grandeurs observées, cl ? celui des équations de condition , 
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en choisissant, parmi les premières, o — <j quantités, rien 
ne nous empêche de les considérer comme nos seules in- 
connues, les autres étant regardées comme des fonctions de 
celles-là, que les équations de condition définissent. Par cet 
artifice, nous rentrerons dans le cas du Mémoire précé- 
dent. 

Néanmoins, quoique cette marche conduise souvent au 
résultat d'une manière assez commode, on ne peut nier 
quelle ne soit moins naturelle, et il est, par consé- 
quent, désirable de traiter le problème sous l'autre forme, 
qui admet, du reste, une solution très- élégante. Il y a 
plus : cette solution nouvelle, conduisant à des calculs plus 
rapides que la précédente toutes les fois que <j est moindre 

que - or , ou, ce qui revient au même, toutes les fois que le 

nombre des éléments que nous avons désigné par p dans le 

Mémoire précédent, est plus grand que il faudra, dans 

ce ras, préférer la solution nouvelle, quand bien même il 
serait facile, par la nature du problème, de faire dispa- 
raître sans ambiguïté les équations de condition. 

2. 

Désignons par y, v' , i/'/etc, les quantités, en nombre tar, 
dont les valeurs sont fournies par l'observation. Supposons 
qu'une inconnue dépende de celles-là et soit exprimée par 
une fonction connue m, de y, \>* ', v h r , etc. Soient /, /', etc., 
i f que dev iennent les quotients différentiels 

du du du 

— « — , ..... 

dv dv f ,/.." 

lorsqu'on y substitue à (/, v" etc., leurs valeurs véri- 
tables. Si l'on substituait à i>, v* , etc., dans la fonc- 
tion n , leurs valeurs véritables, on obtiendrait aussi la 
véritable valeur de u \ mais si les observations sont afl'eo- 
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tées d'erreurs e', e", etc., il en résultera pour u une 
erreur totale représentée par 

le -f- /V4-/V-+-..., 

pourvu que Ton puisse, comme nous le supposerons tou- 
jours (lorsque u n'est pas linéaire) , négliger les carrés et 
les produits des erreurs e, e', e", etc. 

Quoique la grandeur des erreurs e, e', e", etc., soit in- 
certaine , l'incertitude attachée à la valeur trouvée pour u . 
peut généralement se mesurer par Terreur moyenne à 
craindre dans la détermination adoptée. D'après les prin- 
cipes développés dans le premier Mémoire, celte erreur 
moyenne est 

m y m', m' etc., étant les erreurs moyennes des diverses 
observations. Si toutes les observations sont affectées du 
même degré d'incertitude, cette expression devient 

Il est clair d'ailleurs qu'au degré d'approximation auquel 
nous nous arrêtons, on peut remplacer/, /', /", etc., par 
les valeurs que prennent les coefficients différentiels 

tin du dit 

lorsqu'on y remplace v' y v n \ etc., par leurs valeurs ob- 
servées. 

3. . 

Lorsque les quantités y, f', v" , etc., sont indépendantes, 
l'inconnue ne peut se déterminer que d'une seule manière, 
et l'incertitude attachée au résultat ne peut être ni évitée 
ni diminuée. Les observations fournissent une valeur de 
l'inconnue qui n'a, dans ce cas, rien d'arbitraire. 11 en est 
tout autrement lorsque les quantités i'', v", etc., sont 
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liées par des relations nécessaires , que nous supposerons 
exprimées par a équations, 

X = o, Y = o", Z = o , . . . , 

X,Y,Z,elc, désignant des fonctions données des indétermi- 
nées t>, y', v" , etc.; car, à la fonction m, on peut, dans ce cas, 
substituer toute autre fonction U telle, que la différence 
U — u s'évanouisse identiquement en vertu des équations 

X = o, Y = o, Z — o , . . . . 

Si les observations étaient rigoureusement exactes, cette 
substitution ne changerait en rien le résultat; mais, en rai- 
son des erreurs inévitables , à chaque forme adoptée pour u 
correspondra un résultat di lièrent, et Terreur commise, au 
Heu d'être 

l e+ i' e ' /V'-h..., 

deviendra 

Le -f- LV -+ LV +. . 

en désignant par L, L', L", etc., les quotients dilléren- 
tiels 

rfU <TU iiV 
~<U> ' ih' ' ds>"" 

Quoiqu'il soit impossible d'assigner la valeur des diverses 
erreurs , nous pouvons cependant comparer les erreurs 
moyennes à craindre dans les diverses combinaisons. La 
combinaison la plus avantageuse sera celle qui donnera, à 
Terreur moyenne, la valeur minimum. Cette erreur étant 

yJUm* t L"/«' 2 L'"///"' + . . . , 
nous devons chercher à rendre la somme 

L J //< 3 H- L"///- -+ L" W* -+-..., 
aussi petite que possible. 



( 74 ) 
A. 

Les fondions U, en nombre infini, par lesquelles on peul 
remplacer m, ne différeront les unes des autres, dans nos 
recherches, que par les valeurs qu'elles fourniront pour 
L , L', L", etc. : nous devons donc, avant tout, chercher 
les relations qui existent entre les systèmes de valeurs que 
peuvent prendre ces coefficients. Désignons par 

,_ » // 

b, b', b\ . ., 



les valeurs que prennent les coefficients 



</x 


</X 






Ta/' 




dY 




dY 


ITv' 






dZ 


r/Z 


d'I 




57' 





si l'on y substitue pour i>, ^ ', y", etc., leurs valeurs véri- 
tables. Il est clair que si Ton donne à r, \>\ v", etc., des 
accroissements dv, dv' , dv" , etc., qui ne changent pas 
X, Y, Z, etc., et leur laissent, par conséquent, la valeur 
zéro, ces accroissements, qui satisferont aux équations 

o = adv -f- a'dv' 4- a"dv" H- . . . , 
o = bdv 4- b'dv' ~h Vdv" 4- . . . , 
o = c dv 4- c'dv' 4- c"dt>"-h. . . , 



ne changeront rien à la valeur de U — m, et Ton aura , par 
conséquent, 

o = (/ — L) dv 4- [t — V) dv' 4- [I" — V) dv" 4- 

On en conclut facilement que L , L , L", etc., doivent avoir 
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la forme 

L = / -h ax -f- b y -f- cz -f- . . . , 
L' = ï -\- a'x -f- b'y + </» + ..., 
L"= /"-h- 6" 7 -f- c"z 4- . . . , 

• 

or, j, 2, etc., désignant des multiplicateurs déterminés. Réci- 
proquement, il est clair que, pour toutes les valeurs de x, 
y, s, etc., on pourra former une fonction U, pour laquelle 
les valeurs L, L', L", etc., seront précisément celles que 
fournissent ces équations , et cette fonction pourra , d'après 
ce qui précède, être substituée à u. La forme la plus simple . 
qu'on puisse lui donner est 

U = « -h j?X+/Y+ zZ . .4- u' , 

m' désignant une fonction de y, v" , etc., qui s'annule 
identiquement lorsque X, Y, Z, etc., sont nuls, et dont la 
valeur, dans le cas actuel, sera maximum ou minimum (car 
ses dérivées devront toutes s'annuler). Mais Y introduction 
de celte fonction n'apporte aucune différence dans les ré- 
sultats. 

5. 

Il est maintenant facile d'assigner àx,y,r, etc., des 
valeurs telles, que la somme 

L m 7 4- L' 2 m' 7 -h L" m" 7 -+-... 
soit un minimum. 

Il est clair que pour atteindre ce but, la connaissance des 
erreurs moyennes absolues wi, m', m", etc., n'est pas néces- 
saire; il suffit de connaître leurs rapports. Introduisons, en 
effet, au lieu de ces quantités, les poids des observations , p, 
p',p", etc., c'est-à-dire des nombres réciproquement pro- 
portionnels à m'*, m"*, etc. Les quantités x, y, z, etc., 
devront être déterminées de telle sorte que le polynôme 

. . , 

P P 
acquière une valeur minimum. Supposons que a. 0 , y", 
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s°, etc., soient les valeurs déterminées de x y j^ z, etc., 
auxquelles répond ce minimum, et adoptons les notations 
suivantes : 



à 1 












•= (*«), 


— 


4- 


7 

P- 


4- 




4-.. 


P 






P 




au 




n' h' 




n" h" 
u u 




.= («*), 


-h 




4- 




-h.. 




p' 




i) 
r 






P 




ac 




a'c' 




a c 




• = , 




4- 


P 


4- 


•i 


4-. . 


P 






P 








b" 




b'* 




.= (**), 


— 

p 


4- 


r 

P 


4- 


P 


+ . . 


c< 












.= (cc), 




-h 




4- 




4-. . 


p 




P' 




P" 




al 




a'V 




a"l" 




.= («/), 




4- 


* 

P 


4- 




4-. • 


P 






P" 




bl 




b'r 




b"l" 




. = (*/), 




4- 




4- 




4-. . 


7 




p' 




P" 




cl 




c'I' 




c"ï' 




.= (c/), 




4- 


y 


4- 




4-.. 


p 






V 





La condition de minimum exige évidemment que Ton 
ait 

o — (aa)x ù 4- (ab) y" 4- (ac) s° 4-. . . 4- («/), 
o " (ba) x° j°4- (Ac) z° 4- . . .4- (A/), 

o = (ca) x* 4- (c£)/ u 4- [ce) 2° 4- ... 4- (r/), 

Après que ces équations auront fourni les valeurs de a? 0 , y 0 , 
z°, etc., on posera 

a x° -\- b y 0 4- c z" 4- . . . 4- l = L, 
a' x° 4- b' y* 4- C s 0 4- . . . 4- t = V , 
a»* -h b"y« 4- *V 4- . . . 4- 1"~ L", 

et la fonction la plus propre à déterminer notre inconnue, 
à laquelle correspond Terreur moyenne minimum, sera 
celle dont les coefficients différentiels, pour les valeurs con- 
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sidérées des variables, seront égaux à L, L', L", etc. Le poids 
de la détermination ainsi obtenue, sera 

(3) P = 



L» L" I,"' 

H j ~ f~ — J, ~\~ • • • 

P P l> 



c'est-à-dire que ~ sera précisément la valeur que prend le 

polynôme considéré plus haut, pour les valeurs des variables 
x,f, z,etc, qui satisfont aux équations (i). 

6. 

Dans l'article précédent, nous avons montré à déterminer 
la fonction U qui fournit la détermination la plus conve- 
nable- de l'inconnue u. Examinons actuellement quelle 
est la vafeurqui en résulte Désignons-la par K: on l'ob- 
tiendra en substituant, dans U, les valeurs observées des 
quantités f, i^', v" , etc. Soit A la valeur que prend u lors- 
qu'on y fait les mêmes substitutions, soit enfin y. la véri- 
table valeur de cette inconnue, telle qu'on l'obtiendrait 
par la substitution des valeurs véritables de y, v f , v ff , etc., 
soit dans m, soit dans U. On aura 

/ = x + /V -h te" . . 

K=« + Le + LV+LV..., 

et , par conséquent , 

K = k -h (L - /) e -h ( 1/ - /' ) J + ( V - 1" ) > -f. . . . . 

Substituant, dans cette équation, à la place de L — /, 
V — L" — etc., leurs valeurs fournies par (2), et 
posant 

' ae -H a'e -h a" e" 4-. . . = A, 
■ A 1 k ce -h c'e' + c"e"+. . .= G, 



nous aurons 

(5) K = k-hAx* + '&y 9 + Qz 1> -h 
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11 n'est pas possible de calculer <l , £, elc, par le moyen 
des formules (4), car les erreurs <?,e',e", etc., qui y figurent, 
ont des valeurs inconnues , niais on voit facilement que ces 
quantités JW, ift», C, etc., ne sont autre chose que les valeurs 
deX, Y, Z, etc., qui correspondent aux valeurs observées 
de t>, y', v", etc., et alors le système des équations (î), (3), 
(5), forme la solution complète de notre problème. Il est 
clair, en effet que Ton peut appliquer au calcul de a, a', 
a'\ . . . , b , b f , b" , etc., la remarque faite à la fin de l'ar- 
ticle 2, à l'occasion des quantités /, l'\ etc., c'est-à- 
dire remplacer les valeurs véritables de y, v* , v ,r y etc., 
par les valeurs observées. 

7. 

On peut substituer à la formule (3), qui représente le 
poids de la détermination la plus probable, plusieurs ex- 
pressions qu'il est utile d'indiquer-, remarquons, d'abord, 
qu'en ajoutant les équations (2) après les avoir multipliées 



a a a 



par-î — -77Î etc., on aura 
1 P P P 

. . . . , aL a'V a"L" 

laa) x» + (ab)f -+- (ac) z> -h. . . = + — -h -y- + 

Le premier membre est nul *, en désignant donc, d'après 
la notation adoptée, le second membre par (aL), on aura 

(aL) = o, 

et de même 

(6L) = o, (cL) — o,.... 

L L' L" 

Multiplions ensuite les équations (2) par -1 etc. \ 

nous aurons, en les ajoutant, 

IL l'V V'L" L 1 L" L"' 

— H ; 1 j, h — - I 7 I 77 H" • • • » 

P P P P P P 

et , par suite , nous avons cette seconde expression du poids, 

1 



P — 



IL l'V l L 

1 "H 7T 

/' P V 
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Si, enfin, nous ajoutons 1rs mêmes équations (2) après 

/ i< i" 

les avoir multipliées par - , — , - , cte. , nous obtiendrons 
une troisième expression du poids , 



P = T~ 



(«/)*•+ {bl) y* -h {cl)z> + . .. + (//)' 
où l'on a posé, conformément à la notation précédente, 

On passera facilement de là à la quatrième expression 
du poids , 

i = (//) — (««) x"— (bb)y»- (fc)a" — ... 

8. 

La solution générale que nous venons d'exposer s'ap- 
plique principalement au cas où Ton n'a qu'une seule in- 
connue à déterminer. Lorsqu'on cherche, au contraire, les 
valeurs les plus plausibles de plusieurs inconnues, dépen- 
dant des mêmes observations, ou lorsqu'on ignore quelles 
sont les inconnues qu'il faut, de préférence, déduire des 
observations, il convient de procéder d'une manière diffé- 
rente, dont nous allons actuellement nous occuper. 

Considérons x, y, z , etc., comme des indéterminées, et. 
posons 

/ (aa)x -h [ab) y -+- (ac) z -+- . . . = £, 

(6) \(ab)x + (bb)y+{bc)z+ ..= «, 
[ ' ) (ac) x-t-(bc) y -4-(«r)zH-...= Ç, 

• I 

Supposons que l'on en déduise, par l'élimination, 

(aa) Ç H- (a£) »j -H ( ay) Ç 4- . . . = jr, 

(7) Wpa)Ç-h(P?)»4-(p 7 )^+-.- = r, 
4 (ya) Ç -H (7p) » -h (77) Ç -h. . =2; 
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et remarquons, avant tout, que les coefficients placés symé- 
triquement sont nécessairement égaux, c'est-à-dire que 

(fSa) = («?), 
('/*) = ( a ï)> 
(7P) = (P7). 

comme cela résulte de la théorie de l'élimination et comme 
d'ailleurs nous le démontrons plus loin ; nous aurons 

x° ~ — (aa) {al) — («?)(/>/) — (ay)^/) - . . . , 

, 8 x b- = -(P*)(«0-(??)(W)-(p7)( f O — .., 

1 j S s« = ~(7*)(«0-(7p)(W)-(77)(^)-.-.* 

et en posant 

(*a) t k -h (a£) 1)1,4- (ay) C -h. . . = A, 

(j5«) -h(^) ifl»-r-(?y) 3 -h. . .= H, 



(9) 



(ya) X -f- (y t 5) Ui> -+- (yy) G -}-...= C, 



nous obtiendrons 

K = X — À(«/) — B(W) — C(c/)— 

et si nous posons de plus , 

a A -f- 6 B + c C + ...=/> «, 
rt' A + B 4-c'C 4-. . . = /?V, 
<?"A 4- *"B4-*"C -h...=r//V\ 



(,o) 



il viendra 

(n) K = A-/i — /V-/'V 

9. 



La comparaison des équations (7) et (9) nous apprend 
que les quantités auxiliaires A , B, C, etc., sont les valeurs 
que prennent les indéterminées .r,^, z , etc., lorsque Ton 
suppose 

£ = A, , >? — tf> , £ = £,... 5 
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on i'ti conclut évidemment, 

(an) A -h (ab) B+ (ne) C+... = «l, 
\ba) A -h (bb)ft-h(bc)C -h... = i&, 
(en) A -h (c£)BH-(cc)C 



(.*) 



En ajoutant les équations (10) après les avoir multipliées 



a a a 



par — » — » etc., on obtient 
r P P P 

X = aî + a's' +ûV + . • • , 

et de même 

i ^ = i î + è' £ ' + /,'V' + ..., 
G = c s -h c'e' -+- c'V -4- . . . , 



X étant, comme on Ta dit, la valeur que prend X lorsqu'on 
substitue à y, v\ v", etc., leurs valeurs observées; on aper- 
çoit facilement que si Ton applique à ebacune de ces quan- 
tités les corrections — e, — e', — e", etc., la valeur de X de- 
vient égale à zéro, et que de même Y, Z, etc., s'évanouis- 
sent par cette hypothèse. L'équation (i i) prouve aussi que K 
est la valeur que prend u par suite des mêmes substitutions. 

En nommant compensation des observations, l'applica- 
tion des corrections — e, — c', — e", etc. , aux grandeurs 
directement observées, nous sommes conduit à une con- 
séquence très-importante : 

Les observations compensées comme nous l'avons indi- 
qué, satisfont exactement à toutes les équations de condi- 
tion, et font prendre à toute fonction des quantités obser- 
vées la valeur qui résulterait de la combinaison la plus 
convenable des observations non modifiées-, et puisque les 
équations de condition sont trop peu nombreuses pour 
qu'on puisse en déduire les valeurs exactes des erreurs, 
nous aurons au moins trouvé, par ce qui précède, les er- 
reurs les plus plausibles. C'est sous ce nom que les quan- 
tités e, e', e", etc., seront désormais désignées. 

c. 6 



Digitized by Google 



( 8, ) 

Le nombre des observations étant plus grand que celui 
des équations de condition, outre le système des correc- 
tions les plus plausibles, on peut en trouver un nombre 
infini qui rendent les équations de condition exactement 
satisfaites. 

Il importe d'examiner les relations qui lient entre eux 
ces divers systèmes. Soit — E, — E', — E", etc., un pareil 
système de corrections , autre que le système le plus plau- 
sible ; nous aurons 

flE + a'E' ■+- a"E"H-, . . = X , 
bE-hb'E' + £"E"H-...= ift,, 
cE-hc'E'-f-c"E"-f-...= e, 



Multipliant ces équations par A , B, C, etc., et ajoutant, il 
vient, en ayant égard aux équations (10), 

peE 4-pVE' -4-p'V'E" + ...= A Jt. -f- Bift> CG-h .... 

JVlais les équations (i3) combinées de la même manière don- 
nent 

de la combinaison de ces résultats on déduit facilement 
p E' -4- p' E" 4- p" E'" H- ...=/>«' + p' t" + p" g"' -f-... 

et, par suite, la somme 

p E» -{-p' E' 1 -h //' E"- -h . . . 

est, nécessairement, plus grande que 

ce que Ton peut énoncer de la manière suivante : 

Théorème. — Les carrés des corrections qui peuvent 
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concilier les observations avec les équations de condition 
étant respectivement multipliés par les poids des observa- 
tions correspondantes , donnent une somme minimum 
quand on adopte les corrections les plus plausibles. 

On reconnaît précisément le principe des moindres carrés 
dont, au reste, les équations (12) et (10) peuvent facilement 
se déduire. La somme minimum, que nous désignerons 
désormais par S, est égale, d'après l'équation (i4)î à 

AX+B\ii, + Ce + 

■ 

11. 

La détermination des erreurs les plus plausibles étant in- 
dépendante de /, l','l f \ etc., fournit la préparation la plus 
commode, quel que soit l'usage que Ton veuille faire ulté- 
rieurement des observations. En outre, on voit sans peine 
que, pour atteindre ce but, il n'est pas nécessaire d'effec- 
tuer l'élimination indéfinie y c'est-à-dire de calculer (aa), 
(a/3), etc., il suffît de déduire des équations (12), par une 
élimination définie, les quantités auxiliaires A, B, C, etc., 
que nous nommerons, dans ce qui va suivre, les corrélatifs 
des équations de condition 

X = o, Y = o, Z = o, ..; 

ces quantités seront ensuite substituées dans l'équation (10). 

Cette mélhode ne laisse rien à désirer lorsqu'on demande 
seulement les valeurs les plus plausibles des quantités four- 
nies par l'observation. Mais il en est autrement lorsqu'on 
* désire, en outre, le poids de chacune des valeurs trou- 
vées. Quelle que soit alors celle des quatre formules précé- 
dentes que l'on veuille employer, il est indispensable de 
connaître L, I/, L", etc., ou, ce qui revient au même, 
x°, y 0 , z°, etc. ; par cette raison, il sera utile d'étudier, de 
plus près, l'élimination qui fournit ces quantités et d'obtenir 
une méthode plus commode pour la détermination des poids 

6. 
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12. 

Les relations qui lient les quantités dont nous nous oc- 
cupons sont notablement simplifiées par la considération 
de la fonction indéfinie du second degré 

(aa)x i H- 2 (ab)xy -+- 2 (ac) xz . . . 

-h (bb)y 7 -h i(bc)yz -h . . . -+- ( ce) z* -h . . . , 

que nous désignerons par T. 

Cette fonction est évidemment égale à 

(ax + bjr + cz+...y {„' X -\-b'y->re'z+...y • 
(15) - t -, H- ... 

De plus, on a évidemment 

( 16) T rr \x m y -h 4- . . . ? 

et si , enfin , on exprime x , ^, z , etc. , au moyen des équa- 
tions (7) , en fonction de £, ry , f , etc., on aura 

T = + a (ap)&r-h 2(a 7 ) ÇÇ -f- . . . + (pp) 

-4-a(p7)«Ç-f-...+ ( 77 )Ç>-|-.... 

La théorie développée plus haut fournit deux systèmes 
de valeurs déterminées pour les quantités x,y,z, etc., 
I , >? , Ç , etc. Le premier est 

% = -(«/), «=-(*/), ? = ..; 
à ce système correspond la valeur 

T = (")-p, 

ainsi qu'on le voit en comparant à l'équation (16), la troi- 
sième forme du poids P, ou par la considération directe de 
la forme (4). 

Le second système de valeurs est 

x = A, y — B, z = C,..., 
| = A , m = i)l , Ç = C , . . . ; 

la valeur correspondante de T est 

T = S, 
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comme cela est évident par les formules (10) et (i5), et 
encore parles formules (14) et (16). 



13. 

Nous devons, avant tout, faire subir à la fonction T 
une transformation semblable à celle qui a été indiquée 
(Theoria Motus , art. 182), et, avec plus de développe- 
ments , dans les Recherches sur Pallas. 

Posons , à cet effet , 

(if , l) = (6c) _(^Hf£], 

■ ' (aa) 

(bd, 1)_(W) ^y-, 



(ce, 2) — [ce ) — — yj- j 

(aa) (bb f 1) 

(i' si ='(,«)'- ^ - ^4* 1' £4, 

v y (on) (bb,i) (cc t n) 

et, ensuite (*) , 

{bbj 1)/ + (0**, 1) z -+- (M, i)tv-h. . . = 

(ce, 2) 3 H- (crf, 2) (v . .= Ç", 

3) * + = 



on aura 

n' 1 «p w * 
et y}', Ç", <p'", etc., se déduiront de £, y; , f . 9, etc., par les 

(* ) Dans les calculs précédents il suffisait de trois lettres de chaque série 
pour faire apercevoir la loi des formules ; il a paru nécessaire d'en faire ici 
figurer une quatrième, pour rendre l'algorithme plus manifeste. 

( Note de M. Gaiss. ^ 
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équations suivantes : 



K = r, - 



(an) 

_ , (" c ) r l > bc > „' 
' ~ (™) ' (l'b, ') ' 

t - = ¥ _ (fl) 5 _ (Mil) „' _ 



1 



(cd, 2} 

(aa) * (bb y i) " (ce, i) 



et l'on en tirera facilement toutes les formules utiles à notre 
but. Ainsi , pour déterminer les corrélatifs A , B, C, etc., 
nous poserons 

= i)S> — } X , 

( «« ) 



(.8) 



( ri ri ) ( M, I 

(aa) (bb 9 i) {ce, 2) 



et enfin A, B, C, D, etc., s'obtiendront par les formules 
suivantes, en commençant par la dernière : 

(aa) A 4- (ab) B-H (ac) C-h (ad) f) -h . . . = ^ , 

(^,i)b+(/^i)c + i ) r ) -+- . . . = uv, 

( l 9) [ (ce, 7.) C-h (cri, 2 )1)+ .. = ©", 

(r/r/, 3) !)-!-...:=<&'", 



Pour exprimer S, nous aurons la formule nouvelle 

« 

, . c «v a e" 1 , 

(20) b = — — -h -p. — - -1- — — - -+- 



(aa) (bb,i) (ce, 2.) [dd, A) 

enfin , le poids P, qu'il faut attribuer à la détermination la 
plus plausible de la quantité //, sera donné par la formule 



[aiy (bi 9 jy (r/,2)- - (di,3)> 

(aa) \bl,i) (rr,?.) (r7r/\3)~ 
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dans celte formule, on fait 

(««) (ii, l) («, 2) 

Les formules (i 7), (ai), dont la simplicité ne laisse rien 
à désirer, fournissent la solution complète de notre pro- 
blème. 

44. 

. Après avoir résolu le problème que nous avions en vue, 
nous allons aborder quelques questions secondaires qui 
éclaireront davantage l'ensemble de cette théorie. 

Nous chercherons, en premier lieu , s'il peut arriver que 
l'élimination qui fournit y y -z, etc., en fonction de J, 
y} y Ç, etc., devienne, dans certains cas, impossible*. Cela 
aurait évidemment lieu si les fonctions £, y;, f , etc., n'é- 
taient pas indépendantes les unes des autres. Supposons, 
pour un instant, qu'il en soit ainsi, et que l'une d'elles 
puisse s'exprimer en fonction des autres, de telle sorte que 
l'on ait la relation identique 

aÇ H- pu -h 7$ -h . . . = o , 

a, (3, 7, etc., désignant des nombres déterminés. 
On aura alors 

a (an) -h p (ab) -h 7 (ac) -f- . . . = o ; 
ot(ab) 4- $(bb) -h f{bc) . . = o, 
u(ac) -h bc) -h y (ce) -4-, . . = o, 

et si nous posons 

art -h p& +yc +.. 

art' + p6' -4- 7c' -t~. . . =p'e\ 



Digitized by Google 



( «8 ) 

on en déduira 

aG 4- a'®' 4- a" Ç>" 4- . . = o , 

Oe 4- b'e'+ b"%" -+-.. =o, 

«0 -h c'0' 4- c"9" -h. . . = O, 



et, par suite, 

/jO' H- //©'» 4- />"0" J 4- • • = o; 

/>, p', />", etc., étant, par leur nature, tous positifs, cette 
équation exige. 

e = o , e' = o , 0" = o , . . . . 

Si nous considérons les différentielles complètes dX, 
dY, d7j, etc., répondant aux valeurs de y, v\ v'\ etc., 
immédiatement fournies par les observations, ces dilîcren- 
ti elles 

adv 4- «' dv' 4- a" dv" 4- . . . , 
^4- ^'r//4- £"rfp"4-..., 
cdv 4- c' rf/ 4- c" dv" 4- • • . , 



d'après les résultats précédents, seront liées, les unes aux 
autres, de telle sorte qu'en les multipliant respectivement 
par «, (3,y, etc., la somme des produits sera identiquement 
mille, en sorte que, parmi les équations 

X ~ o , Y = o, Z =. o . . , 

il eu est une au moins que l'on peut regarder comme 
inutile , car elle sera satisfaite dès que les autres le se- 
ront. 

En examinant la question de plus près, on voit que 
cette conclusion n'est applicable qu'à des valeurs des 
variables infiniment peu différentes de celles que fournit 
l'observation. Il y a , en effet, deux cas à distinguer : le 
premier est celui où l'une des équations 

X = o, Y = o, Z — o 7 . . . , 
est renfermée dans les autres d une manière générale cl 
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absolue, et peut, par conséquent, être supprimée; le se- 
cond est celui où , pour les valeurs particulières de y, 
v' ', \J\ etc., auxquelles se rapportent les observations, Tune 
des fonctions X, Y, Z, etc., X par exemple, acquiert une 
valeur maximum ou minimum, ou, plus généralement, une 
valeur dont la différentielle s'annule lorsque les autres 
équations restent satisfaites. 

Mais comme nous ne considérons pour nos variables que 
des variations dont les carrés soient négligeables, ce second 
cas (qui dans les applications ne se présentera que bien rare- 
ment) , pourra être assimilé au premier, et l'une des équations 
de condition pourra être supprimée comme surabondante. 

Si les équations restantes sont indépendantes dans le sens 
que nous venons d'indiquer, on peut, d'après ce qui pré- 
cède, être certain que l'élimination est possible. Nous nous 
réservons, du reste, de revenir sur cette matière qui mé- 
rite d'être examinée comme subtilité théorique plutôt que 
comme question d'une utilité pratique. 

Dans le premier Mémoire, art. 37 et suivants, nous 
avons montré le moyen de fixer, à posteriori, d'une ma- 
nière très-approchée, le poids d'une détermination. Si les 
valeurs approchées de xs quantités sont fournies par des 
observations également précises, et qu'on les compare avec 
les valeurs qui résultent pour elles des hypothèses les plus 
plausibles qu'on puisse faire sur les p éléments dont elles 
dépendent, on a vu qu'il fallait ajouter les carrés des diffé- 
rences obtenues, diviser la somme par gj — /o, et que le 
quotient pouvait être regardé comme une valeur approchée 
du carré de Terreur moyenne inhérente à ce genre d'ob- 
servations. 

Si les observations sont inégalement précises, la seule 
modification que Ton doive apporter aux préceptes précé- 
dents, consiste en ce que Ton doit multiplier les carrés des 



d byJttoogle ! 

4 % 



(9") 

différences par les poids respectifs des observations corres- 
pondantes, et Terreur moyenne obtenue de cette manière 
se rapporte aux observations dont le poids est pris pour 
unité. 

Dans le cas actuel , la somme des carrés dont nous par- 
lons se confond évidemment avec la somme S, et la diffé- 
rence u — p avec le nombre a des équations de condition. 
Par suite, pour Terreur moyenne des observations dont le 

poids est i, nous aurons l'expression y/^, et la détermi- 
nation sera d'autant plus digne de confiance que a sera plus 
considérable. 

Mais il e3t bon d'établir ce résultat, indépendamment des 
raisonnements du premier Mémoire; pour y parvenir, nous 
introduirons quelques notations nouvelles. Supposons 
qu'aux valeurs 

Ç = <* , r> = b , Ç = c,. . . , 

répondent 

de sorte que l'on ait 

a = a (aa) b (ap) + c (ay )-+-... , 

P = "(«P) + *(W)H-'(P7) 
7 = «(a 7 ) + b{fa)-hc[- n ) + 

et, en outre, qu aux valeurs 

£ = a', y) = b' y Ç = c , . . . , 

répondent 

x = a', Y — p', 3 = 7%...; 
enfin, qu'aux valeurs 

l = a", »i = Ç = c" , . . . , 

répondent 

j: = a", j = P", z = 7 ", . . . , 
et ainsi de suite. 



• > 
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La combinaison des équations (4) et (9) fournit 

A = «t?+ a!é -h ct."e" -h . . . , 
B=p*+pV-+-p'V-H..., 
C — ye 4-7V + y" e" 4-..., 

et, comme on a 

S= A, A -h Db B -H e C H- . • • , 

on aura 

S=(fl«+«V + «"«"-4-.. •) (ac+«V 4- a'V-K ..) 
.+-(6*4- iV4-tV + ...) (pff -t-pV + p^^-f-...) 
_i_ ( C e -h cV 4- c'V-h. ..) (7e + 7V +7 / V /, H-. . .)-»-•••• 

16. 

La série des observations qui fournissent les quantités 
v, S, v", etc., affectées des erreurs fortuites e, e f , <?", etc., 
peut être considérée comme une épreuve qui ne fait pas 
connaître, il est vrai, la grandeur de chaque erreur, 
mais qui, par le moyen des règles exposées plus haut, per- 
met de déterminer la quantité S, fonction connue de toutes 
les erreurs. Dans une telle épreuve, les erreurs peuvent être 
les unes plus grandes, les autres plus petites; mais plus 
sera grand le nombre des erreurs employées, plus il y aura 
une grande probabilité que S diffère peu de sa valeur 
moyenne : la difficulté revient donc à trouver la moyenne 
de S. 

Par les principes exposés dans le premier Mémoire et 
qu il est inutile de reproduire ici , on trouve pour cette va- 
leur moyenne 

{a a 4- b p -h cy 4- . . . ) m- 4- (*' a -h V p' 4- c' y' 4- . . ) w" 4- . . . . 

En nommant l'erreur moyenne qui coi respond aux obser- 
vations dont le poids est 1 , de telle sorte que Ton ait 

= pm* = p'm" = p"m" 3 = . . . , 
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l'expression précédente peut s'écrire comme il suit : 



Mais on a trouvé 



«a a'*' a" ce" 



+ -7- 4-...= (an) (aa) + (ao) (af) 4- (ac) (a 7 ) 



/» /> 

or, le second membre est l'unité, comme on le reconnaît 
facilement par la comparaison des équations (6) et (7). 



On trouvera de même 



ip b'W b" S" 



c 7 c' 7 ' c" 7 " 



= 1, 



/' />' />" ' 



et ainsi de suite. 

D'après cela , la valeur moyenne de S devient <7jm*, et si 
1 on juge permis de regarder la valeur fortuite de S comme 
égale à la valeur moyenne, on en conclut 



17 

On peut apprécier la confiance que mérite cette détermi- 
nation en calculant Terreur moyenne à craindre, soit 
pour sa valeur propre, soit pour celle de son carré. La se- 
conde sera la racine carrée de la valeur moyenne de l'ex- 
pression 



(§->•)• 



dont le développement s'obtiendra par des raisonnements 
semblables à ceux qui ont été exposés dans le premier Mé- 
moire (art. 39 et suivants). Nous les supprimons pour abré- 
ger, en nous contentant d'indiquer le résultat. 
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L'erreur moyenne à craindre dans la détermination du 
carré s'exprime par 



si 



2f . + v 1 _3£ N 

<T CT 1 



v* étant la valeur moyenne des quatrièmes puissances des 
erreurs dont le poids est l'unité, et N la somme 

(«a+ b$ -hey -h.. (a' <*' + b'$'-hc'y' + . ..)' 

-t- (a'V -+- b"p 4- c"i' -h. ..)'+...• 

Cette somme ne peut pas, en général, se simplifier 5 mais, 
par un procédé analogue à celui dont on a fait usage au pa- 
ragraphe^ du premier Mémoire, on peut montrer que sa 

valeur est comprise entre xs et ^ Dans l'hypothèse sur la- 
quelle nous avions primitivement établi la méthode des 
moindres carrés , le terme qui contient cette somme dis- 
paraît à cause de 

et la précision que l'on doit attribuer à la détermination 



=v1 



est par conséquent la même que si l'on avait opéré sur u 
observations entachées d'erreurs exactement connues, con- 
formément aux préceptes des art. 15 et 10 du premier 
Mémoire. 

18. 

Pour la compensation des observations, il y a, comme 
nous l'avons dit, deux opérations à exécuter: première- 
ment, il faut déterminer les corrélatifs des équations de 
condition, c'est-à-dire A, B, C , etc., qui satisfont aux équa- 
tions (12) -, secondement, substituer ces quantités dans l'é- 
quation (10). La compensation ainsi obtenue peut être dite 
parfaite et complète, par opposition à la compcnsaiion 
imparfaite ou incomplète. Nous désignerons de celte der- 
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nière manière celles qui résultent des mêmes équations (10), 
dans lesquelles on substituera des valeurs de A, B, C qui ne 
satisferont pas aux équations (12)» c'est-à-dire qui satisfe- 
ront à quelques-unes seulement, ou à aucune. Nous ne nous 
occuperons pas ici d'un tel système de corrections, et nous 
ne leur accorderons même pas le nom de compensation. 

Lorsque les équations (10) sont satisfaites, les sys- 
tèmes (12) et (i3) deviennent équivalents, et la différence 
dont nous parlons peut alors s'énoncer comme il suit : Les 
observations complètement compensées satisfont aux équa- 
tions de condition 

X — o , Y = o, Z — o y . . . ; 

les observations incomplètement compensées ne satisfont 
qu'à une portion de ces équations, et peut-être à aucune; 
la compensation à la suite de laquelle toutes les équations 
sont satisfaites, est nécessairement complète. 

19. 

Il résulte de la définition même des compensations, que 
la réunion de deux systèmes de compensations peut en four- 
nir un troisième, et l'on voit qu'il importe peu que les 
règles données pour obtenir une compensation parfaite 
soient appliquées aux observations primitives ou aux ob- 
servations déjà imparfaitement compensées. 

Soient — 0, — 0', — 0", etc., un système de compen- 
sations incomplètes, résultant des formules 

/ 07? r= A 0 * -hWb +C°c -h .., 
) 0>' - A 9 a' 4- B° b' CV +. . ., 
1 ] J &"p" r= A 9 a" -4- B°6" -h CV -h . . , 

Les observations ainsi changées ne satisfaisant pas à toutes 
les équations de condition, soient iJ!>*, C*, etc., les va- 
leurs que prennent X, Y, Z, etc., quand on y substitue les 
valeurs ainsi obtenues pour y, v\ v" etc. On devra cher- 
cher les valeurs A*, B*, C*, etc., satisfaisant aux équa- 
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lions 

X* = A* (an) -h B* (ah) -h C* (ac) -f- . . . , 
i)l>* = A* (ni) + B* (M) -f- C* (bc) -h . . . , 
G* = A* (ac) B* (bc) -h C* (ce) H- . . . , 

et cela fait, la compensation complète des observations 
ainsi modifiées se fera par les nouveaux changements — x, ■ 
— x', — x", etc.; x, x', x", etc., se déduisant des for- 
mules 

/ v.p = A*rt 4- B*£ -f- C*c -h. . ., 
. J x>' = A*«'4-B*A' -f-tfV 

1 ' x"/ = A*a" + B*6" + C*c" +. 

Cherchons commentées corrections s'accordent avec la com- 
pensation complète des observations primitives. Il est clair 
d'abord que l'on a 

JL* = X -aiè-a'e'- a"0"-..,, 
«I* = dî, — Z>0 — £'0' — b"(d" . , 
G* = C — c0 — c'0' — c"0" — 

En substituant dans ces équations, pour 0, 0', 0", etc., 
leurs valeurs fournies par le système (I), pour V, ifc*, 
G*, etc., celles que donne le système (II), il vient 

X = (A°H-A*) (an) H- (B° -f- B*) (ab)-h..., 
.Ub = (A°-*-A*) (ab) 4-(B°-hB*) (bb)+..., 
G = ( A 0 H- A*) ( ac ) + (B° -h B*) ( bc) -h . . . ; 

et il suit de là que les corrélatifs des équations de condi- 
tion (i a) sont 

A = A 0 H~ A*, B = B° + B*, C = C° + C*, . . . , 

et alors les équations (to), (I) et (III) montrent que l'on a 

e=8+x, i' = 8' + /, e" = 0"+x",...; 

et, par suite, la compensation parfaite a la même valeur 
pour chaque inconnue, soit qu'on la calcule directement, 
soit qu'on l'obtienne médiatement en partant d'une coin-" 
pensation incomplète. 
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20. 

Lorsque les équations de condition sont très-nombreuses, 
la détermination des quantités corrélatives A, B, C, etc., 
peut exiger des calculs tellement longs, que le calculateur 
en soit rebuté-, il pourra être avantageux d'obtenir, dans 
ce cas, une compensation complète, à l'aide d'une série 
d'approximations reposant sur le théorème de l'article pré- 
cédent. On partagera, pour cela, les équations de condition 
en deux ou plusieurs groupes, et Ton cherchera d'abord une 
compensation qui rende satisfaites les équations du pre- 
mier groupe. On traitera ensuite les valeurs modifiées par 
ce premier calcul , et on les corrigera de nouveau en ayant 
égard seulement aux équations du second groupe. Ce se- 
cond calcul donnera des résultats qui , en général, ne satis- 
feront plus aux équations du premier groupe, et il faudra, 
si l'on n'a formé que deux groupes, revenir alors au premier 
et y satisfaire à l'aide de nouvelles corrections. Les observa- 
tions seront ensuite soumises à une quatrième compensa- 
tion , dans laquelle on n'aura égard qu'aux conditions du 
second groupe 5 et en opérant ainsi alternativement sur l'un 
et l'autre groupe d'équations, on formera des corrections 
qui seront nécessairement de plus en plus petites. Si le 
choix des groupes a été fait habilement, on arrivera bien 
vite à des valeurs que les corrections ultérieures ne chan- 
geront plus. 

Quand on forme plus de deux groupes, on doit procéder 
de la même manière, les divers groupes étant employés suc- 
cessivement jusqu'au dernier, après quoi ou revient au pre- 
mier pour les reprendre dans le même ordre. Il nous suffit 
d'avoir indiqué ce procédé, dont la réussite dépendra beau- 
coup de l'habileté du calculateur. 

21. 

Il nous reste à donner la démonstration du lemme admis 
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dans l'art. 8. Adoptons, pour plus de clarté, des notations 
plus propres à mettre la démonstration en lumière. 

Soient x f , etc., des indéterminées-, supposons 
que les équations 

n w x° H- n» x* -+- n<" x" -f- . . . = X°, 

« ,0 /« -h /!" X' ■+- «•» X" -f- . . . = X'. 
n njcO -f- x > 4_ „« x " _|_ , . . — X", 



aient donné, par élimination , les suivantes : 

N"X°-+- N 4 «X'-f- N"X"+ ...= *% 
IV'X 0 -f- N" X' -h N i: X"+...= x', 
N 20 X' -h N" X' 4- N M X" -h ... — x", 

En substituant, dans les deux premières équations du se- 
cond système, les valeurs de X°, X', X", etc., fournies par 
le premier, nous obtiendrons deux équations identiques : 

x° = N 00 (« 09 x° -+- x' ■+- x" -j- . . . ) 
H- N 01 («'•x° H- «'«x' H- n"x" H- . . ,) 
■4- N" («" x° H- « 2I x' -h x" -H • . . ) 
-h N" (/î 39 x° fl 31 x' -H « M H- . ■ . } 



x 7 N'° (n°°x° -h ««'x' -h /i«x' / -h. • •) 
-»- N" a' H- /*" x' -f- n 11 x" -f-, .) 
-h N' 1 (*» x° -+- /i» x' H- x" -h . . . j | . 
H- 

Ces équations étant identiques , on peut y substituer telles 
quantités que Ton voudra à la place de x 0 y x f i x", etc. Fai- 
sons dans la première 

x e =N'% a' = N", x" = N",..j 

et dans là seconde 

x« = N", x^N", x" =N°V..; 

En retranchantensuite les deux identités membre à membre, 
o. 7 
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il viendra 

N" — N 8 ' = (N M N" — N" jN°') (/*•' — «'•) 
+ (N»» N" — N' 9 N") (/!" — 
(N« N' 3 — N'» N» 3 ) («•» — /i M ) 



(n»' N" — N" N") (/i n — n") 
+ ( N o. N .» _ N" N M ) (* n — * 3 ') 


+ (N" N' 3 — N' 1 N") (/*" — 

-h 

Ce que Ton peut écrire plus brièvement de la manière sui- 
vante : 

or, (3 désignant deux indices pris au hasard } on conclut de 
là que les égalités 

et, généralement , 
entraînent nécessairement 

et, comme l'ordre des indéterminées est arbitraire, il est 
évident que, dans la supposition admise, on aura généra- 
lement 

22. 

La méthode exposée dans ce Mémoire devant surtout s'ap- 
pliquer utilement aux calculs de haute géodésie, le lecteur 
nous saura gré d'y joindre quelques exemples puisés dans 
cette partie de la science. 

Les équations de condition qui existent entre les angles 
d'un système de triangles, peuvent, en général, se par- 
tager en trois catégories. 

1. La somme des angles horizontaux , formés autour d'un 
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même sommet et embrassant la totalité de l'horizon , doit 
être égale à quatre droits. 

II. La somme des angles de chaque triangle peut toujours 
être regardée comme connue -, car, lors même que le triangle 
est situé sur une surface courbe, l'excès de la somme de ses 
angles sur deux droits peut être calculé avec une telle 
approximation, qu'il est permis de considérer le résultat 
comme absolument exact. 

III. Enfin , on obtient un troisième genre de relations en 
examinant les rapports des côtés dans les triangles qui for- 
ment un réseau fermé. Si, en effet, les triangles sont tel- 
lement placés, que le second triangle ait un côté a commun 
avec le premier, et un côté b commun avec le troisième; 
si le quatrième triangle a deux côtés c et d respectivement 
communs avec le troisième et le cinquième, et ainsi de 
suite, jusqu'au dernier triangle, qui ait avec le précédent 
un côté commun h , et avec le premier de tous un côté com- 
mun /, les quotients 

rc h c d l 
l a b c A 

pourront se calculer au moyen des angles qui leur sont 
opposés dans le triangle dont les deux côtés comparés font 
partie , et comme le produit de ces fractions est évidemment 
l'unité, on aura une relation entre les sinus des divers 
angles mesurés (diminués du tiers de l'excès sphériqueou 
sphéroïdique lorsqu'on opérera sur une surface courbe). 
Du reste, dans les réseaux un peu compliqués, il arrive 
souvent que les équations de la seconde et de la troisième 
catégorie rentrent en partie les unes dans les autres , et que, 
par suite, leur nombre doit être réduit. Au contraire, il 
pourra arriver, mais seulement dans des cas assez rares , 
que l'on adjoigne quelques équations nouvelles à celles de 
la seconde catégorie; c'est ce qui aura lieu lorsque le ré- 
seau contiendra des polygones non divisées en triangles; on 

7- 
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pourra alors introduire des équations relatives aux ligures 
qui ont plus de trois côtés. Dans une autre occasion , nous 
reviendrons avec plus de détails sur ces diverses circon- 
stances, dont l'examen nous éloignerait en ce moment de 
notre but. Nous ne pouvons cependant nous dispenser de 
faire ici une remarque indispensable à ceux qui voudraient 
faire l'application rigoureuse de notre tbéorie : nous suppo- 
sons toujours que les quantités désignées par y, v 1 v", etc. , 
ont été observées immédiatement, ou déduites d'observa- 
tions telles, que leurs déterminations soient indépendantes 
les unes des autres, ou, au moins, puissent être regardées 
comme telles. Dans la pratique la plus ordinaire, on observe 
les angles que l'on peut regarder comme étant les éléments 
^, i>', y", etc., eux-mêmes. Mais on ne doit pas oublier que 
si le système contient, en outre, des triangles dont les angles 
n'aient pas été directement observés et aient été déduits de 
ceux que Ton connaissait, par des additions ou soustrac- 
tions, ces angles ne devront pas être mis au nombre des 
grandeurs déterminées par l'observation, et l'on devra les 
faire entrer dans le calcul comme des fonctions des éléments 
qui ont servi à les former. Il en sera autrement si l'on 
adopte la métbode d'observations de M. Struve [j4strono~ 
mische Nachrichten, II, page 43 1), qui consiste à détermi- 
ner toutes les directions autour d'un même sommet, en les 
rappor^nt toutes à une seule et même direction arbitraire. 
Les angles mesurés ainsi seront pris alors pour v , v\ t>", etc. , 
et les angles des triangles se présenteront tous comme des 
différences. Les équations de la première catégorie devront, 
dans ce cas, être supprimées comme superflues, car elles 
seront identiquement satisfaites. Le procédé que j'ai suivi 
moi-même dans les triangulations exécutées pendant ces 
dernières années , diffère des deux méthodes précédentes -, 
on peut cependant l'assimiler, quant au résultat, avec le 
procédé de M. Struve, en ce sens que, dans chaque station , 
on doit regarder y, v\ v" , etc., comme les angles formés 
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par les directions qui en partent, avec une même ligne 
arbitrairement choisie. 

Nous donnerons deux exemples : le premier se rapporte 
au premier mode d'opération, et le second est relatif a 
des observations faites d'après la seconde méthode. 

23. 

Le premier exemple nous sera fourni par l'ouvrage de 
M. Krayenhof : Précis historique des opérations trigono- 
métriques faites en Hollande. Nous chercherons à cow- 
penserla partie des observations relatives au terrain com- 
pris entre Harlingen, Sneck , Oldcholtpade , Ballum , 
Lecuwardcn,.Doekuin , Drachten , Oosterwolde et Grônin- 
gen. Entre ces points , on a formé neuf triangles numérotés, 
dans l'ouvrage cité, 121, 122, 123 ,124, 126, 127, 1 r >8 , 
1 3 1 , i32, Les angles observés sont les suivants : 

Triangle 121. 

0. Harlingen 5o"58' i5",238 

1. Leeuwarden.... . 82.47 v i5,35i 

2. Ballum 46 >4- 27,202 

Triangle 122. 

3. Harlingen 5i . 5.39,717 

4. Sneek 70.48. 3*3,445 

5. Leeuwârden.. 58. 5.48,707 . 

Trianclk 123. 

6. Sneek 49 3o.4o,o5i 

7. Drachten 4 2 «52. 59, 38a 

8 Leeuwarden 87 . 36 . 2 1 ,057 

TwiAXGLE 124. 

t) Sneek 4^'36. 7»49 2 

10. Oldcholtpade 67. 5s. 0,048 

11. Drachten 66 3i.5G,5i3 



Digitized by GcjPfgle 



( M ) 

Triangle 125. 

t2. Drachten 53° 55' z^'ilfô 

i3. Oldeholtpade 47-4 8 -5^,58o 

\\. Oosterwolde 78.15.42,347 

Triangle 127. 

i5. Leeuwarden 59.24. o,645 

»6. Dockum 76.34. 9,021 

17. Ballum 44- i.5i,o4o 

Triangle 128. 

18. Leeuwarden 72. 6 32,o43 

19. Drachten 4^ -53. 27,1 63 

ao. Dockum 6i. o. 4;494 

Triangle 1 3 1 . 

21. Dockum 57 1.55,292 

22. Drachten 83.33. i4,5i5 

23. Groningen 39.24.52,397 

Triangle i32. 

y>4 • Oosterwolde 8 1 . 54 • 1 7 , r \ 47 

25. Groningen.. 31.52.46,094 

26. Drachten 66 . 1 2 . 57 , 246 

La considération de ces triangles montre que les vingt-sept 
angles, directement fournis par l'observation, ont entre eux 
treize relations nécessaires, savoir : deux de la première 
espèce, neuf de la seconde, et deux de la troisième. Mais 
il n'est pas utile d'écrire ici toutes ces équations sous 
leur forme finie, car pour le calcul nous avons besoin seu- 
lement des quantités, désignées dans la théorie générale, 
par X, a, a', a f \. . . , , b, />', //', . . . , c'est pourquoi 
nous écrirons immédiatement les équations (i3), qui met- 
tent ces quantités en évidence. Au lieu de g, e', e", etc., 
nous écrirons simplement ici (o), (1), (2), etc. De cette ma- 
nière, aux deux équations du premier genre répondent les 



Digitized by Google 



( «-3 ) 

suivantes : 

(O ■+■ (5) -h (8) H- (i5) + (18) = - 2", 197, 
(7 ) + (» 0 + l 12 ) + («9) + ( 22 ) + ( 26 ) = ~ °">436. 
Nous trouvons ensuite, pour les excès sphéroïdiques des 
neuf triangles : i", 7 49ô »%«47» l "M*\ i",6 9 8j o",8 7 3j 
i",i675 i",i°4 -, a", 161 *, i",4«3- Nous aurons alors l'équa- 
tion de condition du second genre : 

„<•) 4- «,<•) + pO _ i8o°.o' i",749 = o, 

et ainsi des autres, et nous avons les neuf équations sui- 
vantes : 



(0) + 


(0 + ( 2 ) 




-3",958, 


(3) + 


(4) + (5) 




-1-0,722, 


(6} + 


(7)+ (8) 




— 0,753, 


(9) + 


(10) H- (11) 




-4- 2,355, 


(ia) + 


(i3)+ (i4) 




— 1 ,201 , 


(i5) ^ 


(«6) + (.7) 




— 0 ,«46 1 , 


(«B) + 


(19)-+- (20) 




+ 2,596, 


(ai) + 


(22) H- (23) 




-h 0,043, 


( 2 4)+- 


(25) -4- (26) 




— 0,616. 



Les équations de condition du troisième genre s'expriment 
plus facilement par le moyen des logarithmes : la première 
est 

log sin - o",583) - log sin (*» — o",583) - log sin (p 5 - o",382) 

4- log sin — o",382) - log sin (p 6 -o",4 1 4) 
-4- log sin (,' — o",4i4) — logsin(p"~ o",38 9 ) 
* -Mogsin(('"-o",389) — logsin («'"- o",368) 
-h logsin(p î0 -o",368)=:o. 

11 semble inutile de développer l'autre sous forme finie. A 
ces deux équations répondent les suivantes, dans lesquelles 
les coefficients se rapportent à la septième décimale des lo- 
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garithmes vulgaires (*) : 

17,068(0) — 20,174(2)— 16,993 (3) -h 7,328(4)— 17,976(6) 
+ 22,672(7)— 5,028(16)4-21 ,780(17)— 19,710(191 
H- 11 ,671 (20) = — 371, 

17,976(6)— 0,880(8) — 20,617(9)4-8,564(10) — 19,082(13) 
4- 4,3 7 5 (.4) 4- 6,798(18)- 1 1 ,671 (20) h- i3,65 7 (2i j 
— 25,620(23) — 2,995 (24) 4- 33,854 (25) = -h 370. 

Aucune raison ne nous portant à attribuer des poids iné- 
gaux aux diverses observations, nous supposerons 

p* = /?' =r p> = . . . = I . 

En désignant les corrélatifs des équations de condition dans 
Tordre même où ces équations ont été écrites, par 

A,B,C,D,E,F,G,H,I,K,L,M,N, 
nous les déterminerons par les équations suivantes : 

— 2", 197 = 5A + C + D-f E + H-HlH-5,9i7N, 

— o,436 = 6 B + E + F + G4-I + K + L + 2,962 M > 

— 3,958 = A 4- 3C — 3, 106 M, 
4- 0,722 = A -h 3 D — 9,665 M, 

-o, 7 53 =A + B + 3E-+- 4,696»! 4- 17,096 N. 
-+- 2,355 — B 4- 3F — i2,o53N, 

— 1,201 = B4- 3G — i4»7°7 N > 
^ 0,461 = A 4- 3 H -h 16,752 M, 

— 2 , 596 = A 4- B 4- 3.1 — 8,o3gM -4,8741V, 
4- o,o43 — B-+-3K — 11 ,963 N, 

— 0,616 = B + 3L + 3o,859N, 

— 37 1 ~ 2,962 B — 3, 106 C — 9,665 D 4,696 E 

4-16, 75?. H — 8 , 039 I 4- 2902 , 27 M 

— 459,33 N, 

h 370 = 5,917 A 4- 17,096F — i2.o53 F — 14,707 G 

— 4,8741 — ii,963K*4- 3o,859L 

— 459,33 M 4- 3385, 96N 

(*) Ces coefficients ont été tous multipliés, après la diflerentiation , par 
10', et divisés par 206264,8= l8o, ^° .P, pour convertir les erreurs en 
secondes. J. B. 
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Nous en déduisons par l'élimination : 



A = — o,5 9 8 
B = — o , 255 

C = — 1,234 

D = 4- 0,086 
E = - 0,447 
F = 4- f,35i 



H = 
I = 
K = 
L = 
M = 
N = 



4- o,659 
4- 1 ,o5o 
H- 0,577 

— 1 ,35i 

— 0,10979a 
4- o, 1 19681, 



G = 4 0,271 

Les erreurs les plus plausibles sont enfin données par les 
formules 

(0) = C+ 17,068 M 

(1) = A + C 

(2) = C — 20, 174 M 

(3) = D— 16,993 M, 

et nous obtenons les valeurs numériques suivantes , aux- 
quelles nous adjoignons, pour qu'on fasse la comparaison, 
les corrections adoptées par M. Krayenhof : 





de Kn. 






DE KR. 


(0) =- 3",io8 


— 2",090 


(i4) 


= H-o',79 5 


4- 2",4°o 


(l) =r — 1,832 


4- O, I 16 


(i5) 


= 4- 0,061 


4- 1,273 


(2) =4- 0,981 


— i>9 82 


(16) 


— 4- T ,21 I 


4- 5,q45 


(3) =4- 1,952 


4- 1 »722 




= — 1 ,732 


~ 7,674 


(4) = — 0,719 


4- 2,848 


(18) 


= 4- 1 , 265 


4- I ,876 


(5) = — o ? 5i2 


— 3,848 


{'9) 


— 4- 2,959 


4- 6,?.5i 


(6) =4- 3,648 


— 0,137 


(20) 


= — 1 ,628 


— 5,53o 


(7) =— 3,221 


4- 1 ,000 


(ai ) 


= 4- 2,211 


4- 3,486 


(8) =— 1,180 


- 1,614 


(22) 


— 4- 0,322 


-3,454 


(9) =~ ! > 116 


0 


(23) 


= - 2,489 


0 


(10) = 4- 2,376 


4- 5,928 


(*4) 


= - 1,709 


4- 0,400 


(n) =4- 1,096 


— 3,570 


(25) 


= 4- 2,701 


4- 2,o54 


(12) = 4- 0,016 


4- 2,4*4 


(26) 


— — 1 ,606 


- 3,077 


(i3) = — 2,oi3 


— 6,014 


> 

! 







La somme des carrés de nos corrections est97,8845 5 l'erreur 
moyenne, telle que l'indiquent les 27 angles observés, est 
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par conséquent , 

La somme des carrés des corrections de M. Krayenhof , 
est de 34i)4 201 • 

24. 

Les triangles dont les sommets, dans la triangulation du 
Hanovre, ont été placés à Falkenberg, Breithorn, Hausel- 
berg, Wulfsode et Wilsede, nous fourniront un second 
exemple. 

On a observé les directions suivantes : 

A la station de Falkenberg. 

0. Wilsede i%*) 0 fa' 3o" 3i i 

1. Wulfsode 225.9.39.676 

2. Hauselberg 266. 1 3. 56. 239 

3. Breithorn 274 • i4- 4^.634 

A la station de Breithorn. 

4- Falkenberg 94. 33. 4°. 7^5 

5. Hauselberg i22.5i.23 o54 

6. Wilsede i5o. 18. 35. 100 

A la station de Hauselberg. 

7. Falkenberg 86.29. 6.872 

8. Wilsede 154.37. 9 624 

9. Wulfsode 189. 2.56.376 

10. Breithorn 302.47.37.732 

A la station de Wulfsode. 

11. Hauselberg 9. 5.36.593 

12. Falkenberg 4^- 2 7 33.556 

13. Wilsede 1 1 8 . 44 . 1 3 . 1 59 

A la station de Wilsede. 

14. Falkenberg 7 . 5 1 . 1.027 

15. Wulfsode 298.29.49.519 

16. Breithorn 33o. 3. 7.392 

17. Hauselberg 334*25. 26.746 



Digitized by Google 



( ll >7 ) 

Ces observations permettent de former sept triangles. 

Triangle I". 

Falkenberg 8° o' 47 "*3g5 

Breithorn . 28 . 1 7 . l\i , 299 

Hauselberg. ...... 1 4^ • 4 1 • 2 9» ! 4° 

Triangle H. 

Falkenberg 86 . 27 . 1 3 , 323 

Breithorn 55 . 44 • 54 » 345 

Wilsede 37.47.53,635 

Triangle III. 

Falkenberg 4 1 4- I 6>563 

H auselberg 102. 33 . 49 , 5o4 

Wulfsode 36 . 2 1 . 56 ,963 

Triangle IV. 

Falkenberg 78.26 25,928 

Hauselberg 68. 8. 2,75?. 

Wilsede. 35 25.34,281 

Triangle V. 

Falkenberg 37.22.9, 365 

Wulfsode ........ 73. i6.39,6o3 

Wilsede.. 69.21.11,508 

Triangle VI. 

Breithorn 27 27.12,046 

Hauselberg 148.10,28,108 

Wilsede 4- 22 * I 9>354 

Trianci.e VII. 

Hauselberg 34.25.46,752 

Wulfsode 1 09 38 . 36 , 566 

Wilsede * 35.55.37,227 

Nous avons ici sept équations de condition du second 
genre (il n'y a pas lieu évidemment d'eu former du premier 
•genre) \ pour les former, nous devons chercher, avant tout, 
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les excès sphéroïdiques des sept triangles, et pour cela il est 
indispensable de connaître la longueur d'un côté. Celui 
qui réunit Wilsede à Wulfsode est 22877™, 94. On en con- 
clut, pour les excès sphéroïdiques des divers triangles : 
I...o",2oa; II. . .-/',442; III. . . i ff ,a575 IV . . . î^ojcj ; 
V...i",i)57; VI...o",32i\ VII... i",2 9 5. 

Si Ton désigne par i> (0) , y (î) , i' (3) , etc., les angles qui 
déterminent les directions indiquées plus haut,' et mar- 
quées des mêmes indices, les angles du premier triangle 
seront 

,,(3) — p<o , „(0 _ , 360 + 1»<»> — c( ,9 >, 
et la première équation de condition est par conséquent 
~ m(î) 4- „(3) __ p (o + + „('o) 4. 179059' 59", 798= o. 

Les six triangles restants fourniront six équations analo- 
gues, mais un peu d'attention montrera que ces équations 
ne sont pas indépendantes 5 la seconde est en effet identique 
avec la somme de la première, de la quatrième et de la 
sixième; la somme de la troisième et de la cinquième est 
identique avec celle de la quatrième et de la septième : c'est 
pourquoi nous négligerons la seconde et la cinquième. Au 
lieu des équalions restantes sous forme finie, nous écrirons 
ici les équations correspondantes (i3), en subsituant aux 
notations £, e', e", etc.. (o), (1), (2), etc. : 

- i",3G8=-.(2)+(3) -(4)+ (5) H- (7) -(10), 
-4- 1 ,773 = -(1 ) + (a) ~f 7 )-+-(9) - (11) + (1a), 
+ 1,042 = -(o) + (2)-( 7 ) + (8) + (i4) — (17), 

- o,8i3 = - ( 5)+ (6) - (8) + (10;- (16) + (17), 
-o,75<> = — (8) + (y) -(n)+ (i3) — (i5) + (i 7 ). 

On peut obtenir, au moyen des triangles du système, 
huit équations du troisième genre, et pour cela il est per- 
mis de combiner trois des quatre triangles I , Il , 1 V, VI , ou 
des triangles III, IV, V, VII; cependant un peu d'attention m 
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montre qu'il suffit d'en considérer deux appartenant res- 
pectivement aux deux systèmes de triangles et que celles-là 
comprendront toutes les autres. 

Nous aurons ainsi , pour sixième et septième équation de 
condition, 

logsin( p(» — — 0V67) — Iogsin (p( j > — i>( 4 > — o",o67) 
-h Iogsin (o(" ' — ,»(»>_ o",64o) — logsin(p( î > — o",64o) 
4-logsin(f( 6 > — yf*)— o",io7)— logsin(Mf ,î )— p( ,6 '—o",i07)=ro, 

logsin(^>~ c(" _o",4i9;— Iogsin (<,(»>—<,(">— o",4 19) 
+ Iogsin (p( v (" ) — o w ,64o) — log sin (<>(» — »»(•> — o",64o) 
+ Iogsin (f< ,3 >— <>(">— o",439.) — Iogsin (p(">— o ,/ ,43?.)=o, 

auxquelles répondent les équations 

+ ,5 = -h4,3i (o) - i53,88(2) + 149,57 (3) 

4-39,11(4)- ;9,64(5)+4o,53(6) 

-h 3.1,90 (i4) 4- 275,39(16) — 307,29(17^, 

- 3 = + 4,3i (o) -24,16(1) + 19,85(2) 
+ 36,n (11) — 28,59 (12) — 7,52(1 3) 
+ 3 1, 90(14) + 29,06(15) — 60,96(17). 

Si nous attribuons la même certitude aux diverses direc- 
tions, en supposant p {0) = p (i} = p (i) ... = 1 ,les corrélatifs 
des sept équations de condition étant désignés par A, B, 
C, D, E, F, G, leur détermination dépendra des équations 
suivantes : 

— i,368 = + 6 A — 2B— 2C — 2D + 184,72 F — 19,85 G, 
-+-1 ,773 =— 2A+6B + 2C+ 2E — 1 53,88 F —20,69 c, 

4- I ,o42 =— 2A + 2B+6C — 2D — 2E -H l8l,00 F 

4- 108, 4o G , 

— o,8i3 = — 2A — 2C + 6D4- 2E — 46?.,5i F — 60,960, 

— 0,750 = 4- 2 B — 2C + aD + 6E — 307,29F — 1 33,65 G 

+ 25 = 184,72 A — i53,88B + i8i,ooC — 462,51 D 

— 307,29 E 4- 224868F + 16694,1 G, 

— 3 = — 19,85 A— 20,698 + 108,40 C — 60,960 

- 133,65e + 16694,1 F + 8 7 52,3 9 G. 



Digitized by google I 



( 110 ) 

Nous en déduisons, par élimination, 

A = — 0,225, 

B = -f- o,344, 

C = — 0,088, 

D = — 0,171, 

E = — o,323, 

F = -h 0,000215915, 

G — — 0,005474620 , 

et les erreurs les plus probables sont données par les for- 
mules 

(0) = — C-h4,3iF H-4,3iG, 

( 1 ) — B — 24,i6G, 

( 2 )=:-A + B + C- i53,88F -+- i 9 ,85G, 



d'où Ton déduit les valeurs numériques suivantes : 



(0: 




+ o",o65 


(9) 


= -f- o",o2 1 


(>; 




0,212 


(10) 


= -h 0,054 






-H o,33g 


(") 


=r— 0,219 


(3) 




— 0, 193 


(12) 


= H- 0 , 5o I 


f4) 




-H o,?.33 


(-3) 


3= — 0 , 282 


(5; 




— 0,071 


(>4) 


= — o,256 


(6) 




— 0,162 


(.5) 


— -h 0,l64 


(7) 




— o,48i 


(.6) 


= -f- o,23o 


(8) 




-f- 0,406 


(■7) 


= - o>»3 9 . 



La somme des carrés de ces erreurs est égale à 1,2288; l'er- 
reur moyenne résultant des 18 directions observées est, par 
conséquent, 

/ 1 ,2288 , 
y/ — — = 0^,4190. 

25. 

Afin de donner un exemple de la dernière partie de notre 
théorie, cherchons dans quelle précision les observations 
compensées déterminent le côté Falkenberg-Breithorn , 
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au moyen du côté Wilsede- Wulfsode. La fonction m, par 
laquelle il est exprimé , est, dans ce cas, 

_ fi m sin („(■») _ „<■») — p",652). sin (p("> — p("> — o»,8i4) 

sa valeur déduite des observations corrigées est 

26766 m ,68. 

La ditférentiation de cette équation fournit, en exprimant 
<iV (0) , etc., en secondes, 

du = [ o™, 1 6991 ( rfp(«) ) — ^(-) ] 4. o m ,o8836 (rf^î — 

— [o-,o3899(rA'i»>) — f/p(*>] H- o m , 16731^»)— «/„<••>); 

on déduit de là : 

(«/)=— o,o8836, 
(*/)=-h o,i3o 9 2, 
(c/) = — 0,00260, 
+ 0,07895, 
(<•/) =r -h o,o38 9 C), 
(//)=-4o,i3i5, 
[gl) = + 10,9957, 
( //) = 4- o,i3238, 

Les méthodes indiquées plus haut donnent, en prenant le 
mètre pour unité de longueur, 

^=0,08329 ou P= 12,006. 

On en conclut que l'erreur moyenne à craindre dans la va- 
leur du côté Falkenberg-Breithorn est o m ,iSS6 m (m dé- 
signant Terreur moyenne à craindre dans les directions 
observées, cette erreur étant exprimée en secondes), et, 
par conséquent, si nous adoptons la valeur de m annoncée 
plus haut, cette erreur moyenne à craindre est o m ,i209. 

Au reste, l'inspection du système de triangles montre 
immédiatement qu'on pouvait complètement laisser de côté 
la station Hauselberg, sans rompre le réseau qui réunit les 
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quatre autres. Mais il ne serait pas permis pour cela de 
supprimer les opérations qui se rapportent à ce point, car 
elles contribuent certainement à augmenter la précision de 
l'ensemble. Pour montrer plus clairement quel accroisse- 
ment de précision en résulte, nous terminerons en faisant 
de nouveau le calcul, après avoir exclu tous les résultats qui 
se rapportent au point Hauselberg. Des dix-huit directions 
mentionnées plus haut, huit cessent alors de servir, et les 
erreurs les plus plausibles sur celles qui restent , sont 



(o) = 


-*-o",327 


(12) = 


-f- o",2o6 


(') = 


— o , 206 


(i3) = 


— o , 206 


(3) = 


— O, 121 


(•4) = 


+ 0, 32; 


'4) = 


-h O, 121 


(i5) = 


-h 0, 206 


(6) = 


— 0, 12 1 


(16) = 


H- 0, 121 



La valeur du côté Falkeuberg-Breithorn devient alors 
26766™ ,63, résultat peu diiïerent de celui qui a élé obtenu 
plus haut. Mais le calcul du poids donne 

p = 0,1 3o82 ou P =. 7,644» 

et Terreur moyenne à craindre est, en mètres, 

0,36169m = o m ,i5i5. 

On voit que par l'adjonction des opérations qui se rap- 
portent à Hauselberg, le poids de la détermination du 
côté Falkenberg-Breithorn est augmenté dans le rapport 
de 7,644 à 12,006, c'est-à-dire dans le rapport de l'unité 
à 1,571. 
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NOTES. 

NOTE L 

EXPOSITION DE LA MÉTHODE DES MOINDRES CARRÉS. 
(Extrait du Theoria Motus Corporum cœleslium.) 

Abordons maintenant une recherche beaucoup 

plus générale et des plus fécondes dans toute application 
du calcul aux phénomènes naturels. Soient V, V, V, etc., 
fonctions des v inconnues g, r, $, etc. , et supposons 
que des observations directes aient donné, pour ces fonc- 
tions, les valeurs 

V=M, V'=M', V // =M ,/ ,.... 

En général, le calcul de ces inconnues constituera un pro- 
blème indéterminé , déterminé ou plus que déterminé, sui- 
vant que Ton aura 

f*<v, f»=v, ou p>v (*). 

Nous ne nous occuperons ici que du dernier cas, dans lequel 
évidemment il ne serait possible d'obtenir une représenta- 
tion exacte de toutes les observations , que si ces observa- 



(*) Si, dans ce troisième cas, /t + i- v des fonctions V, V, V, etc., 
pouvaient être regardées comme des fonctions de toutes les autres, le pro- 
blème deviendrait plus que déterminé relativement à ces fonctions , mais 
indéterminé relativement à p, q, r, s, etc. On ne pourrait pas en déduire 
les valeurs de ces dernières , même si les valeurs des fonctions V, V, V, etc., 
étaient d'une exactitude absolue : mais nous excluons ce cas particulier tic 
nos recherches. 

o 8 
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tions n'étaient affectées d'aucune erreur. Mais comme cela 
n'a jamais lieu dans la nature, on devra regarder comme 
possible tout système de valeurs des inconnues p, q, r, 
s, etc., desquelles résultent, pour les fonctions V — M, 
V — M', V" — M", des valeurs qui ne surpassent pas les 
limites des erreurs que Ton peut commettre dans les obser- 
vations, mais on ne doit pas regarder tous ces systèmes 
possibles comme jouissant du même degré de probabilité. 

Supposons d'abord, dans toutes les observations , un état 
de choses tel, qu'il n'y ait pas lieu de regarder l'une d'elles 
comme plus exacte qu'une autre, c'est-à-dire, que l'on 
doive regarder des erreurs égales dans chacune d'elles 
comme également probables. La probabilité qu'une erreur A 
soit commise dans l'une des observations sera une fonction 
de A, que nous nommerons <p (A). Quoique cette fonction ne 
puisse être assignée d'une manière précise, on peut du 
moins affirmer qu'elle doit devenir maximum pour A =r o, 
avoir dans la plupart des cas la même valeur pour des va- 
leurs de A égales et de signes contraires , et, enfin, s'éva- 
nouir quand on donne à A une valeur égale ou supérieure 
à l'erreur maximum ; <p ( A) doit donc , à proprement parler, 
être rapportée à la classe des fonctions discontinues, et, si 
nous nous permettons, pour la facilité du calcul , d'y sub- 
stituer une fonction analytique, il faudra que cette dernière 
soit choisie de telle sorte qu'elle tende rapidement vers o à 
partir de deux valeurs de A, Tune supérieure, l'autre infé- 
rieure à o , et qu'en dehors de ces deux limites on puisse la 
considérer comme nulle. Or la probabilité que l'erreur soit 
comprise entre A et une quantité A d& qui en diffère 
infiniment peu, sera exprimée par <p (A) .r/A, et, par suite, 
la probabilité que l'erreur est comprise entre D et D', par 




Cette intégrale, prise depuis la plus grande valeur néga 
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tire de A jusqu'à sa plus grande valeur positive, ou plus 
généralement depuis A = — oo jusqu'à A = oo , devra né- 
cessairement être égale à i . On aura donc 

y (A) dA =r i. 

■oo 

Supposons donc qu'on ait un système déterminé de valeurs 
des quantités /?, q, r, s, etc. : la probabilité que l'obser- 
vation donnera pour V la valeur M, sera exprimée par 
<p (M — V) , après qu'on aura substitué dans V les valeurs 
de q, r,5, etc.-, de môme 9 (M' — V), <p (M" — V"), etc., 
exprimeront les probabilités pour que les observations don- 
nent aux fonctions V, V", etc., les valeurs M', M*, etc. 
C'est pourquoi , tant qu'on pourra considérer toutes les ob- 
servations comme des événements indépendants les uns 
des autres , le produit 

?(M— V) ? (M'-V') ? (M"-V")... -=12, 

exprimera la probabilité que toutes ces valeurs résulteront 
en même temps des observations. ' 

2. 

De même qu'en se donnant des valeurs quelconques des 
inconnues, il en résulte, avant toute observation , une pro- 
babilité déterminée pour un système de valeurs des fonc- 
tion V, V ', V", etc., de même , après que l'observation aura 
donné pour ces fonctions des valeurs déterminées, il en 
résultera pour chaque système de valeurs des inconnues 
qui en découleront , une probabilité déterminée : car il est 
clair qu'on devra considérer comme les plus probables les 
systèmes qui donnent à l'événement observé la plus grande 
probabilité. L'appréciation de cette probabilité peut s'ob- 
tenir par le théorème suivant : 

Si> en adoptant une certaine hypothèse H, la probabilité 
d'un événement déterminé Eest h, mais qu'en adoptant 
une autre hypothèse H', exclusive de la première et ayant 

8. 
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à priori la même probabilité, la probabilité du même évé- 
nement soit h' : je dis que lorsque T événement E aura eu 
lieu y la probabilité que H soit la vraie hypothèse sera à la 
probabilité que H ; soit la vraie hypothèse comme h est à h'. 
Pour le démontrer et afin de distinguer toutes les circon- 
stances d'où peut dépendre, soit que l'hypothèse H ou H', 
ou toute autre ait lieu, l'arrivée d'un événement E ou d'un 
autre événement, formons un système des cas différents qui 
peuvent se présenter et que nous regarderons comme éga- 
lement probables à priori (c'est-à-dire tant qu'il y a doute 
si c'est l'événement E ou un autre qui aura lieu). Ces cas 
peuvent être ainsi distribués : 



NOMBRE DES CAS. 


HYPOTHÈSE 

propre à ces cas. 


ÉVÉNEMENT 

qui doit en résulter. 


m 


H 


El 


n 


H 


Différent de E 


m' 


H' 


E 


n' 


H' 


Différent de E 


m" 


Différente de H et de H' 


E 


n" 


DifTérentedcHet de H' 


Différent de E 



On aura d'après cela 



m -f- n m -+- n' 

Or, avant l'arrivée de l'événement, la probabilité de l'hypo- 
thèse H était 

m -+- n 

m -f- n -4- m' -f- n' +- m" -+- n"' 

Après l'événement qui exclut n -f- n' -f- n" cas, parmi ceux 
qui sont possibles, cette probabilité sera 

m 

m -f- /// -f- m" 
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De même les probabilités de l'hypothèse H' avant et après 
l'événement sont respectivement 

m' -h »' m' 



et 



m -h n H- /m' -h «' H- jw" H- «" w 4- i»' H- m" ' 

mais, comme on a supposé que les hypothèses H et H' avaient 
avant l'événement la même probabilité, on aura 

m -f- n = m' H- n' y 

d'où résulte immédiatement la vérité du théorème. 

Si maintenant on suppose qu'on n'a , pour déterminer les 
inconnues, que les observations 

V=M, V' = M', V" = M", 

et que tous les systèmes de valeurs des inconnues étaient éga- 
lement probables avant ces observations, il est visible que 
la probabilité d'un certain système, après ces observations, 
sera proportionnelle à ft. C'est-à-dire que IQdpdqdr. . . 
exprimera la probabilité que les valeurs des inconnues 
soient respectivement comprises dans les limites infini- 
ment voisines p et p -h dp, q et q -f- dq, r et r -h dr, . . . , 
X représentant une quantité indépendante de p, q, r, 
y, etc.; et l'on aura évidemment 

-=z I I j ' • • H dp dq dr. , . • 

«/_ ce J— oc J — ce 

3. 

De là résulte naturellement que le système le plus pro- 
bable des valeurs dep, q, r 9 etc., correspondra au maxi- 
mum de lî, et se tirera des v équa lions 

dLl di\ dS2 

— = o, — = o, -~=o, • • y 

dp dq dr 

si Ton pose 

V-M = ,>, V'-M'=/,V''_M'W',... et 
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ces équations prendront la forme suivante : 



De là résulte qu'on pourra obtenir par Félimination une 
solution pleinement déterminée du problème, dès que. la 
nature de la fonction y' sera connue. Mais comme cette 
fonctioa ne peut être définie à priori , abordons la question 
à un autre point de vue et cherchons une fonction acceptée 
tacitement comme base, en vertu d'un principe simple et 
généralement admis. Or on a coutume de regarder comme 
un axiome l'hypothèse que si une quantité a été obtenue par 
plusieurs observations immédiates, faites avec le même soin 
dans des circonstances semblables, la moyenne arithmétique 
des valeurs observées sera la valeur la plus probable de cette 
quantité, sinon en toute rigueur, du moins avec une grande 
approximation, de telle sorte que le plus sûr soit toujours 
de s'y arrêter. Si donc l'on pose 

V = V'=V".,.=:^ 

et 

M + M'+M" + ... 

,= 

on devra avoir en général 

?' ( M -p) + ?' ( M' - p ) 4- M" -/>) + ... = o 

pour toute valeur entière et positive de u. Faisant en- 
suite 

M'=M"...= M- P N, 

» 

on aura généralement 

f [((«'- i)B] = (i -j*)f'(-N), 
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d'où Ton tire facilement que ff -~ doit être en général une 



J (A) 

constante h. On aura donc 

log? (A) = y H-const. = |/ A 1 -Mogx, 



d'où 



Or on voit facilement que la constante h doit être négative, 
pour que 12 puisse devenir maximum : posons donc 



l k = _ /<> 

2 



et comme, d'après un élégant théorème de Laplace, on a 

/°° -*»A» , 4 sfi 
ao n 

notre fonction deviendra 

4. 

La fonction que nous venons de trouver ne peut pas ex- 
primer, en toute rigueur, la probabilité des erreurs, puis- 
que les erreurs possibles étant toujours renfermées entre 
certaines limites, la probabilité d'erreurs plus grandes de- 
vrait être toujours nulle, tandis que notre fonction a toujours 
une valeur finie. Cependant ce défaut, que présenterait 
également toute autre fonction analytique, n'a aucune*im- 
portance dans les applications, parce que la valeur de notre 
fonction décroît si rapidement, pour peu que A A ait une 
valeur considérable , qu'on peut, en toute sûreté, la re- 
garder alors comme équivalente à o. D'ailleurs, la nature de 
la question ne permettra jamais d'assigner les limites des 
erreurs avec une rigueur absolue. 

Au reste, la constante h peut être regardée comme ser- 
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vant de mesure à la précision des observations. Si en effet la 
probabilité de Terreur A dans un système d'observations 
est exprimée par 

h — a 1 

et dans un autre système d'observations plus ou moins 
exactes que les premières par 

h' _A'*A* 

T e 

la probabilité que dans une observation du premier sys- 
tème l'erreur soit comprise entre les limites — à et -H à, 
sera exprimée par 



x 



' A — A* A* . 
— e f/A, 



et de même la probabilité que l'erreur d'une observation du 
second système soit comprise entre les limites — $' et -+-(?', 
sera exprimée par 

,<r h' — h'* A* 



X 



- c dà ' 

-S' \/7t 

or ces intégrales sont manifestement égales lorsqu'on a 
Si , par exemple, on a 

une^errcur double dans le premier système sera commise 
aussi facilement qu'une erreur simple dans le second , de 
sorte que les dernières observations, pour nous servir 
d'une expression consacrée par l'usage , jouissent d'un degré 
de précision deux fois plus grand. 

5. 

■ 

Voici maintenant quelques conséquences de celte loi. Il 
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est clair qu'il faut, pour que le produit 

devienne maximum , que la somme 

■+■ p * H- p'" -H . . . 

devienne minimum. Donc le système de valeurs des incon- 
nues/?, <jr, r, s, etc., le plus probable correspond au cas où 
les carrés des différences entre les valeurs observées et les 
valeurs calculées des fonctions V, V, V, etc., donnent la 
somme la plus petite possible, pourvu que toutes les obser- 
vations soient également présumées précises. 

Ce principe, qui est de la plus grande utilité dans toutes 
les applications des mathématiques à la philosophie natu- 
relle, doit être regardé comme un axiome, au même titre 
que le principe qui nous fait adopter la moyenne arithmé- 
tique des valeurs observées d'une même quantité comme la 
valeur la plus probable de cette quantité. 

Le principe s'étend sans peine au cas d'observations 
d'une précision inégale. Car si les précisions des observa- 
tions par lesquelles ou a trouvé 

V = M, V' = M', ¥" = 11",..., 

sont représentées respectivement par h, h\ h", etc., c'est- 
à-dire si Ton suppose que des erreurs réciproquement pro- 
portionnelles à ces quantités puissent être commises avec la 
même facilité, il est clair que cela revient au même que si, 
pardesobservations d'une égale précision (représentée pan), 
les valeurs des fonctions /iV, /t'V, A"V", etc., avaient été 
trouvées égales à AMj.A'M', A "M", etc.} c'est pourquoi le 
système le plus probable de valeurs des quantités p, q , r, 
5, etc., sera celui où la somme 

ccsl-à-dirc oit la somme des canes des différences entre 
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les valeurs observées et calculées y multipliés respectivement 
. par les carrés des nombres qui expriment le degré de pré- 
cision^ devient un minimum. Par là, il n'est pas même né- 
cessaire que les fonctions V, V, V", etc., se rapportent à 
des quantités homogènes, niais elles pourront représenter 
des quantités hétérogènes (par exemple des secondes d'arcs 
et de temps); pourvu que Ton puisse estimer le rapport 
des erreurs qui, dans chacune de ces grandeurs , peuvent 
être commises avec la même facilité. 

6. 

Le principe exposé dans l'article précédent se recom- 
mande aussi par cela qu'il réduit le calcul numérique des 
inconnues à un algorithme très-expéditif , quand les fonc- 
tions V, V, V", etc., sont linéaires. Supposons 

V — M = v = — m -\- a p ~\- b q-hcr-\-ds-\-..., 

V — M' = v' = — m' -h b' e r + tts + . . . , 
V" — M" ^v"—~ m" -h a"p -h b" q -+- c"r+A+.. . , 



Posons 

c ç + c ' v' -+- c" v" -f- . . . = R, 
rfi>H-rfV -h iî'v N -h. . .= S, 



alors les v équations de l'art. 3, qui déterminent les valeurs 
des inconnues , seront 

P=o, Q = o, R = o, S = o,..., 

si nous supposons les observations également bonnes, cas 
auquel nous pouvons ramener tous les autres , comme nous 
l'avons montré dans l'article précédent. On a ainsi autant 
d'équations linéaires que d'inconnues : on les résoudra par 
la méthode ordinaire. 



Digitized by Google 



( '^3 ) 

Voyons mainteuant si cette élimination est toujours pos- 
sible ou si elle peut donner une valeur indéterminée ou im- 
possible. Il résulte de la théorie de l'élimination que le se- 
cond ou le troisième cas aura lieu si , en laissant de côté une 
des équations 

P~ °> #Q ~ °> R = o , . . . , 

ou peut déduire des équations conservées une équation 
identique ou contradictoire à celle que l'on a omise , ou, ce 
qui revient au même, si Ton peut assigner une fonction li- 
néaire 

aP + + ■+-••• > 

qui soit identiquement nulle ou qui ne contienne aucune 
des inconnues. Supposons donc que Ton ait 

aP-h PQ -h 7R -h. . .= x, 

on a l'équation identique 

( p + ,„) „ + (/ ,„') _|_ („* + + . . . 

= pV q Q -f- rR H- sS -h 

Si Ton suppose qu'en faisant 

p = otx y 7 — (3.r, r= 7 x,.. , 

les fonctions \>, *>', y", etc., deviennent respectivement 

— w-flx, — /iT-j-Vx, — m" Y'x, . . , 

on aura l'équation identique 

(X' + V'-hX"'-*-. ..)** — ()»+V«' +lV+...)*=xï f 
et, par suite, 

V + V J + V' , +... = o, x -h X//i -h \'m' -h X"m" -f-. . . = o y 
d'où résulte 

X — o, X'=o, X" = o,..., 

et, par suite, 

/ = 0; 
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c'est-à-dire que les fonctions V, V, V", etc., devraient ne 
pas changer si p, q, r, s, etc., reçoivent des accroisse- 
ments quelconques proportionnels aux nombres a , (3 , y, 
à, etc. Un pareil cas, dans lequel la détermination des in- 
connues ne serait pas possible, même si l'on donnait les 
vraies valeurs des fonctions , n'appartient pas à notre sujet, 
comme nous en avons averti plus haut. 

Au reste, on peut réduire facilement tous les cas à celui 
où les fonctions V, V, V", etc., sont linéaires. Désignons 
par 7r, p, a y etc., des valeurs approchées des inconnues 
p, <7, r, 5, etc. (que nous obtiendrons en faisant usage de v 
équations prises parmi les p. équations 

V=M> V' = M', V":=MV..), 

et posons 

p = n 4-//, <] = ■/ + r/, r = p + r', * = * + ..; 

il qst clair que ces nouvelles inconnues seront si petites, que 
leurs carrés et leurs produits seront négligeables, et que les 
équations deviendront linéaires par suite des substitutions 
indiquées. Que si, à la fin du calcul, on trouve contre 
toute attente que les valeurs de p\ ç\ r', s', etc., qu'on 
en tire soient trop considérables, et qu'il paraisse peu sui- 
de négliger leurs carrés et leurs produits, on remédiera 
à cet inconvénient en répétant la môme opération (mais 
en prenant pour 7:, ^, p, 7, etc., les valeurs corrigées de 
p, <7, r, 5, etc.). 

7. 

Tant qu'on n'a qu'une seule inconnue p , pour la déter- 
mination de laquelle on a trouvé que les fonctions 

ap ~\-n, a'p ri, a" p -+- n" , . . . , 

prenaient respectivement les valeurs 

M , M', M", . . , 

» 

et cela par des observations également exactes , la valeur la 
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plus probable de p est 



am -h a' m' -h -f- . . . 

a'+^+a" 1 -h. . . 



A = 
en posant 

m = M — n , m' = M' — /«" = M" — 



Pour apprécier le degré de précision qu'on doit attri- 
buer à cette valeur, supposons que la probabilité d'une 
erreur A, commise dans les observations, soit exprimée 
par 

h _/,t A t 

il en résultera que la probabilité que la vraie valeur de p 
soit A+p', sera proportionnelle à la fonction 

e -h*[(ap- my +• («> - m')* -f- (a" p - «-)» +. . . J 

dans laquelle on aura fait 

= A -h/?'. 

L'exposant de cette fonction peut être réduit à la forme 
— h* (a 7 + a" 4- a"' H )(/>*— 2/?A-f-B), 

dans laquelle B désigne une quantité indépendante de p : la 
fonction sera par suite proportionnelle à 

On voit que le degré de précision qu'il faut attribuer à la 
valeur de A est le même que si cette valeur avait été 
trouvée par une observation immédiate dont la précision 
serait à la précision des observations primitives comme 



h \Ja- -h «" 4- n'" 1 -h . . . est à /i, 

ou comme 

y/a» -4- n ,J a"* -f- . . . est à i 
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8. 

Avantde rechercher, dans le cas de plusieurs inconnues , 
le degré de précision qu'on doit attribuer à chacune d'elles , 
il importe d'étudier plus attentivement la fonction 

*»" H- v" i H- . • • , 

que nous désignerons par W. 

I. Posons 

et 

P ,x 

W — — = W, 
a 

il en résulte évidemment 

p' = Pî 

comme on a 

a. dp * 

on voit que la fonction Wsera indépendante de/?. Le coeffi- 
cient 

sera toujours évidemment une quantité positive. 

II. De même, posons , 

i dW 

1 dq 

et 



on aura 

1 dq a dq a 

et 

dW 

— t— -o. 

Donc la fonction W" est indépendante à la fois de p et 
de q. Ces circonstances n'auraient plus lieu si Ton pouvait 
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avoir 

f' = o. 

Mais il est clair que W se déduit de 

^ + c' , + c" , -f-.,. , 

en remplaçant dans y, v' , v h etc., la quantité p par sa va- 
leur tirée de l'équation 

// = o : 

donc (3' sera la somme des coefficients de q* dans v % 

etc., après cette substitution. Mais ces coefficients sont 
tous des carrés et ne peuvent s'évanouir tous à la fois, si ce 
n'est dans le cas , que nous excluons de nos recherches , où 
les inconnues seraient indéterminées 5 donc jS' doit être po- 
sitive. 

III. Si l'on pose , enfin , 

1 </W" 



2 dr 

et 



r 1 



AV" « -, = W", 



on aura 

a p 

et W" sera indépendant de yu, de (7 et de r. On prouvera 
comme plus haut que le coefficient 7" doit être positif. On 
voit, en effet, facilement que f est la somme des coeffi- 
cients de r* dans y*, y'*, v"% etc., après que les quantités // 
et q ont été éliminées de v', v h \ etc., à l'aide des équa- 
tions 

p' = 0, q' — o. 
IV. De la même manière en posant 

1 riVf* s'* 
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on aura 

a P 7 

W" sera indépendant de p, ^, r et 5, et ô m une quantité 
positive. 

V. S'il y a un plus grand nombre d'inconnues, on conti- 
nuera de la même manière et l'on aura enfin 

W = l -p" -h ?" H- y, r" ■+■ — j" + . . . -4- CODSt., 

expression où a, (3', y", d f ", etc., désignent des quantités 
positives. 

VI. On a déjà vu que la probabilité d'un système de va- 
leurs de p, q, r, $, etc., était proportionnelle à la fonction 

<? — : par conséquent, la valeur de p restant indétermi- 
née, la probabilité d'un certain système de valeurs de (7, 
r, 5, etc., sera proportionnelle à l'intégrale 



x 



00 -A'W 



qui est égale, d'après le théorème de Lapl ace , à 

et cette probabilité sera proportionnelle à la fonction 

e 

- 

De même, si Ton considère de plus q comme indéterminé, 
la probabilité d'un système de valeurs de r, s, etc., sera 
proportionnelle à 



x 



*e-*' w 'rf ït 
— co 



c'est-à-dire à 

/l 'jj' 'ttV V/' ^ /> 
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et, par suite, proportionnelle à e~"* IW \ De même si r 
est aussi regardée comme indéterminée, la probabilité 
d'un système de valeurs déterminées de s, etc., sera pro- 
portionnelle à e~ h * " , et ainsi de suite. Supposons que 
le nombre des inconnues se réduise à quatre 5 les conclu- 
sions seraient les mêmes dans le cas général. La valeur la 
plus probable de s sera 



r 

"F' 

et la probabilité qu'elle différera de a de la véritable va- 
leur sera proportionnelle à 




r m 

_ y 



d'où nous concluons que \f¥' mesure la précision relative 
à cette détermination , en prenant pour unité la précision 
des observations primitives. 

9. 

Par la méthode du paragraphe précédent un certain 
degré de précision a été assigné à la seule inconnue qui, 
dans le travail de l'élimination , a été gardée la dernière. 
Pour éviter cet inconvénient, nous allons calculer d w d'une 
autre manière. 

Des équations 




p= P >, 



o. 



1 '■pr v w 



. ( .3o ) 

on tire, en les résolvant par rapport à p', q' y r\ s', 
/>'=P, 

/ = S + ©"R •+- Uï>" Q 4- X"P, 

de sorte que .A , A', *l>", ui>', ifl>", sont de# quantités 
déterminées. On aura donc (en restreignant à quatre le 
nombre des inconnues) • 

r x" e" i 

d'où résulte la conséquence suivante : Les valeurs des in- 
connues p, q, r, s, etc., que Ton doit tirer des équations 

P=o, Q^-o, R = o, S = o, 

sont évidemment exprimées par des fonctions linéaires de 
P, Q, R, S, etc., savoir : 

/> = L -h A P-f-B Q + C R-f-D S -h ... , 
7 =L'+A'P + B'Q+C'R + D'S + l .., 
r = L" -f- A" P + B" Q -h C" R D" S -f- . . . , 
/ = L w + A m P -h B w Q -h C w R -f- D w S + . . . , 



Cela posé, les valeurs les plus probables de ces incon- 
nues sont respectivement L, L', L", etc. Les degrés de 
précision qui doivent être attribués à ces déterminations 
sont respectivement 

i i i i 

vTÂ' y'B 7 ' j<?' W'"* 

en prenant pour unité la précision des observations pri- 
mitives ; car ce que nous avons dit plus baut de l'incon- 
nue s ^pour laquelle D'" répond à ~\ s'applique aux au- 
tres inconnues par une simple permutation. 
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Pour éclaircir par un exemple les recherches qui pré- 
cèdent, supposons que par des observations pour les- 
quelles une égale précision doit être présumée, on ait 
trouvé 

p — q ■+• 2r = 3, 
3p 4- i(f — 5r = 5, 
kf 4- q -h 4 r = 21 > 

mais que, par une observation à laquelle une précision 
égale à - doit être attribuée, on ait trouvé 

— 2 p 4- 67 -h 6r = 28. 
A cette dernière nous substituerons la suivante , 

— p -h 3q -\-'6r = 

que nous supposerons provenir d'une observation aussi 
précise que les premières. De là on tire 

P =27/? 4- 6q — 88, 

Q = 6p -f- l5q 4- r — 70, 
R = 94-54^—107, 

et par l'élimination 

19899/? — 49*54 4-809P— • 324 Q 4-6R, 
7377= 2617— 12P4- 54Q — R, 
6633 r = 12707 4- 2P — 9Q 4- I23R. 

Les valeurs les plus probables des inconnues seront donc 

p=z 2,470 y = 3,55i, r= 1,916, 

avec des degrés de précision égaux respectivement à 

9 
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11. 

Le sujet que nous avons traité jusqu'ici donnerait lieu 
à d'élégantes recherches, auxquelles nous ne nous arrête- 
rons pas, pour ne pas trop nous écarter de notre objet 
principal. Par la même raison, nous réservons pour une 
autre occasion l'exposition des artifices qui permettent de 
réduire le calcul à un algorithme plus expéditif. Qu'on 
nous permette seulement d'ajouter une seule observation. 

Lorsque le nombre des fonctions ou des équations pro- 
posées est considérable , le calcul est surtout rendu pénible 
par cette circonstance, que les coefficients par lesquels on 
doit multiplier les équations primitives pour obtenir P, Q, 
R, S, etc., sont presque toujours des fractions décimales 
compliquées. Si l'on ne croit pas important dans ce cas de 
calculer ces produits avec le plus grand soin à l'aide des 
Tables de logarithmes, il suffira le plus souvent de leur 
substituer des nombres plus simples qui eu diffèrent peu. 
Il ne peut en résulter d'erreurs notables qu'autant que la 
précision des inconnues devient moindre que la précision 
des observations primitives. 

12. 

Au reste , le principe d'après lequel la somme des carrés 
des différences entre les quantités observées et les quan- 
tités calculées doit être un minimum, peut encore s'éta- 
blir sans recourir au caleul des probabilités, comme il suit. 

Lorsque le nombre des inconnues est égal au nombre 
des observations, on peut déterminer les premières de 
manière qu'elles satisfassent aux secondes. Mais lorsque 
le premier nombre est le plus petit des deux , on ne peut 
obtenir un accord absolu lorsque les observations ne sont 
pas douées d'une précision absolue. Il faut donc dans ce 
cas chercher à établir l'accord le plus satisfaisant, c'est-à- 
dire à faire en sorte que les différences soient atténuées 
le plus possible. Mais cette idée a par elle-même quelque 



Digitized by 



( «33 ) 

chose de vague. En cÛèt, quoiqu'un système de valeurs des 
inconnues doive être sans aueun doute préféré à un» autre 
système où toutes ces différences seraient respeetivement 
plus grandes, le ehoix entre deux systèmes dans l'un des- 
quels l'accord serait plus satisfaisant pour quelques-unes 
des observations, mais moins satisfaisant pour d'autres, est 
en quelque sorte arbitraire, et Ton peut évidemment pro- 
poser plusieurs principes par lesquels la première condi- 
tion soit remplie. En désignant par A. A', A", etc., les 
différences entre le calcul et les observations, ou satisfera à 
cette condition, non-s< ulcmcnt si A*-f- A'*-f- A" s 4-. • 
devient un minimum (ce qui est notre principe), mais en- 
core si A v -h A'* 4- A"*h-. . ., ou A 4 -h A /6 -+- A"«H-..., 
ou généralement une somme de puissances paires, devient 
un minimum. Mais de tous ces principes le notre est le plus 
simple, tous les autres nous entraînant dans des calculs 
extrêmement compliqués. Au reste, ce principe, dont nous 
avons fait usage dès l'année 17ÇP, a été donné dernière- 
ment par Legendrc dans ses Nouvelles méthodes pour la 
détermination des orbites des comètes, Paris, 1806; on 
trouvera dans cet ouvrage plusieurs conséquences que le 
désir d'abréger nous a fait omettre. 

Si l'exposant de la puissance paire dont nous venons de 
parler était infini, nous serions ramené au système dans 
lequel les plus grandes erreurs sont moindres que dans tout 
autre système. 

Laplace se sert, pour la résolution d'équations linéaires 
en nombre plus grand que les inconnues, d'un autre prin- 
cipe, proposé d'abord par Boscovich, savoir, que la somme 
des valeurs absolues des différences devienne minimum. On 
peut facilement démontrer que le système des valeurs des 
inconnues trouvé par ce seul principe doit Nécessaire- 
ment (*) satisfaire à autant d'équations, prises parmi les 



(*) Excepté quelques cas spéciaux où il y a indétermination. 
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proposées, qu'il y a d'inconnues, de sorte que les autres 
équations ne sont employées que pour décider le choix que 
Ton doit faire. 

Si , par exemple, l'équation V = M est du nombre de 
celles qui ne sont pas satisfaites, le système des valeurs 
trouvées par le principe en question ne serait pas altéré si 
au lieu de M on avait observé une autre valeur N telle, que, 
n étant la valeur calculée, les différences M — n et 3N — n 
fussent de même signe. Au reste, Laplace tempère en 
quelque sorte ce principe en y ajoutant cette nouvelle con- 
dition, que la somme des différences, prise avec leurs 
signes, soit nulle. Il en résulte que le nombre des équations 
satisfaites est- moindre d'une unité que le nombre des in- 
connues : mais l'observation que nous venons de faire sub- 
siste encore lorsqu'il n'y a que deux inconnues. 

NOTE IL 

APPLICATION DE LA MÉTHODE DES MOINDRES CARRÉS A LA 
CORRECTION DES ÉLÉMENTS DE LA PLANÈTE PALLAS. 



1. 

M. Gauss a donné, dans le tome 1 des Mémoires de 
Gottingue, l'application de sa méthode à la correction des 
éléments de la planète Pallas. L'illustre géomètre ayant 
développé sur cet exemple l'algorithme indiqué plus briè- 
vement dans son grand ouvrage Theoria Motus Corporum 
cœlestium (uoi'rla Note précédente) , nous avons cru devoir 
traduire ici cette portion de son Mémoire. La première 
partie exigeant la connaissance approfondie de la théorie 
du mouvement des planètes, nous nous dispenserons de la 
reproduire, et nous prendrons pour point de départ les 
douze équations auxquelles les corrections des six éléments 
de l'orbite doivent satisfaire. 
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Eu désignant ces corrections par 

dL, d7, diz, df, dQ, di, 

les équations obtenues par M. Gauss sont les suivantes : 

o = - . 83",93 -f o , 79 363 d L 4- 1 43 ,66 d7 -h o , 3o4 9 3 d n 

4- 0,95920 d f — o,i8856dQ -h 0,17387 d/; 

o = — 6",8i ~o,o2658diL-h46,7id7 4-o,o?.658d7r 

— o,2o858d? 4-0,15946 d$ -+- 1, 25782 d/; 

0=3— o",o64-o,5888od J L4- 358,i2d7 4-o,262o8d7r 

— o, 85234 d<p 4- o, 1491 2 d Q 4- o, 17775 d/ ; 

o = — 3",094-o,oi3i8dL-f-28,39d7 — o,oi3i8d7r 

— 0,07861 d<p + 0,91704 d Q o,54365d/; 
o = — o",o2 4- i,73436dZ,+ 1846, 17 d7 — o, 546o3 d;r 

— 2,o5662 d<p — o, i8833 d Q — o, i7445d/; 

o = — — °3 i26o6dZ. — 227 ,42 d7 4- o, 12606 d n 

— o,38q39<] cp H-o,i7i76d Q — 1 , 35.(4 1 dt; 
o = — 2",3i -ho,99584dZ-t- i579,o3d7-f- o,o6456d;r 

4- 1 ,99545 d ? — o,o6o4od$ — o,33 7 5od/; 

o = -f- 2",47 — 0,08089 d L — 67 ,22d7 -f- o, 08089 dn 

— o,0997od<p — o, 4635g d Q 4- 1 ,22803 di ; 

0=4- o",oi 4-o,653n dZ,4- 1 329,09 d7 4- o, 38994 drc 

— o, 08439 dcp — o,o43o5 d Q 4- o,34268d< ; 

0 = 4- 38", 12 — o,oo2i8dZ+ 38,47 d7 -4-0, 00218 dir 

— o,i8 7 iod.p -+-0,47301 dQ — 1,14371 di; 
0 = — 3i 7 ",73h- o ,69957 d L 4- 1719,32 d7 -+-o,i2 9 i3d7r 

— 1 ,38787 d<j>H-o,i7i3od£J — o,o836o di ; 

0=4-1 17",97 — o,oi3i5 dL — 43>84 ^7 4- o,oi3i5 dn 

4- 0,02929 d<p H- 1 ,o2i38d Q — 0,27187 di. 

D'après la nature des observations qui ont fourni la 
dixième de ces équations, elle inspire trop peu de confiance 
pour qu'on juge utile de la faire intervenir, et c'est d'après 
les onze autres seulement que l'on déterminera les six in- 
connues. 
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Les explications suivantes sont littéralement traduites du 
Mémoire de M. Gauss. 



2. 

Dans l'impossibilité où nous sommes de satisfaire exac- 
tement aux onze équations proposées , c'est-à-dire d'an- 
nuler tous les seconds membres , nous chercherons à rendre 
la somme de leurs carrés aussi petite que possible. 

On aperçoit facilement que si Ton considère les fonc- 
tions linéaires 

n ap -\- bq -+- cr -h ds -+-... = w, 
ri -+- a'p -h b'q +Jr -\-dfs . . = w', 
„"+ a " p + b"q+ c"r+ d" s H- . . . = 



les équations qu'il faut résoudre pour rendre 

un minimum, sont 

aw -f- a! -f- a" vS' -f- . . . = o , 
bw -+■ b'w'+ b"w" + ... — o, 
cw -h d a/ c" w" . . . = o , 



ou, en posant, pour abréger, 

an -+- a' ri -f- a" n" -f- . . . = (an ) , 
à 1 -W -f-a'" -»-.. . = (aa), 
ab +.a'b' + a"b"-h... = (ab), 



b*+V* -h*"' + = (bb), 

bc-\-b'c'-^ b"c"-+- = {bc), 



p y </, r, s, etc., devront se déterminer par les équations suî- 
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vantes : 

♦ 

(an) -h (aa)p -\~ (ab)q +(flc)r + .,. = o, 
(bn) -h (ab)p ~+-(bb)q (bc) r . . = o, 
(<vi) -\-(ac)p + (bc)q -\-(cc) r-*-. . , = o, 

L'élimination, très-pénible lorsque le nombre des incon- 
nues est considérable, peut se simplifier notablement de la 
manière suivante. Outre les coefficients (an), (aa),etc, 

[dont le nombre est - (P-h 3i) , si le nombre des incon- 

nues est i], supposons que Ton ait calculé la somme 

/i 1 «'» H- n" 3 -h . . . = (nn) ; 

on voit facilement que Ton a 

a = (nn) -t- z(an)p n(bn) q -h 2(c/i) r -f-. . . 
(aa)p 3 -f- z(ab)pq -+- 2 (<?<?)/>/• -h. . . 

-h (ce) r 1 -H 2 (c</) rj H- . . . ;* 

et, en désignant 

(an)-\- (aa) p -h (ab) q-{- . . . 

par A , tous les termes de -, r qui contiennent le facteur p, 

r \ aa ) 

se trouvent dans l'expression ft, et, par suite, 

A' 

(aa) 

est une fonction indépendante de p. C'est pourquoi, en 
posant 

(nn ) — )— ^~(nn t i , 
/ 1 x (an)(bn) 

w — (airr= {bn ^ 

I \ ( an ) i cn ) I \ 
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(«fl,i) + 2(è/»,i)7 + 2(c»,i)r+ 2(rf«,i)5 .. 

+ (èi,i) 9 1 + a(^, i)?r + 2(M, i)<jr* 
H- (ce, i) r 1 -»- 2 (c<tf, i) 



nous désignerons cette fonction par Cl'. 
De même, en posant 

(**,i)-h(M, i)f + (Ac,i)r ..= B, 
la différence 



sera indépendante de ^; nous la représenterons par &" . 
En posant , de même , 

et 

(c«, 2) -h (ce, 2) r -h (rrf , 2) j H- . . . = C, 
la différence 

(«r, 2) 

sera une fonction indépendante de r. 
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En continuant ainsi, nous formerons une suite d'expres- 
sions 12, Cl\ Cl" , etc., dont la dernière sera indépendante 
des diverses inconnues, et représentée par (nn,[i), si fi 
désigne le nombre de ces inconnues ; nous aurons alors 

A ' B ' C» DI . ' \ 

(aa) (bb,i) (cc,i) (dd,ô) 

On prouvera facilement que Cl étant une somme de carrés 

w J -h w" -H w" 3 4- ■ 

et ne pouvant devenir négative, les diviseurs (aa), (bb, i), 
(ce, a), etc., sont tous positifs. (Nous supprimons, pour 
abréger, le détail de la démonstration. ) D'après cela , la va- 
leur minimum de Cl correspond évidemment aux valeurs 
des inconnues, pour lesquelles 

A = o , B — o , C — o , . . . , 

et., en commençant à résoudre le système par la dernière 
équation, qui ne contient qu'une inconnue, on trouvera 
les valeurs de p, q, r, s, etc., sans avoir aucune élimination 
à effectuer. La méthode donne, en même temps, la valeur 
minimum de Cl, qui est (nn , p). 



3. 

Appliquons ces principes à notre exemple , dans lequel p, 
g, r, 5, etc., sont remplacés par dZ, d7 , drc, d<p, d Q , dz. 
J'ai trouvé, par des calculs exécutés avec soin : 



(nn) 




i 48848, 


(ae) 




— 0,09344, ■ 


(an) 




— 371,09, 


(nd) 




— 2,285i6, 


(bn) 




— 58oio4, 


(ae) 




— 0 , 34664 f 


(en) 




— 1 1 3 , 45 , 






— 0,18194, 


(dn) 




+ 268 , 53 , 


(bb) 




H- lo834225, 


(en) 




-1-94,26, 


(bc) 




— 49,o6, 


ifn) 




— 3i,8i, 


(bd) 




— 3229,77, 


(an) 




+ 5,91569, 


(be) 




— i9 8 » 6 4, 


(ab) 




-h 7203,91 


(¥) 




— i43, o5, 
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(ce) 


— 


4-0,71917 


(cd) 




4- 1 , 13382, 


(ce) 


: 


4- 0,06400, 


(</) 


= 


4- 0 , 2634 1 , 


(dd) 




-h 1 2,oo34o , 


[de) 




— 0,37137, 


(df) 




— 0, I 1762, 


(ce) 




4- 2,28215, 


(ef) 




— o,36i36, 


( ff\ 




-h 5,62456. 


où I on 


déduit : 




— . 


-h 125569, 


(bn,i) 


__ 


— i38534, 


(c/i,i) 


— 


— 1 *9>3i , 


(dn,i) 




— 125, 18, 


{en,\) 




+ 7 2 >52, 


0,0 




— 43,22, 


(bb,i) 




4- 2458225, 


(bc,i) 




H-62,i3, 


(bd,i) 




— 5io,58, 


(be,x) 




4- 2i3,84, 


(*/,0 




4- 73,45, 


(ce, i) 






(cd,i) 




4- 1 ,09773, 


{cc,i) 




— o,o5852 , 


(</>0 




4- 0 , 26054 , 


(dd,J) 




-h 1 1 , 1 2064 , 


(de,i) 




— o,5o528, 






— 0,18790, 






-h 2,26185, 






— 0,37,202, 


(//,0 




-h 5,6i9o5. 



De la même manière : 

(nn 9 2) = 4- 1 17763, 
(en, 2) = — 1 i5,8i , 
(dn y 2) — — i53,95, 
(e/1,2) = -h 84, 5 7 , 



o ) 



0,2) 


— 


— 39,o3, 


(CC 9 7.) 


• — • 


-h 0,7 161 2 , 


(«/,*) 




-f- 1 , 1 io63, 


(«?, 2) 


= 


— 0,06392, 


('/,*) 




+•0,2.5868, 






4- 1 1 ,01466 


V ac > -'•J 






(<*/,») 




— 0, 17265, 


2) 





4- 2 , 2432.5. 


('/>*) 
V ' / 





— 0,3*7841 , 


(//,2Ï 





4- 5,6i686. 


u ou : 






(/?/*, 3) 


— 


H- 99 o3 4, 


(''«, 3) 


: — ■ 


4- 25,66, 


(<-//, 3) 


: — 


4- 74 y ?-3 , 


e/«,3) 


= 


4- 2,75, 


(rW, 3) 





-+-9,29213, 






— O, OOI 70 , 


MA 3) 




0,57384, 


(et-, 3) 


— ■ 


4- 2 ,23754. 


(</,3) 


; 


— 0, 35632, 


(>7,3) 





4- 5,52342. 


Lie même : 




V 




— 


4- 90903, 


(on ft\ 




4- 7^,20, 


(>,4) 




4-4,33, 


v*-** » 




-4- 2.29 3/l(> 


(</> 4) 





— 0,37766, 


(//, 4) 


— 


-h 5,48798. 


D'où : 






(un, 5) 




4-96418, 


(/«, 5) 




-h I7,M, 


(/A 5) 




4- 5,42383. 


D'où enfin 


• 
• 




(«*, 6) 




4-96364. 
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Nous avons donc les six équations suivantes : 

o = -f- 1 7", 1 1 5 ,42383 di, 

o = -h 75", 2.3 4- 2, 2?.346 d Q — o, 37766 d/, 

o r= -+- ?5",66 -+- 9,2921 3 d? — o,36i75dQ — o,57384d/, 

o — — 1 i5",8i -f- 0,71612 dit -f-i,no63dç — 0,06392 dQ 
4- o, 25868 d/, 

o = — i3854" 4- 2458225 d 7 -4-62, i3dw ~o,5io58d<p 
4-2i3,84dQ 4-73, 45d/, 

0 = — 37i",o9-h5,9i569dZ, -4- 7203,91 d7 — o,oo344d 7r 
— 2.2o5i6d«p — o,34664d Q— o, i8i94d/; 

d'où l'on déduit : 

di = — 3",i5; 

'dQ = - 34",3 7 ; 

d<p = — 4%2 9 ; 

d ?r = -4- i66",44 » 

d? =+ o ff ,o54335; 

d L — — 3",o6. 

Telles sont les corrections qu'il faut apporter aux élé- 
ments trouvés d'abord pour la planète. 



NOTE III. 

MÉMOIRE SUR LA DÉTERMINATION DE LA PRÉCISION 

DES OBSERVATIONS, 

PAR M. GAtJSS. 

( Extrait du Zeitschrijt fur Astronomie und Verwandte Wissenschaften , 

tome I , page i85.) 

i. 

Pour établir les principes de la méthode des moindres 
carrés, nous avons admis que la probabilité d'une erreur 
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d'observation A soit exprimée par la formule 

h —h* A 1 

où tt représente la demi -circonférence, e la base des loga- 
rithmes hyperboliques, et h une constante que l'on peut con- 
sidérer [Thcoria Motus Corporum cœlestium, art. 178 (*)] 
comme la mesure de l'exactitude des observations. Il n'est 
pas nécessaire de connaître la valeur de h pour déterminer, 
à l'aide de la méthode des moindres carrés, les valeurs les 
plus probables des quantités dont les observations dépen- 
dent, le rapport de l'exactitude des résultats à l'exactitude 
des observations est également indépendant de h. 

Toutefois , comme la connaissance de la quantité h est 
très-intéressante et instructive, je vais montrer comment 
les observations peuvent servir à la déterminer. 



Commençons par quelques remarques qui éclairciront 
la question, et représentons par 0 (t) l'intégrale définie 



Quelques valeurs particulières de cette fonction donneront 
une idée de sa marche : 





0,47^0,363 




P 


et = 


o,5 




0 , 5g5 1 1 6 1 




0 X 1 ,247 790 


®e = 


0,6 


t = 


0 ,732 8691 




p X 1 ,5366i8 


et = 


°'»7 


t = 


0,906 1939 




p X 1 1900 o32 


et = 


0 , 8 




1 




p X 2,096716 


et = 


0 , 842 7 


t = 


1 , 1630872 




px 2,348664 




°>9 


/ =r 


1 ,821 3864 




p x 3,818930 


0/ = 


°>99 




2,327 6754 




pX 4,88o 4 7 5 


et = 


°>999 




2,751 o654 




0x5,768 204 


et =7: 


°>9999 








00 


et = 


1 



(*) Voyet pape 1 19 de ce rolurae. J. B. 



Digitized by Google 



( i43 ) 

La probabilité que l'erreur dune observation soit com- 
prise dans les limites 4- A et — A-, ou, abstraction faite 
des signes , qu'elle ne surpasse pas A , sera égale à 

elle sera le double de l'intégrale en question , lorsque celle- 
ci est prise depuis x = o jusqu'à x = A, et, par consé- 
quent, elle sera égale à 0 (h A). 

Ainsi la probabilité que Terreur ne soit pas moindre que 

| est égale à ~, ou égale à la probabilité du cas contraire. 

J'appellerai donc cette quantité | V erreur probable, et je 

la désignerai par r. 

Au contraire , la probabilité que l'erreur excède 
rx 2,43866*4 n'est que de — ; la probabilité que Terreur 
surpasse r X 3, 81 8 390 n'est que de ^ ; et ainsi de suite. 

3. 

Supposons maintenant que les erreurs de m observations 
réellement faites soient a, /3 , 7, etc., et cherchons quelles 
conséquences on en peut tirer relativement aux valeurs 
de A et r. 

En faisant deux hypothèses sur la valeur exacte de A, et 
la supposant égale à H ou égale à H', les probabilités qu'elles 
soient entachées des erreurs a, (3 , 7, etc., seront, pour les 
deux cas , dans le rapport de 

à 

H' e -"'•« , X H'»-"""' X H' e - H 'V X . ., 
c'est-à-dire comme 

Il est évident que les probabilités que H ou H' soient les 
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véritables valeurs de h, sont dans le même rapport [Théo- 
ria Motus Corporum cœlestium, art. 176 (*) par consé- 
quent la probabilité d'une valeur quelconque de h est pro- 
portionnelle à 

et la valeur de h la plus probable est celle pour laquelle 
cette fonction devient un maximum. Mais on trouve par 
les règles connues que h est alors égal à 



m 



a (a' -h p»4- 7 » + . . .) ' 
donc la valeur de r la plus probable sera alors 



> 



2 (a J -f- p a + 7 7 + • • •) 



ou 



0,6744897 X y/^( a ' -t-P'-H 7'-+-- • •)• 
Ce résultat est général, que m soit grand ou petit. 

4. 

Il est facile de comprendre que les valeurs trouvées 
pour h et pour r sont d'autant moins certaines que le nom- 
bre m est plus petit. 

Développons maintenant le degré d'exactitude que Ton 
doit attribuer aux valeurs de h et de r lorsque m est un 
nombre considérable. 

1 

Désignons par H la valeur de h la plus probable que nous 
avons trouvée , 



m 



et remarquons que la probabilité que H soit la véritable 
valeur de //, est à la probabilité que (H -f- X) soit cette véri- 

( 11 ) Voree page 1 16 de ce volume. J. B. 
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table valeur dans le rapport de 



ou comme 



Le second terme ne sera sensible par rapport au premier 
que si ^ est une petite fraction , et, dans ce cas , nous pour- 
rons remplacer le rapport indiqué par 

/'m 

H* * 
lie 11 

Ce qui veut dire : La probabilité que la valeur véritable 
de h soit comprise entre (H -i- 1) et (H -h A dl) est ap- 
proximativement égale à 

Ke »' d\, 
où K est une constante telle, que l'intégrale 



/ 



Ke H * dl 



prise entre les limites admissibles de ^, devienne égale 
à l'unité. 

Comme dans le cas actuel , à cause de la grande valeur 



de m , e u * devient excessivement petit lorsque ~ cesse 

d'être une petite fraction, il sera permis de prendre l'in- 
tégrale depuis — oo jusqu'à -h oo , et l'on obtient 



H V 7T 



Par conséquent, la probabilité que la véritable valeur de h 
g. io 
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soit comprise entre (H — A) et (H -f- sera égale à 



I ^ 

e t sera égale à - lorsque — \m = jo- 
li y a donc un contre un à parier que la véritable valeur 
de h soit entre 

U ( \ ^=r\ et h( H--4r^ , 

ou que la véritable valeur de r soit entre 

R R 

et , 



où R désigne la valeur la plus probable de r trouvée dans 
l'article précédent. Ces limites peuvent s'appeler les limites 
probables des véritables valeurs de h et de r. 11 est évi- 
dent que nous pourrons admettre ici pour limites probables 

»(-*)-»(■**)• 



5. 

Dans la discussion précédente, nous avons considéré a, 
(3, y, etc., comme des quantités définies et données, afin d'é- 
valuer la probabilité que la véritable valeur de h ou de r 
soit comprise entre certaines limites. 

On peut envisager la question sous un autre point de 
vue, en admettant que les erreurs des observations soient 
soumises à une loi déterminée de probabilité; on peut alors 
évaluer la probabilité pour que la somme des carrés de m 
erreurs d'observations tombe entre certaines limites. La- 
place a déjà résolu ce problème dans le cas où m est un 
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nombre très-grand, ainsi que le problème de déterminer 
la probabilité que la somme de m erreurs d'observations 
tombe entre certaines limites. 

Il est facile de généraliser cette recherche ; je me bor- 
nerai à indiquer ici le résultat. 

Désignons par 9 (x) la probabilité d'une erreur d'obser- 
vation x , de manière que 



r 



y x . dx ~ 1 . 



Désignons encore par K„ la valeur de l'intégrale 

<p x . x n dx. 



L 



Soit ensuite 

S„=a»+P"-l-7 n -|- .., 

où a, (3, 7, etc., représentent m erreurs quelconques d'ob- 
servation j les termes de cette somme seront tous pris posi- 
tivement, même si n est impair. 

mK n sera alors la valeur la plus probable de S„, et la 
probabilité que la véritable valeur de S„ tombe entre les 
limites (mK„ — X) et (mK rt -HA) sera égale à 

X 

par conséquent , les limites probables de S„ seront 



mK n — p s/zm(K, n — K' ) 

Ce résultat s'applique, d'une manière générale, à toute 
loi des erreurs d'observation. En l'appliquant au cas parti- 
culier où 



<f* = —p= 'G , 



10. 



- 1 ■ 
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nous trouverons 

n •;■(*-!) 

&n — r=z > 

où le signe caractéristique ET est pris dans la signification 
des Disquisitiones générales circa sérient infinitam (*) 
(Comm. no\>. Soc. Gotting., tome III; M. 5, aii. 28.) 
Ainsi 



R = 




K - 1 


K,= 


1 




K,= 


1.3 


1.9 


K 6 = 


1 .3.5 
8/< a ' 


i.*»..3 
K?= 777=' 



par conséquent la valeur la plus probable de S„ sera 

/nn{(fl-i) 
h* \[k 

et les limites les plus probables de la véritable valeur de S„ 
seront 

n — 1 

mil 



A-v^ T p V«L[n + («-i)p "Ml 



et 



/1 — 1 



(*) D'après la notation adoptée par un grand nombre de géomètres, on 
aurait 

n(n) = r(m-i:. 

J. B. 
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Donc, si nous posons comme ci-dessus 



P 

où v représente l'erreur probable <l observation , la valeur 
la plus probable de 




-m 



sera évidemment /•; et les limites probables de la valeur de 
cette quantité seront 



IL 

et 



r) y — - 
n 




//% / /// 



■ I »>- , I /. 



Il y a donc aussi un contre un à parier «pie /' soit compris 
entre les limites 

V •"(•:■■) l V" W'-?' J 

et 

Pour » = n , ces limites seront 



et 




s*? 
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ce qui s'accorde parfaitement avec celles que nous avons 
trouvées art. \. 

En général on aura, si n est pair, les limites 



^'V^//,. i.3. 5. 7 . («-î; 

( «y /»L i .3.5. . (« — i) | ) 



et 



, & d — >— ■ -, ' 

r * Va/ii.3.5.7...(« — i) 

-li t n/ 2 r %+o-(^-^3)...(2/»-ij i) 
x I » V « L r3.5...(n-i) 'Jr 

et si « est impair, les limites suivantes : 



y/ «t.i 2. 3. . . 



i f"' - 3 * 5 -7 - (a* — -| 

( «V »L [2.4.6.... i 



et 



1.2.3 



-m 



/i Ti .3.5.7. ( n — ')* II 



6. 

J'ajoute ici encore les valeurs numériques pour les cas les 
plus simples : 
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Limites probables de r. 

.. 0,84534,3**, .^^Biy/ 

m. o^^x^.^-M); 
iv. o,5„5o I7 xy/|.( i:F 0 v 5 -^ ,86 ) ; 

V. o, 465 553a x ^ | . ( . q= fi^j£»° ) , 

On voil encore, par cetle comparaison , que la deuxième 
manière de déterminer r est la plus avantageuse; car cent 
erreurs d'observation, traitées d'après cette formule, don- 
uent un résultat aussi sur que 

1 14 erreurs d'après la formule I; 

loi) » ». » » III; 

i33 « • » » IV; 

178 V; 

25 1 » » ». » VI. 

Cependant la formule 1 présente l'avantage de se prêter le 
mieux au calcul numérique; et comme son degré d'exacti- 
tude est peu inférieur à celui de la formule II, on peut tou- 
jours s'en servir, à moins que l'on ne connaisse déjà la 
somme des cariés des erreurs ou qu'on désire la connaître. 



Le procédé suivant est encore plus commode, mais beau- 
coup moins exact. 
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Ordonnons par ordre de grandeur les valeurs absolues 
des m erreurs d'observation , et désignons par M le terme 
du milieu si leur nombre est impair, ou la moyenne arith- 
métique entre les deux termes du milieu si leur nombre 
est pair. 

On peut démontrer (ce que nous ne ferons pas ici) que, 
pour un grand nombre d'observations, la valeur la plus 
probable de M estr, et que les limites probables de M 
sont 



et 

ri i-heP 



ou que les limites probables de la valeur de /' sont 



ci 



ou, en nombres, 

m/,_Q,75™974\. 

Ce procédé n T est , par conséquent, pas beaucoup moins 
exact que l'application de la formule VI, et il faudrait 
mettre en compte 249 erreurs d'observation pour obtenir 
un résultat aussi certain qu'en appliquant la formule II à 
cent erreurs d'observation. 

8. 

L'application de ces méthodes aux erreurs commises dans 
48 observations des ascensions droites de l'étoile Polaire, 
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par Bcssel (*), a donné 

S,= 60", 46, 
S 2 = i io",6oo, 
S 3 = a5o" , 34 1 1 1 8. 

De là on a déduit les valeurs les plus probables de /• : 

D'après la formule I. 1 ",065 , Erreur probable = ± o",o68 , 

» II. i",024, » = db 0^,070, 

» III. i",ooi, » =dbo ,/ ,072, 

et d'après l'article 7. i",o45 , » = ±o",i i3 : 

concordance de résultats qu ? on pouvait à peine espérer. 
Bessel donne i",o6j et semble, par conséquent, avoir cal- 
culé d'après la formule I. 

NOTE IV. 

APPLICATION DU CALCUL DES PROBABILITÉS A UN PROBLÈME 

' DE GÉOMÉTRIE PRATIQUE. 

(Extrait d'une Lettre de M. Gaiss à M. Schumacher; Astronomischc 
Nachrichien , tome 1, page 80. ) 



D'après votre désir, je vous envoie les règles relatives à 
l'emploi de la méthode des moindres carrés dans la solution 
du problème suivant : 

Déterminer la position d'un point d'après les angles 
horizontaux observés de ce point entre d'autres points 
exactement connus. 

Cette question, très-élémentaire, ne peut embarrasser 
ceux qui ont bien saisi l'esprit de la méthode des moindres 
carrés. Je développerai néanmoins les formules auxquelles 



(*) Bode , Annuaire astronomique pour 1818, page t3f\. 
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elle conduit, en fa\eur des personnes qui auraient à traiter 
Ja question pratique sans avoir été à même d'étudier cette 
théorie. 

Soient a et b Jes coordonnées de F un des points donnés; 
nous supposerons que Ton compte les x positifs du nord 
au sud et les y positifs de l'ouest à Test; soient x et y les 
coordonnées approchées du point inconnu et dx, dy les 
corrections encore inconnues qu'il faut leur attribuer. Dé- 
terminons deux quantités <p et r par les formules 

b — r a — x b — r 

tang«) = -, r — — = . • , 

n — x cos ? sm <j> 

■ 

r f étant pris dans un tel quadrant, que la valeur de r soit 
positive. 

Posons de plus 

206265" (£— y) 0 2o6265' (a-x) 

a — _— — , p= 

I /azimut du premier point est alors, pour un observateur 
placé au second (en prenant o pour azimut d'une parallèle à 
Taxe des .r), 

les deux derniers termes étant exprimés en secondes. 

Soient <p', a', j3' les quantités analogues à a , /3 et re- 
latives au deuxième -des points donnés, a", fi' celles qui 
se rapportent au troisième, et ainsi de suite. Supposons 
que, pour les mesures angulaires prises au point dont la 
position est inconnue, on ait fait usage d'un théodolite sans 
répétition , dont on dirige successivement la lunette vers 
les points connus , sans changer la place de l'instrument 
lui-même. Si à, h\ h" sont les azimuts observés, on aurait, 
en supposant les observations rigoureusement exactes, et 
r/.r, dy exactement connus, 

— h -f- y.dx +> pdy rrf- // -|- a! dx -f- p' dy 

fdy.... 



(-) 
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Si donc on écrit que trois de ces différences ont la même 
valeur, on trouvera des valeurs approchées pour dx et dy ; 
si Ton n'a observé que trois points, il n'y a rien de plus à 
faire; mais si le nombre des points considérés est plus 
grand, les erreurs seront le mieux compensées, en prenant 
la moyenne de ces diverses expressions (1), égalant à zéro 
la différence entre chacune d'elles et cette moyenne, et 
appliquant à ces équations la méthode des moindres carrés. 

Si toutes les mesures sont indépendantes les unes des 
autres , chacune d'elles fournit une équation entre dx et dy, 
et jl faut combiner ces équations par la méthode des moin- 
dres carrés, en ayant égard , si l'on veut, à l'inégale préci- 
sion des observations. 

Soient, par exemple, i l'angle compris entre le premier et' 
le deuxième point, i' l'angle compris entre le second et le 
troisième, et ainsi de suite, en comptant toujours de gauche 
à droite-, on aura les équations 

? '- f - l >( a '- tf ) ( fc + (P'-p)rf ; = o, 
? "- ? ' - i'H- (a"- *')<lx + ($"- ?)dy = O. 

Si les diverses mesures ont le même poids , on déduira de 
ces équations deux équations normales, en les ajoutant, 
après avoir successivement multiplié chacune d'elles par le 
coefficient de dx ou par celui de dy. 

Si , au contraire, les mesures des angles sont d'une exac- 
titude inégale, et, par exemple, la première soit fondée 
sur n et la seconde sur fi f répétitions , il faut que, dans les 
deux cas, et avant leur addition, les équations soient encore 
multipliées par fx', etc. ; on trouve ensuite dx, dy, etc., 
par l'élimination entre les deux équations normales ainsi 
trouvées. 

(Les préceptes qui précèdent sont seulement destinés aux 
personnes auxquelles la méthode des moindres carrés est 
encore inconnue et pour lesquelles il sera peut-être bon de 
rappeler que, dans les multiplications, les signes de a' — a, 
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p' — 13, etc., doivent être rigoureusement conservés. Enfin, 
je remarque encore que l'on a seulement en vue de com- 
penser les erreurs commises sur les angles , les coordonnées 
des points donnés étant supposées exactes. ) 

Appliquons les indications qui précèdent aux observa- 
tions que nous avons faites en commun sur le Holkensbas- 
tion , à Copenhague. Je dois prévenir que les résultats ne 
peuvent cire ici d'une exactitude rigoureuse. Les points ob- 
servés étant très-rapproebés de la station , une inexactitude 
de quelques dixièmes de pied sur leur position peut exercer 
une influence beaucoup plus grande que les erreurs habi- 
tuellement à craindre sur les angles mesurés. On ne s'éton- 
nera donc pas que la meilleure compensation des angles 
laisse subsister des différences beaucoup plus grandes que 
celles que Ton peut admettre comme possibles dans les ob- 
servations de même nature. Cette application doit être prise 
comme un exemple de la marche à suivre dans d'autres cas. 

Angles mesurés du Holhensbastion. 

Friedrichsberg — Pétri 7 3° 35' 22" , 8 ; 

Pelri — Erlosersthurm 1 o4 . 57 . 33" , o ; 

Erlbsersthurm — Friedrichsberg . 1 8 1 . 27 . 5" ,0 ; 

Friedrichsberg — Fraueuthurm . . 80 . 37 . 10", 8; 

Frauenthurm — Friedrichsthurm. 101 . 1 1 . 5o",8; 

Friedrichsthurm — Friedrichsberg. 178.11. i",5. 

Coordonnées des divers points en pieds de Paris, l'origine étant 
a l'Observatoire de Copenhague, 

Pétri 487 ,7 -h 1007 ,7 ; 

Frauenthurm 710,04- 684,2; 

Friedrichsberg .... 243o ,6 -h 8335 ,0 ; 

Erlosersthurm 2 94° > 0 — 3536 , o ; 

Friedrichsthurm . . . 3o5o, , 3 — 223 1,2. 

Les coordonnées approchées du Bastion sont 

x= +2836,44, 
.r=-f- 444 ,33, 



Digitized by Google 



( «57 ) 

cl nous trouvons ainsi les azimuts : 

Pétri * i66°3o'42",56+ 19,92^4- 83, o^dy, 

Fruenthurm 1 73 . 33 . 5o" , 54 -h 10 , 80 dx 4- 95 , 78 dy ; 

Friedrichsberg 92 . 56 . 29" ,46 4- 26 , 07 dx 4- 1 , 34 dy ; 

Erlosersthurm 271 .29. 25", 38 — 5 1,791b; — i y 35dy; 

Friedrichsthurm . . . 274. 45. 4 1 ">4^ — 76,56<te — 6,38^j. 

L'angle sous lequel on voit la distance de Pétri à Friedrichs- 
berg est, par suite, 

7 3° 34' 3" , 1 o — 6, 1 5 dx 4- 8 1 . 7 o dy ; 
en l'égalant à l'angle observé, on a 

~~ 79" >7° ~~ 6 » »5 dx 4- 81 ,70 dy = o. 

On obtient, de la môme manière, les équations sui- 
vantes : 

69", 82 — 71,71 rte — 84,39^=0, 
9", 08 -h 77,86<te-h 2,69^7 = 0, 
o",a8 — 15,27 dx 4- g4»44'(r ==0 9 
o",o4 — 87,36<te — 102,1 6 tfj =5 o, 
— 3",4 2 4- 102, 63 «te 7,72^ = 0. 

En supposant les observations également précises, on déduit 
de là les équations normales 

29640 rfr + 1 4o33 dy = 4 1 68", 
i4o33<te 4- 33a 19^ = 12383", 

et, par suite, 

dx — — o,o5, r//=o,4o; 
les coordonnées du Bastion sont donc 

2836,39, 
444,73. 

Les différences entre les valeurs observées des angles et 
celles que Ton calculerait d'après ces résultats sont trop 
grandes pour qu'on puisse les attribuer aux erreurs d'ob- 
servation-, elles tiennent, comme nous l'avons dit, à un dé- 
faut de précision dans la déterminalion des points connus. 
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Les coordonnées x et y, adoptées comme premières ap- 
proximations, ont été déduites directement du quatrième 
et du cinquième angle. Quoique la méthode directe doive 
être considérée comme un sujet presque épuisé, j'indiquerai 
néanmoins , pour compléter, la méthode que j'ai l'habitude 
d'employer en pareil cas. 

Soient a et b les coordonnées du premier point connu; 
celles du second seront de la forme 

o + RcosE, b -+- R sin E , 

et celles du troisième, 

«4-R'cosE', A-f-R'sinE'. 

Soient 

a -+- p cos « , b -+- p sin i 

les coordonnées cherchées du point d'où l'on observe \ soient 
M l'angle observé (toujours de la gauche à la droite) entre 
le premier et le deuxième point, et M' l'angle observé entre 
le premier et le troisième (en supposant qu'on en ait, s'il y 
a lieu, retranché 180 degrés); soient 

R _ R' __ , 

sinM-"' sinM'""' 
E — M = N , E' — M' = N': 

on a les deux équations 

p = n sin (« — N) , 
p =ri sin(e — N') : 

qui, écrites de la manière suivante, 

sin (« — N) 

„> 

— gin («_ N')' 

6e résoudront par la méthode exposée ( Theovia Motus 
Corporum cœlestium, page 82). 

L'une des solutions exposées en cet endroit conduit à la 
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règle suivante. Soit n' plus grand, ou au moins pas plus 
petit que n , ce qui est évidemment permis, puisqu'on peut 
choisir le second point arbitrairement; posons 



n 

tang^N ' — N ) 
tang ( 45« — Ç ) 



=: tang : 



on aura 



i 



6 = _( N '4-N)H-4s 



s étant connu, Tune des équations, ou même encore toutes 
deux, fourniront la valeur de p. 

Dans notre exemple , si nous considérons Frauen- 
thurm comme le premier point, Friedriehsberg comme le 
deuxième et Friedrichsthurm comme le troisième, nous 
aurons : 



a = 


710,0, 


b — 


684,2, 


E = 


77 0 19' 3i ",92. 


E' = 


3o8°5i / 45",77, 


logR = 


3,8o,44 2 °5> 


logR' = 


3,5733549, 


M = 


99 0 22' 5o",20, 


M' = 


ioi°n'5o",8o, 


N = 


337 0 56 , 42' , , 7 2 ) 


N' = 


207039. 54", 97, 




3,9002650, 


log n' — 


3 ,5817019. 



n' étant plus grand que w, nous changerons ici l'ordre et 
nous poserons 

N= 207°3 9 / 54 , \97, N' = 33 7 °56'42",72, 
log n = 3,5817019, log n' = 3,9002650; 

on en déduit 

ç= «9039' 3",8 7 , 
\ — 80. 45. 3 .",69, 
e = 353 . 33 . 5o" , 53 , 
log p =3,3303990; 

et les coordonnées du Holkensbastion , 

2836, 44 ! > 444>33o. 
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NOTE V. 

SUR LA DÉTERMINATION CHRONOMÉTRIQUE DES LONGITUDES. 
( Astronomische Nachrichten, tome V, poge 227.) 



Soient 0, Ô', 0", etc., les époques (au nombre de w), 
auxquelles un chronomètre a accusé les différences a, a', 
a", etc., avec les temps des lieux dont les longitudes sont 
x, x f , x", etc., 0, 0', 0", etc., étant supposées réduites au 
temps d'un seul et même lieu et u désignant l'avance jour- 
nalière du chronomètre ; on aurait, si l'instrument était 
parfaitement régulier, les équations 

a — Qu — x = a' — 9'« — x' = a" — 9" u — x" — 

Pour que ces équations suffisent à la détermination des in- 
connues x, x ( , x")..., m, il faut, d'une part, considérer l'une 
des longitudes comme donnée, et, d'autre part, il est néces- 
saire que Ton ait observé au moins deux fois dans le même 
lieu; de telle sorte que deux au moins des inconnues x , x\ 
x", etc., soient égales entre elles. Si parmi ces quantités il 
n'y en a que deux qui soient identiques , le problème est 
entièrement déterminé; dans le cas contraire, il devient 
indéterminé, et l'on doit faire en sorte que les équations 

O = a — a' -+- ( 9' — 9) u — x -h x', 
o = a' — a"-4-(ô"— 9') m — x'-f-x", 
o = a"— a'" -h ( 0 m — 9") * — x"+ x m , 

soient satisfaites aussi exactement que possible, car les 
imperfections inévitables du chronomètre ne permettront 
jamais de satisfaire rigoureusement à toutes. Mais on ne 
doit pas accorder à ces équations des poids égaux; car les 
quantités 

a — a' -f- (9'— 9) u — x -h x', 
a '—a" + (Q"— 9')« — x'+x", 
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représentent l'accumulation de tous les écarts de marche des 
chronomètres dans les intervalles B r — 0, 0" — etc., et 
s'il s'agit d'un bon chronomètre, à qui l'on puisse réelle- 
ment attribuer une marche moyenne sans variation crois- 
sante dans un seul sens, la valeur moyenne à craindre pour 
une pareille somme peut être considérée comme propor- 
tionnelle à la racine carrée du temps écoulé. 

On devra donc, dans l'application de la méthode des 
moindres carrés, regarder les équations précédentes comme 
ayant des poids inversement proportionnels aux différences 
0"_0', Q»>— 0% etc. 

La solution n'a alors aucune difficulté, et fournira les 
valeurs les plus probables de .r, x', x", etc., ainsi que le 
poids de chaque détermination. 

J'ajouterai pourtant quelques remarques. 

1. Si la première et la dernière observation ont été faites 
au môme lieu, la valeur la plus plausible de u est celle qui 
résulte simplement de Ja comparaison de ces observations 
extrêmes. Le calcul devient alors très-simple, car, en vertu 
d'un théorème bien facile à démontrer, on peut, dans les 
équations, remplacer u par sa valeur la plus plausible, 
ou, ce qui revient au même, on peut employer celte valeur 
supposée exacte pour corriger les observations et les rame- 
ner à ce qu'elles seraient avec un chronomètre fictif dont 
l'avance serait nulle. 

IL Si l'on a attribué simplement aux diverses équations 
des poids égaux à 

i i i 

l'unité de précision pour les poids obtenus sera l'exactitude 
de celle que l'on obtiendrait a l'aide du même chronomètre 
observé deux fois seulement , et à un jour d'intervalle ; mais 
pour pouvoir comparer les résultats obtenus à l'aide de 
divers chronomètres inégalement précis, il faut encore in- 
o. I i 
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troduire dans le résultat un facteur qui dépende de la plus 
ou de la moins grande perfection de chaque chronomètre 
employé. 

Pour y parvenir, je suppose que les expressions 
a — a' -+- (9' — 9 ) u — x -h x\ 
a' — a" -+- (9" — 9') u — x 1 -h 

.•••••-«•••» 

deviennent respectivement X, X', X", etc., quand on substi- 
tue aux inconnues leurs valeurs les plus plausibles, soit 

F^TÔ + S*=? + F=F + S ' 

si v est le nombre des inconnues et que l'on pose 

/ S 

le facteur spécifique relatif à chaque chronomètre est pro- 
portionnel à-^-ouà *"~*~"' » et Ton peut regarder m 
comme l'écart de la marche moyenne qui est a craindre 
pour la marche d'une journée. 

III. Les règles qui précèdent sont relatives à un chrono- 
mètre dont la marche n'est soumise à aucune irrégularité 
sensible et croissante avec le temps. Si celte hypothèse n'é- 
tait pas permise, on pourrait supposer, lorsque les observa- 
tions n'embrassent qu'une période qui n'est pas excessive- 
ment considérable, une variation proportionnelle au temps 
dans l'avance journalière de l'instrument, en introduisant 
ainsi une inconnue de plus. 

Les équations prendraient alors la forme suivante : 

o = a — a' -h (G' — 9) « + (9" — 9») v — x 4- 

IV. Pour ce qui concerne la résolution des équations 
d'après la méthode des moindres carrés, il n'est peut-être pas 
inutile de rappeler qu'on doit commencer, dans le plus grand 
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nombre des cas , par calculer les valeurs approchées des in- 
connues et appliquer ensuite la méthode à la détermination 
des corrections très-petites que les valeurs doivent subir. 

Il nous a paru utile de rappeler ce conseil général , parce 
que beaucoup de calculateurs ont paru l'oublier, et ont été 
conduits à des calculs plus laborieux et peut-être moins 
exacts . 

J'ai déterminé la marche des cinq chronomètres suivants : 



Greenwich, 2û Juin 3 . ii 
25 Juill. 1. 15 
28 » 3i3 
2 Août 1 . ih. 
[2 » io. 18 
2Û » 7.37 



in Sept. 7 . /jo 



Hclt;oland, i Juill. 3 \o 

22 » I 3 . /|Q 



5 Août iJfc 
U » liL 9_ 
3_o_ » 19.30 
Q Sept. ÎL fi 
j » 8.4a 



Alloua 



fi Août 5 . 53 



1 '.33 



9_&2 



Brëiue , 



/J Sept. 22^L2 

i_l Août û—2 



s '7 ' < 



2^ 



11 .28,48 



13.5q,4o 



'3-47,98 



1 5. 3.4 ,47 



4o. 8,00 



42. 3,0 2 



43i8,n 



45.53,08 



4G.5i.5G 
46.38,72 



- 5 i.38, 95 
5i .57,35 



:V| . io,H3 
54-39, 1 G 



— 47 .50,65 



?.'-»7 



1 .32, 1 5 
j .36, 96 



2_ G , 2 4 
2.i5,84 



a.4o, 3fi 
30.26,84 



3o. 3,89 
•J9 .43,^3 



39-33,43 
29 7>qG 
28.38,94 



28.56,71 
37.55,76 



37.5o,o3 
3j . 2 1 , 3o 

3y. l5,31 

— 33. 16,49 



nu toi ET 

3o5G. 



„ I u 

4-30.59, 



7_3 4-4^ '■)■?<> 



3o_i 50,07 



3o.3i ,78 
39-35,69 



-';>■ 10,48 



— o. 2o,34 



1 . 10,24 
i .32,75 
3.40,67 



3_i 4>53 



— 9.28,50 



9.38,8i 
io.56,68 



il . i5,36 



BARRA CD 

904. 



4»-4o,a4 



48.58 ? 8 7 



4o- f »7,83 



3^L J7, I_3 



4-16.47,39 



'7 37,5i 
liL 1 ,3o 
'9- ï 7><>3 
'9-43, 80 



9-38,48 



9. 4o>3o 
1 i . 5,93 

ll2 4>49 



1/1 



,86 



KESSEL 
1352. 



zb5ji.32,ii 



50.39,69 
5o.53, 14 
5i.38,66 
3a . 3,43 



4-18.48,39 



'9-^>77 

'9-47 »^ 
20.47,68 

ai. 6,56 
4-11 . 16,25 



12.25,96 
13. 48, LQ 

-liL 5,83 



Prenons pour exemple le chronomètre de Breguet 3o56. 
Soient o la longitude de Helgoland, — x cellede Greenwich, 
y celle d'Allona. Je ne fais pas entrer celle de Brème dans 
le calcul, car n'ayant pour cette ville qu'une seule obser- 
vation , il est impossible de la contrôler. Je compte le temps 
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à partir de la première comparaison du chronomètre n° i 
(Grcenwich, 3o juin, 3 h 22 m ). En substituant au chrono- 
mètre de Rreguet un instrument fictif dont l'avance soit 
nulle, nous trouvons : 

6 



22,4 


-+■ 00,20; 


20 ,0 


+ ! 949. 6 ° — x y 


20 ,0 


-h »t)30 ,07 — .T , 


32,9 


-h 1950,29 — x; 


35 ,9 


-h 59,08; 


37,1 


- 434,98 -h /; 


40,4 


— 433,49 -h;; 


4**4 


-h 5 9 ,88; 


48,3 


-+- 1949,60 — x ; 


56,2 


+ 1952,74 — x; 


6i,6 


-h 6i,32; 


62,2 


— 432,53 -h y i 


66,8 


— 434>9 8 + .r; 


68,0 


H- 60,19. 



Dans les équations ci-dessus, les inconnues x et y sont 
séparées, ce qui facilite leur détermination ; nous trouvons 
pour x quatre déterminations : 





POIDS. 


II 

1889,40 


2 ,6 = °' 38; 


1891,21 


3o = o,33; 


1889,78 




1891,42 


3,4 = °' ,9; 



d'où I on déduit : 

x = i89o",36 i ,07 ; 

on trouve de même : 

y= 494",. a 3,83. 
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D'après ces valeurs, le chronomètre fictif marquera, en 
temps de Helgoland , 



22,4 
25,0 
28,0 

32,9 
35,9 

4o,4 
42,4 
48,3 
56,2 
61,6 
62,2 
66,8 
68,0 



60,20; 

59,24 — 0^96; 
60, 5i + 1 ,27 ; 
59,93 — o,58; 
59 , 08 — o , 85 ; 
59, l4 -h 0 , 06; 

60, 63 -4- 1 ,47 ; 
59,88 — 0,75; 
59,24 — 0,62; 
62,38 + 3, i4; 
61 ,32 — 1 ,06; 
61,59 + 0,27; 
59,14 — 2,45; 
60,19+ 1 ,o5; 



d'où l'on déduit : 



S = 6,00, 



m 



v 



■ 



et Terreur moyenne à craindre est : 

pour x, o",75; 

pour x » o",4°- 

•■ 

Les résultats fournis par les cinq chronomètres donnent : 





x = 1890,36 


ERREUR MOYENNE 

à craindre. 


POIDS. 


Breguet 


0,75 


1,78 


Kessel 


1893,39 


0,67 


2,23 


Barraud .... 


1892 ,32 


°»49 


4,l6 




1892,39 


0,43 


5,4l 


4 


1 892 , 52 


o,35 


8,l6 


Moyenne. . 


x= 1892,35 




21,74 



Digitized bty Google 



( K» ) 

on trouve de même, pour y : 

Breguet y = 494» 12 °»'° 6,25 

Kessel. tïfîify o,36 7,72 

Barraud. ... 49^67 0,21 >4>79 

' 49 3 »9 8 °' 2 9 1 

4- 494» , 6 °» a 4 » 7 , 36 

58, 01 

Le nombre placé sous le nom de poids dans la dernière 
colonne, est l'inverse du carré de Terreur moyenne à 
craindre, en prenant ainsi , pour poids unité, celui qui cor- 
respond aux observations qui donnent une erreur moyenne 

à craindre égale à 1", de sorte que, pour Alloua, Terreur 

," 

moyenne à craindre est - = o", i3} mais il vaut mieux 

v58,oi 

considérer les nombres de la dernière colonne comme indi- 
quant seulement des rapports , et déduire la précision ab- 
solue de la différence entre les valeurs des derniers résultats 
trouvés pour x et y au moyen de chaque chronomètre. La 
précision trouvée de cette manière sera toujours un peu trop 
grande, puisque les déterminations de temps à Greenwich, 
à Helgoland et à Altona n'ont pas une précision absolue, 
et que, par conséquent, quel que soit le nombre des chro- 
nomètres, les erreurs qui proviennent de cette source se 
feront toujours sentir dans chaque résultat final. 

On peut obtenir, de la manière suivante, la longitude 
de Brème. 

Soit z cette longitude à Test de Helgoland -, la comparai- 
son du chronomètre de Breguet donne la position du chro- 
nomètre fictif, 

— i65",52 + 2, 

et Ton déduit de la comparaison avec les résultats précé- 
dents, 

POIDS. 

z = 325*40 — 0,7 i 

1 >H 
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( '67 ) 

les autres donneraient 



POIDS. 

2=234,76 -ç-j; = 0,2 

225,24 0,9 

Le poids 0,9 doit être multiplié par g-j^i les cinq ohrono 
mètres donnent : 

Breguet 226,24 i,5 

Kessel 225,84 1 >9 

Barraud 225,3g 3,6 

1 226,04 2 >9 

4 224,86 4,3 

225,42 14,2 



La longitude de Brème, qui serait, d'après cela, 268", 54 
à l'ouest d'Altona , est naturellement affectée des erreurs 
sur la détermination du temps à Brème , et cette différence 
semble trop petite de plusieurs secondes. D'après mes trian- 
gulations, la tour de Ànsgarius est de 273", 5 1 en temps à 
l'ouest de Gôttingue, et l'observatoire d'Olbers à 271", 9. 
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